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Summary (*). — Two experiments are reported, one performed at Cervinia 
(3500 m a.s.l.) and one performed at Rome, on the zenithal distribution 
of penetrating particles and mesons. It is found that the penetrating com- 
ponent at sea level can be represented by the law /(0) = I, cos” 9 with 
n = 2.14 for E = 500 MeV whereas at 3500 m a.s.l. the distribution 
does not obey this law. At sea level, moreover, the mesonic component 
with energy between 0.5 an 1.9 GeV follows very closely the law cos*.°% 9. 


(*) Editor’s translation. 


1. — Introduction. 


Several measurements of the intensity of the penetrating component at 
sea level have been made in the past. One of the most complete and accurate 
determinations of zenithal dependence of cosmic radiation at several altitudes 
was made by GREISEN (!) who found that the zenithal dependence could be 
expressed by the equation 


UO) == Tg Cos” 0! 


(4) K. GREISEN: Phys. Rev., 61, 212 (1942). 
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where J, is the vertical intensity, 0 the zenith angle and n is a constant deter- 
mined from the experiment; for « penetrating component» (300 MeV for 
mesons) n was found to be 2.1 while for «total radiation » Greisen found 
n = 2.3 at sea level. 

Notwithstanding these accurate measurements some doubts remain on this 
question. In the past everybody has assumed the penetrating component of 
cosmic rays to be completely composed of u-mesons; at the same time 
several authors defined penetrating particles as the particles with a range 
greater than 10 or 15 em of lead while we know that several electrons have 
penetrations greater than 10 em of lead. 

Probably the greatest source of error is in the assumption that all pene- 
trating particles are mesons. Many authors (?) have shown that at high alti- 
tude protons are present together with mesons, although the intensity of pro- 
tons at sea level, and at mountain altitudes is not yet known. 

Recently KRAUSHAAR (*) and ZAR (*) found that zenithal spectrum depends 
on the type of technique used for the measurement; the exponent n was found 
by KRAUSHAAR with a coincidence technique to be 2.09 for mesons of mo- 
mentum greater than 200 MeV/e or for other particles of the same range; 
with an anticoincidence technique he found n equal to 3.3 for particles with 
momentum between 200 and 330 MeV/e and he found the same value with 
delayed coincidences selecting only mesons with momentum between 25 and 
170 MeV/c. ZAR with similar technique found that n is of the order of 2 in 
coincidence measurements and of the order of 3 in delayed coincidence or 
anticoincidence measurements. He found that his results are consistent 
with the hypothesis that the exponent for zenith angle distribution decreases 
as the meson energy increases, but statistical errors and also the fact that only 
two angles were investigated make a satisfactory conclusion impossible. 

Moreover VoISIN (5) with a delayed coincidence experiment found an expo- 
nent n— 3.0 for mesons of momentum between 300 an 510 MeV/c at 30° zenithal 
angle; on increasing the zenithal angle, however, n decreases (2.26 + 0.07 at 60°, 
1.71 40.06 at 80° for mesons of momentum between 300 and 410 and 
2.59 + 0.08 at 60°, 2.11 + 0.06 at 80° between 410 and 510 MeV/c). Accord- 
ing to Voisin’s experiment zenithal distribution can be expressed by equation: 


1(0) = 1,(1 — a sin? 6), 


(2) M. Conversi: Phys. Rev., 76, 444, 849, 851 (1949); T. C. MErKLE, E. L. 
GOLDWASSER and R. B. BRODE: Phys. Rev., 79, 926 (1950); E. L. GOLDWASSER and 
T. C. MERKLE: Phys. Rev., 83, 43 (1951). 

(5) W. L. KRAUSHAAR: Phys. Rev., 76, 1045 (1949). 

(4) J. L. Zar: Phys. Rev., 83, 761 (1951). 

(5) A. G. VorsiN: Phys. Rev., 84, 850 (1951). 
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where a~1 and b=1.47 + 0.12 in the momentum range 300-410 MeV/e 
and b = 1.61 + 0.15 in the momentum range 410-510 MeV /e. 

In addition it has been found recently in photographic plate work (°) that 
the zenithal distribution of mesons of energy greater than 40 MeV at 2000 m 
above sea level can be roughly expressed by the law: cos®:50. 

For the reasons previously explained we think it is of some interest to com- 
pare the results of our measurements performed at « Laboratorio della Testa 
Grigia» (L.T.G., Cervinia, 3500 m 
above sea level, 400 50’ North geogra- 
phical latitude) in 1951, with those 
performed in Rome in 1952. 


2. — Apparatus. 


The apparatus used at L.T.G. con- 
sisted of two winding magnetic cir- 
cuits, one placed over the other, in 
which by means of a direct current, it 
was possible to create a magnetic field. 
A telescope of three 450 em? coun- 
ters A, B, C defines a solid angle 
through the two magnetic circuits as 
shown in fig. 1. Anticoincidence coun- 
ters prevent spurious counts due to 
lateral showers. The apparatus could 
rotate about an horizontal axis. This 
apparatus has been used for positive 
excess measurements (7) but we shall 
now refer to data collected for the total component (positive + negative) in the 
vertical direction, at 45°, 60° both in the North-South and in the Bast-West 
planes. 

The total thickness of iron was 40 em which corresponds to about 500 MeV 
for mesons, but three-fold coincidences (ABC) recorded also other penetrating 
particles such as protons; we will refer to the total radiation penetrating 40 em 
of iron as «T ». 

With the experiment in Rome the geometry of the counters was changed and 


See 2e 6 
&) © 2 SE Se En 


(6) G. Baroni, G. CORTINI, A. MILONE, L. SCARSI and G. VANDERHAEGHE: Nuovo 
Cimento, 9, 867 (1952). 
(7) I. F. Quercıa, B. RıspoLi and 8. Scruri: Nuovo Cimento, 7, 899 (1950). 
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in addition a set of counters D was 
placed under 20 em lead (fig. 2). The 
« central counters» A,A,A,; B,B,B,; 
C,C.C, were connected to a three-fold 
coincidence circuit «7,» and a four- 
fold coincidence « T,» was formed with 
counters D. An electronic device was so 
arranged that no counts were recorded 
if more than one counter of each set 
A, B, C was discharged. The cut off 
energy for «7,» results of 500 MeV 
and for «7,» of 800 MeV for me- 
sons. Moreover three-fold coincidences 
(A,B,C;); (A, B.C2); (A43B3C3); (A,B,C,) 
are mixed together and recorded as 
«M,». «M,» in coincidence with the 
set of D counters (?) is recorded as 
€ M,». It should be noted that M-coin- 
cidences refer only to mesons because 
protons of the same range have a mi- 
nimum energy of 1.8 GeV, so that the 
magnetic field is inefficient to curve their 
paths. For this reason we refer to these 
counts as « M », or meson-component. 


3. — Experimental Results. 


In table I we report the ratios 
T(0)/T(0°) as function of zenithal angle 
in East-West plane at L.T.G. and at 
Rome (?). These data have been obtained 
by averaging the results obtained in 
the eastern and the western directions. 


(5) A more detailed description of ap- 
paratus will be given elsewhere. 

(9) It must be noted that the ratio 
T(6)/T(0°) is a quite good approximation 
for the ratio between the corresponding 
absolute intensities. The error involved is 
of the order of 0.4%. 
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TABLE I. — Relative intensity of total penetrating component in East-West plane (T(6)/T,;). 


votes Energy 5 
een (MeV) range 150 300 450 | 60° | 
m a.s.l. f ER | | 
or mesons | | 

«Ts» >500 0.915 40.008 0.716 -0.007/0.502 -0.004.0.227 + 0.002 

50 «T,» = 800  0.892-L0.008 0.708 £0.007 0.510 £0.003 0.241 +-0.002 

« T,—T, » 500 < E < 800 1.050 10.020 0.771 £0.015 0.423--0.009.0.133+0.003 

3 500 T,» > 500 = = 0.541 0.003 0.209 +0.001 


Data of table I and fig. 3 show that integral zenithal spectra at sea level 
for particles which have range greater than 40 em Fe (500 MeV for mesons) 
and for particles whose range is greater than 40 cm Fe + 20 em Pb (800 MeV 
for mesons) are in a very 
good agreement with the Fr 


+ 


law: 


T(0) = I, cos" 0, 


where 7 is found to be 2.14 
for the former case and 2.09 
for the latter case. 

On the contrary the dif- 
ferential spectrum of part- 
icles whose range is equal 
to that of mesons lying bet- 


ween 500 and 800 MeV at oat 3 25 

the sea level and the integral Fig. 3. 
spectrum at 3500 m are not AU at 5. 1. 

in agreement with a cos” + T,— T, at sl. 

law, as one can see from © T; at 3500 m a.s.l. 


fig. 3 where as abscissa we 
have plotted log sec 4 and as ordinate log 7(9)/T(0). Dotted lines represent 
cos? # and cos® 4. We can say that # is a function of zenithal angle and the 
average value is 2.40 for integral spectrum at 3500 m and 3.06 for differential 
component in Rome. 

In table IT and in fig. 4 we report the ratios M(#)/M(0°) which refer only 
to mesons recorded in the East-West plane at sea level. 

From fig. 4 one can see that the exponent #7 seems to decrease with in- 
creasing zenithal angle with an average value 2.55 (2.61 and 2.49 respectively 
for the two ranges of energy). 
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Tage II. — Intensity of meson component in Bast-West plane at sea level (M(6)/M,). 


à 

| | Energy range 15° 300 | 459 | 600 
| 

| | 


| | | 
«My» | 500 <E <1900 | 0.77 +0.02 | 0.58 0.02 | 0.365-+0.012 | 0.163 40.006 | 
«My» | 800 <E <1900 | 0.7404-0.031 | 0.7004-0.030 0.377 0.017 | 0.184 +0.008 


In table III we report the results obtained in North-South plane; unfortu- 
nately the data collected are 
not so complete as the pre- 
ceding ones, but they do not 
show any significant diffe- 


rence. 


Present results are in 


very good agreement with 
measurements of ZAR and 
KRAUSHAAR for the T com- 
ponent. For M component 


we have found an exponent 
n lower than theirs but if 
we interpret most of the 7 
component stopped in 20 em 
of lead as meson component 
we have for mesons with 
energy between 500 and 800 MeV, n = 3.06 and for mesons with energy between 
500 and 1900 MeV, 2.61 and between 800 and 1900 MeV, 2.49. From these 
data we can conclude that n decreases when energy of mesons increases and 
so Kraushaar’s value 3.3 for mesons with momentum between 25 and 170 MeV 
is not in disagreement with the present work. In a next work we shall at- 
tempt an interpretation of present data. 


TABLE III. Intensity of T and 


: : = = = 
Altitude | Energy (MeV) a | } 
m a.s.l. | range for mesons | Bp | oP wy 
| > 500 | 0.745 40.007 | 0.522 40.002 | zu 
T(6)/T, 50 > 800 | 0.736 40.006 | 0.512 40.003 = 
500 <E<800 | 0.801 +0.015 | 0.477 +0.007 — 
3 500 > 500 | — 0.551 +0.004 | 0.228 +0.002 
50 500 < H <1900 | 0.641 +0.030 | 0.374 +0.010 | = 
| | 800 <E < 1900 | 0.602 +0.030 | 0.353 40.010 = 
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RIASSUNTO 


Vengono riportati due esperimenti, uno eseguito a Cervinia (3500 m s.l.m.) ed uno 
a Roma sulla distribuzione zenitale delle particelle penetranti e dei mesoni. Si trova 
che la penetrante al livello del mare pud essere rappresentata dalla legge I(6) = 1, cos" 6 
con n = 2,14 per H= 500 MeV, mentre a 3500 m la distribuzione non segue tale legge. 
Al livello del mare, inoltre, la componente mesonica con energia compresa fra 0,5 e 
1,9 GeV ha un andamento assai prossimo a cos?-60. 
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Riassunto. Si calcola il contributo alla componente elettronica della 
radiazione cosmica generato, lungo l’atmosfera, dal decadimento del 
mesone 7° in due fotoni, usando le equazioni di diffusione con una sorgente 
di fotoni il eui andamento in funzione dell’energia & stato desunto dai 
dati sperimentali del gruppo di Bristol, e la cui variazione con la quota 
© regolata da un cammino di assorbimento di 130 g-em-?. L’integrazione 
delle equazioni di diffusione si fa con un metodo di approssimazioni suc- 
cessive tenendo conto sia della perdita di energia per ionizzazione sia 
della variazione delle sezioni d’urto di produzione di coppie e di brems- 
strahlung con Venergia. Si normalizza il risultato sullo spettro integrale 
della componente elettronica totale quale risulta dai più recenti risultati 
sperimentali e si determina il rapporto dei mesoni neutri a quelli carichi 
ottenendo r’/n= — 0,53. Dal confronto dei valori assoluti cosi ottenuti 
con quelli forniti dal gruppo di Bristol, si conclude con la compatibilita 
dei due dati entro i limiti delle notevoli incertezze relative alla cono- 
seenza di aleuni processi necessaria a detto confronto. 


1. — Introduzione. 


E noto che per rendere conto della componente elettronica totale della 
radiazione cosmica nell’atmosfera, occorre postulare la esistenza, oltre a quella 
degli elettroni di knock-on e di decadimento delle particelle u, anche di un’altra 
sorgente (1) che, con la scoperta della disintegrazione del mesone 7°, & stata 
individuata nello spettro dei fotoni che da esso derivano. 


(1) B. Rossi: Phys. Rev., 61, 121 (1942). 
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Calcoli per verificare tale ipotesi sono stati fatti da BUDINI (?) e da CAL- 
DIROLA (?). Entrambi questi Autori ottengono la sorgente di fotoni da consi- 
derazioni di ordine teorico sul meccanismo di produzione di mesoni neutri 
da parte di nucleoni primari, e nel calcolo della cascata elettrofotonica che deriva 
dai fotoni di disintegrazione dei r° introducono delle semplificazioni alquanto 
grossolane per le basse energie, dove adottano lo schema di ARLEY (*) che tien 
conto delle sole perdite di energia per jonizzazione ed assume per tutte le 
sezioni d’urto la loro forma asintotica, o le traseura per valori dell’energia 
minori di quella critica. 

Col presente Javoro intendiamo riprendere tale caleolo considerando anzi- 
tutto la cascata elettrofotonica che nasce da una sorgente di fotoni di cui 
abbiamo desunto l’andamento dai dati sperimentali pit attendibili; si evitano 
cosi tutte le jpotesi sui meccanismi di collisione dei nucleoni; inoltre si cerca 
di tenere conto in modo più accurato di tutti quei fattori che influenzano l’evol- 
versi della cascata, sia non trascurando mai le perdite di energia per ioniz- 
zazione, sia prendendo in considerazione la variazione delle varie sezioni d’urto 
con Venergia. 


2. — Sorgente di fotoni e sezioni d’urto. 
Precisiamo qui i dati del nostro calcolo. Assumiamo per la sorgente di 
fotoni la forma 


(1) I(W,t) = NS(W) exp [At] , 


dove N & una costante di normalizzazione, W Venergia dei fotoni misurata 
in GeV, S(W) lo spettro energetico che è stato desunto dai dati sperimentali 
ottenuti mediante lastre nucleari dal gruppo di Bristol (5) e cosi approssimato: 


W-2s per W>1 GeV, 

ws 02<W<1, 
S(W) = 

4,26 W- 0,075 < W <0,2, 


1,65-105 W25 0,01<W<0,075. 


L’andamento di S(W) per W>1 GeV, che non & dato nel lavoro sopra 
eitato, & stato uniformato su quello della sorgente dei mesoni 7 carichi che 


(?) P. Bupını: Nuovo Cimento, 7, 868 (1950). 
(3) P. CALDIROLA e P. GULMANELLI: Nuovo Oimento, 8, 229 (195]). 
- (4) N. ARLEY: Proc. Phys. Soc., A 168, 519 (1938). 
(5) A. G. CARLSON, I. E. Hooper e D. T. Kina: Phil. Mag., 41, 701 (1950). 
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deriva dai dati di SANDS corretti alle alte energie in base all’andamento dello 
spettro delle particelle  sottoterra. 
Infine l’andamento con la quota, rappresentato dall’esponenziale nella (1), 
è stato scelto assumendo un cammino di assorbimento di 130 g-cm~? pari a 
quello della componente nucleonica di media energia, tra 1 e 10 GeV (§), 
(A, = — 0,33); t è misurato in unita di radiazione che in aria risulta di 43 g-em°?. 
Per comodita, nei calcoli che seguiranno abbiamo posto: 


[ S(W)=S(W) per W>1GeV, 
| Sı(W) = SEW) per W>02, 
= | SW) = S(W) per W>0,075, 

| S.(W) = S(W) per W>0,01. 


Le varie sezioni d’urto per unita di radiazione, per le quali abbiamo usato 
le notazioni di ROSSI e GREISEN (7), sono state approssimate nel seguente modo: 
sezione d’urto differenziale di bremsstrahlung 


{ (Ww\aw U wiaw 
== » | = a 
(7) 7 " w h- | a per U>0,2 GeV, 


W\ dw U U\| dW 
ue. | is = | > a 2 
| AL a a hr h- (7) v per U< 0,2, 


dove U & Venergia totale dell’elettrone misurata in GeV ed 


aw 
CNT Men 


f (2) —0,2 mn a — 0,6 
le 


€ una funzione che abbiamo introdotto per tenere conto della variazione della 
sezione d’urto con l’energia dell’elettrone; 
sezione d’urto differenziale di produzione di coppie 


2 T 
Wo PE — Zus h La + 7 ay per W> 0,2 GeV, 


w/w WwW Ww Ww 
> OO) 
Be rai au 2W  109\aU 
LM = vf In the er) PetignW—2.0;275 


in entrambi i casi b = 0,012 ed h = 4/3 + 2b = 1,357. 


() J. G. Witson: Progress in Cosmic Ray Physics (Amsterdam, 1952), p. 317. 
(7?) B. Rossı e K. GREISEN: Rev. Mod. Phys., 13, 240 (1941). 
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Abbiamo assunto costante.la perdita di 
energia per ionizzazione ed uguale a 


p = 0,1 GeV per unita di radiazione . 


Per le sezioni d’urto totali di produzione 


ws 
di coppie e Compton, infine, dato che il mi- SS ae u 
nimo nel caso dell’aria non & molto pro- Ba Er 
fondo, abbiamo posto Fig. 1. — Confronto fra la espres- 
sione approssimata (- — —) della se- 
zione d’urto differenziale di brems- 
On + Ge = Go = 0,773 . strahlung per unita di radiazione 


e l’espressione esatta (——). 


Inoltre, poiché nel nostro calcolo avremo 


sempre a che fare con energie > 0,01 GeV, abbiamo trascurato l’energia 


ET - 
in o (W) 
1.2 

3 
af 1.04 


dal 
“| 7 
x 
a 
Se 0 = — — — 
° ner a SE el 3 ONG 10°? 102 vo 1 WGw 
Fig. 2. — Confronto fra la espressione ap- Fig. 3. — Sezioni d’urto totale di pro- 
prossimata (-—-) della sezione d’urto diffe- duzione di coppie e Compton per 
renziale di produzione di coppie per unitä unita di radiazione. 


di radiazione e l’espressione esatta (——). 


u, di riposo dell’elettrone rispetto a quella cinetica Æ e cio@ abbiamo posto 
Ox E. 

Il confronto fra le effettive sezioni d’urto e quelle semplificate che abbiamo 
adottato risulta dalle fig. 1, 2, 3. 


3. — Equazioni di diffusione. 


. Indicando con x(E,t)dEdt e y(W, t)dW dt rispettivamente gli spettri dif- 
ferenziali degli elettroni e dei fotoni, ed attenendoci alle precedenti notazioni 
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ed approssimazioni, le equazioni di diffusione si scrivono: 


co E 


onl Et Nip dw aw 
(4) m 2 | WW, t)y(W, B) ra an hc 
J : 
= dE' on(E, t 
a2 | n(E', t)g(£", E'—E) E + Di ? 
3 
an) a | a(E, Op, W) ee o,y(W,t) + ZW, t), 


con le eondizioni iniziali: 


(5) n(B,0)=—0, y(W,0)=0. 

Abbiamo integrato questo sistema per approssimazioni successive e preci- 
samente abbiamo assunto per approssimazione zero e per E, W >1 GeV, la 
soluzione della (4) con 5 = 0, che rappresenta il caso in cui le perdite di 
energia per ionizzazione sono trascurabili; gli spettri che cosi si ottengono 
sono gli stessi dati da BUDINI (?) per energie maggiori di quella critica degli 
elettroni in aria. Abbiamo ricavato l’approssimazione successiva considerando 
il sistema (4) per energie H, W > 0,2 GeV e risolvendo prima l’equazione in y, 
ponendo al posto di x la approssimazione zero 7; la prima approssimazione y; 
cosi ottenuta si è quindi posta nell’equazione in z al posto della y. In modo 
analogo abbiamo proceduto per le successive approssimazioni: la seconda per 
HL, W >0,075 GeV e la terza per E, W > 0,01 GeV, come si vedra meglio 
nella descrizione pit dettagliata seguente. 


Approssimazione zero (valida per BE, W >1 GeV). 
Ponendo nella (4) 6 = 0 si ottengono soluzioni del tipo 
{ x(H, t) = Np(t)S,(B) , 


(6) 
| yo(W, t) = Ng(t)S,(W) , 


con le funzioni p(t) e g(t) che devono soddisfare al sistema 


dp(t) Br, 
mp + Bel), 
| VO _ opt) g(t) + exp [At] , 
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la cui integrazione, con le condizioni iniziali 
ALU HOUR 


direttamente ricavate dalle (5) e (6), si fa nel modo consueto ottenendo 


2 


p(t) = Dıpı exp [Ai], 


0 


git) = > g; exp [Ait], 


ARE a A Eu Bs eus 
a N N om PRES AT ay) Ua T7 m nf 


= te ca A mens rt a 
BA A) TG Da) #7 


SRH. 


con le lettere A, B, C, À, A, abbiamo indicato i valori che assumono le fun- 
zioni A(s), B(s), C(s), A,(s), A,(s) riportate nel citato lavoro di Rossı e 
GREISEN per s = 1,8. 


Prima approssimazione (valida per H, W > 0,2 GeV). 


Seguendo il metodo sopra descritto abbiamo ricavato gli spettri differen- 


ziali x, e y, risolvendo prima l’equazione in y, in cui per la x si & posto la a 
di cui abbiamo esteso la validità anche per # < 1 GeV, e scrivendo, per la (1'), 


»(W, t) = NS,(W) exp [ht] . 
Si ha allora 


5 a 
N _ WOp(t)8(W)— aoys(W, t) + NSW) exp [A], 


il eui integrale generale, con la condizione iniziale y,(W, 0) — 0, & dato dalla 


(7) yi(W, t) = N[kd)S,(W) + m(t)A,(W)] , 
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ove 
kit) = Eee Lot] = An or 2 
Co + Ay 
Be exp [Ait) — exp [— cit] 
mh Ges 2: Re On À à 
Re Au, 1) A7 de Y17-2,8 
H,(W) = | Rpt) PU M) = OW 


Si vede immediatamente che k(t) + mt) = g(t) e quindi lo spettro yı(W, t) 
si riduce a y,(W,t) per W >1 GeV. Osserviamo che questo spettro fotonico 
tiene conto della variazione della forma della sorgente con l’energia ma non 
della perdita di energia per ionizzazione; questo si fara nelle approssimazioni 
successive. 

Poniamo ora y, al posto di y nella equazione in z e facciamo la posizione (8): 


T (E' t E \a@:50+1 
(8) Ma ’ al | 1 


we, t) I NE 


ove s,(H,t) è una funzione di # e ti cui valori vanno da 1,8 per H > 1 GeV 
e t qualunque, ad un valore compreso fra 0,9 ed 1,8 per 02 <E<1 GeV. 
In mancanza di indicazioni più precise, e che potrebbero essere ottenute solo 
conoscendo gia la soluzione z,(Z,t), abbiamo posto: 


| 0:68 HEA 17 02 


E < 
s1(E;, t) = 
LAS E > 
cioè abbiamo scelto per s, una espressione lineare in Æ ed indipendente da # 
e tale che per H = 0,2 GeV essa abbia un valore medio fra gli estremi 1,8 e 0,9. 
Con la posizione (8) l’equazione in 7, diventa: 


om (Hi, t) 


dw 
(9) re 


+ AB, = 2 | OP, PO, B) 7 


Co 


dove con A,(#) abbiamo indicato la funzione 


E © 

; dw "(Hf \% dH’ 
4,0 = | ecw, [jee 2-2, 
E 


(5) N. DALLAPORTA e E. CLEMENTEL: Nuovo Oimento, 3, 235 (1946). 
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che & la Alsı(P)] di Rossi e GREISEN (7); nei calcoli l’abbiamo approssimata 
con la espressione 


0,32 + 1,16 0,2 <'E 
A,(E) = 
1,46, E 


Per integrare la (9) facciamo il cambiamento di variabili seguente 


[E=n-Pr, 
(10) 
| ¢ = ihe 


per cui essa diventa 


om, (n, T) 


Cha 


+ AB =2| 7 VAT UE Be) , 


n-Br 
con la condizione iniziale 7,(7, 0) = 0, ed il suo integrale generale è allora: 


g 


i : 17 ‚dw 
am) = 2 | ag exp | fau eras | EAU 


: 
+ n—BC 


Tornando, dopo alcuni passaggi, alle notazioni in Æ e t tramite le (10), la 
espressione definitiva dello spettro differenziale di prima approssimazione % 
risulta data da 


E+ßt 


su(E,t) de exp 
“il 


Ponendo in questa al posto di y, la (7), si ottiene una espressione non inte- 
grabile in termini finiti, che abbiamo calcolato numericamente. Il risultato 
(fig. 4) mostra in particolare che 7,(1,t) è inferiore a z,(1,t) di circa 18%, il 
che rappresenta l’effetto a quella energia della perdita di energia per ionizza- 
zione; per E ~ 2 GeV le due approssimazioni 7, e % praticamente coincidono. 

Abbiamo approssimato lo spettro x,, per poterlo usare nelle approssima- 
zioni successive, con la espressione analitica seguente 


aura | f(r B=) roro 


€ 


Es E>2 GeV, 
(11) m(E, t) = Np(t) | 0,8128», 05<H<2, 


(t > 6), 0,707 E-»> 02<E<05, 
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che rispetto ai valori calcolati numericamente e affetta da un errore massimo 
del 7°%, finché t< 6; per valori inferiori di t non & possibile approssimarlo 
con variabili E e t separate senza introdurre degli errori troppo forti. 


Seconda approssimazione (valida per E, W > 0,075 GeV). 


Seguendo il solito metodo, integriamo l’equazione in y del sistema (4) con 
le espressioni (2) e (3) per le sezioni d’urto ponendo al posto di x la (11), ed 
estendendo anche all’intervallo 0,075 < E< 0,2 GeV l’andamento in — 2,3. 

Si ha allora, come dai calcoli della precedente approssimazione, 


(12) 4(W, t) = N[K(t)S.(W) + m(t)H.(W)] , 
con 
(ma, t) 
H,(W) =| aay Me MO a 


Questo spettro y, tiene conto della perdita di energia per ionizzazione in 
quanto dipende da m,; per tale ragione per W > 0,2 GeV esso da dei valori 
minori dello spettro y, nello stesso intervallo. 

Analogamente alla (8) poniamo ora: 


72 (E", t) = er 
72 (E, t) EH É 


dove la s;(E,t), assunta indipendente da t, & tale che per H 20,2 GeV essa 
assume i valori dati dagli esponenti dello spettro 7,, che sono noti, e per 
E = 0,075 GeV un valore medio fra 0 e 0,9. Con questa ipotesi abbiamo cal- 
colato numericamente la funzione 


E foo} 
A,[s,(E, 9] = | 9, W)  — | (=) Y(E', E— BE) 
0 £ 
che è stata approssimata colla forma analitica più semplice 


| 1,46 H>2 GeV, 


tia 
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Allora Vequazione in 7, si scrive 


n,(E, t On,(E,t à 7 
ep RE 4 ire 1) 2 rw. opor, m, 


et à CE | 


E 


e la sua soluzione, ottenuta analogamente a quella del caso precedente, risulta 
E o 
1 7 [LE +Bt— Ww 
| deexp F | Aste) ae || AL u) vw, à TT. 


Ponendo in questa la espressione (12) di y,, abbiamo calcolato numericamente 
(EH, t). Si trova che 7,(0,5,t) ~ 2,(0,5, t) mentre per valori di H< 0,5 GeV 
si ha sempre x, < x (fig. 4). 

Nella approssimazione seguente abbiamo rappresentato la 7,(£,t) con la 
espressione analitica 


E-28 E>2 Ge, 
0,812 E25 05<E<2 

(13) n,(E, t) = Np(t) ! 0,575 E- 02<E<0,5, 
0,416 Es 0, < 0,2, 
0,131 E 13 0,075 <E<0,1, 


che rispetto ai valori calcolati numericamente è affetta da un errore massimo 
del 6%. 


Terza approssimazione (valida per E, W > 0,01 GeV). 


Calcoliamo anche qui dapprima lo spettro fotonico integrando la equa- 
zione in y col porre al posto di x lo spettro 7, nella forma analitica (13) (a dif- 
ferenza dei casi precedenti perd poniamo ora 7,(E < 0,075, t) = costante per 
il fatto che la sorgente in questo intervallo diminuisce ed il seguire l’anda- 
mento della (13) condurrebbe a dei valori in eccesso). 

Abbiamo quindi 


ya(W, t) = N[A(t)S,(W) + m()H:(W)] , 


con 


Ji Nuovo Cimento. 
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Facciamo ancora la posizione 


nt(E', t) ” (= SE, 1) 41 
7;(E, t) cé oO ’ 


e scegliamo per s, una espressione dipen- 
dente solo da Æ e tale che per E == 0,01 
essa assuma il valore 0 e per H > 0,075 i 
valori che si ottengono con la (13). 

La funzione A,(E) risulta poi espressa 
dalla 


A:[s(E, t)] = 


E oO 
aw fe ig oe TE 
Je n° | (qe) 9 BE) > 


e la equazione in zt 


en,(E, t) ém(E, t) 


+ ß SE -"A,(E)n;(E, te 


et 
10° 10" 1 E GeV a = 
dW 
=, W, t)y(W, E) — 
Fig. 4. — Spettro elettronico diffe- rk ae Ww’ 
renziale nelle varie approssimazioni “ 
e spettro energetico della sorgente 
di fotoni. la cui soluzione 
E+ßt E © 
2f iL if 1 ? „E+ßt—e e dw 
Zu) = a de exp B A,(e') de m) y(W, €) ww’ 
E é € 


da lo spettro elettronico finale. Esso & stato tabulato per vari valori dit; si 
trova che 7,(0,2, t) = 2,(0,2, 2). 

Nella fig. 4 abbiamo riportato gli andamenti delle varie approssimazioni 
No, Wy, 7%, %s per aleuni valori di t. 


4. — Spettro integrale e normalizzazione. 
Partendo dai dati di 2,(H,t) abbiamo calcolato numericamente lo spettro 


integrale //(t)/N degli elettroni per valori di Z > 0,01 GeV e per varie quote; 
i risultati sono riportati nella tab. 1. 
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TABELLA I. 
5 : 
| t 6 10 20 24 
| IIt)/N 6.31 2,59 0,130 3,6-10-2 


Possiamo ora ottenere il valore della costante N di normalizzazione con- 
frontando lo spettro integrale //(t)/N, dato dalla precedente tabella, con lo 
spettro integrale sperimentale //(t) degli elettroni di decadimento dei neutretti 
alle varie quote. 

Abbiamo rieavato quest’ultimo nel seguente modo: dalla componente elet- 
tronica, ottenuta per sottrazione dalla componente totale delle componenti 
mesonica e protonica totali secondo i pin recenti dati sperimentali (*), ab- 
biamo tolto: 


a) gli elettroni di decadimento delle particelle u. Questo dato é stato 
desunto dal recente lavoro di STROFFOLINI (**): 


b) gli elettroni di knock-on, dedotti dai dati di JANossy (**), e di Rossi e 
K£APwANx (27) per il valore al livello del mare, e per le quote più alte seguendo 
Vandamento dei mesoni veloci. 

Rimane in tal modo il solo contributo degli elettroni di decadimento dei =°, 
che confrontato coi valori della tab. I permette di determinare la costante di 
normalizzazione N. 

La tab. II riassume l’insieme dei dati 


TABELLa II. 


t Compon. 


à (a) (6) I(t) IT(t) N 
re elettronica = 
ur. | gem? 10 10 10 10 N 102 
10 430 700 110 19.5 570,5 2,59 2,20 
16 690 160 32,5 19,5 117 0,43 (*) 2.72 
20 860 62 18 7,4 36.6 0,13 2,82 
24 | 1033 24.5 10 55 9 0,036 2.50 


Il valore (*) & stato ottenuto per interpolazione grafica dai dati della tab. L 


) G. Puri: in corso di pubblicazione sul Suppl. al Nuovo Cimento. 
(5) R. SteorroLısı: Nuoro Cimento, 10, 300 (1953). 

(4) L. Jasossy: Cosmic Rays (Oxford, 1948), p. 243. 

(2) B. Rossı e S. Kiapwan: Phys. Rer., 61, 414 (1942). 
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Come si vede i valori ottenuti per N, mediante gli spettri sperimentali alle 
varie quote, risultano abbastanza concordi con un valore medio di 


N = 2,6-10-. 


5. — Rapporto 2°/z=. 


La normalizzazione dello spettro calcolata nel modo precedente ci con- 
sente di ottenere una determinazione sperimentale del rapporto x°/x>. 

Consideriamo lo spettro di decadimento dei 7°, sorgente dei fotoni, inte- 
grato su tutta l’atmosfera: 


24 
(14) 6,(W)dW = aw | 2(W, t)dt= 7,73-10-2 W-28 dW, (W>1 GeV). 


0 
Se ora assumiamo il seguente schema di decadimento del mesone neutro: 
Re 


e ammettiamo, in prima approssimazione, che l’energia dei fotoni generati sia 
quella media, cioè meta di quella del mesone neutro, lo spettro dei x° risulta: 


"an. 


“a 


(15) GW) AW = als 
cioè, per la (14), 


Ce AVON = DL Siac IS ed (W >2 GeV). 
Lo spettro dei mesoni carichi a sua volta si puö ottenere da quello delle 
particelle u. secondo la relazione ben nota: 
W\dW 

G,+(W)dW= @ | —. )—— , 

„=(W) os) 1,25 ? 


e se assumiamo per lo spettro delle particelle u quello di Sands, corretto alle 
alte energie in base ai dati della componente u. sottoterra e integrato su tutta 
Vatmosfera, otteniamo per i mesoni x carichi lo spettro (°) 


G,+(W)dW = 2,55:10-W-»s dW, (W >2 GeV), 


—— 
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-1 


e da questo, facendo il rapporto: 


2 = G(W)dW ; 
1 = - = 0,53. 
(16) nt @,2(W)aw 9° 


Il valore trovato per il rapporto mesoni neutri e mesoni carichi risulta in 
accordo soddisfacente con i pit recenti dati sperimentali in proposito; per 
citarne qualeuno: 

CARLSON e coll. (*) danno 0,45 + 0,10, BRADT e coll. (12) 0,39 + 0,13, 
SALVINI e coll. (#) 0,36, Lovartı e coll. (5) 0,42, GREGORY e coll. (1%) 0,3. 
Questo ultimo pero si riferisce al rapporto fra mesoni neutri e particelle pene- - 
tranti e quindi rappresenta un limite inferiore del rapporto in questione. 


6. — Confronto con i dati sperimentali del gruppo di Bristol. 


Puö essere interessante, per finire, tentare un confronto dei dati finora 
raccolti con quelli forniti dai lavori con lastre nucleari dal gruppo di Bristol. 

Consideriamo dapprima lo spettro integrale della sorgente di fotoni, defi- 
nito tramite la (15), 


Tr, te) at =¢ | sue to} aw. 


7» 


ove 7, & uguale a 0,25 GeV e f, a 50/43 u.r., il primo essendo il limite infe- 

riore dei dati sperimentali forniti dal gruppo di Bristol ed il secondo la quota 

a cui gli stessi sono stati misurati. I caleoli danno 

(17) WIS, te) — 210 a — 1,29-10-* mesoni neutri con energia 
>0,25 GeV a 50 g-cm- per 
cm*-s-sterad , 


con o, densita dell’aria standard = 1,29-10-* g- cm”. 
Dal lavero di CAMERINI e coll. (77) invece si ha uno spettro differenziale di 


(2) H. L. Brapt, H. F. Kıarpox e B. PETERS: Helv. Phys. Acta, 23, 24 (1950). 

(4) G. Saryısı e Y. Kim: Phys. Rer.. 86. 592 (1952). 

(5) A. Lovarı, A. Mura, G. TAGLIAFERRI e S. TERRANI: Nuovo Cimento, 9, 946 
(1952). 

(5) B. P. Grecory e J. H. Trxtot: Phys. Rev. 81. 675 (1951). 

‘(7) U. Camerısı, P. H. Fowrer, W. O. Lock e H. Murrueap: Phil. Mag., 41, 
413 (1950). 
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mesoni carichi emessi da stelle ottenuto da misure su lastre nucleari esposte 
in pallone a 50 g-em?. Da tale spettro differenziale, che & definito da 0,25 GeV 
a 1,4 GeV e che abbiamo approssimato con la espressione analitica 


W- per W<1 GeV, mesoni carichi per GeV 
a*(W, to) = 213-107) } d’intervallo per stella in 


| W-28 per W>1 GeV, emulsione a 50 gem”, 


abbiamo calcolato Vintensita integrale da 0,25 GeV in su: 


IT* (yo, to) = 4,13:10-1  mesoni carichi di energia > 0,25 GeV 
per stella in emulsione a 50 g:em ?. 


Da questo, usando del rapporto r°/7* — 0,45 dato da CARLSON e coll. (5), 
otteniamo 


(18) IT°(yo, to) = 1,86-10-1 mesoni neutri di energia > 0,25 GeV 
per stella in emulsione a 50 g-cm?. 


Per avere un dato confrontabile con la (17) dobbiamo moltiplicare la (18) 
per il numero », di stelle prodotte in un cm? di aria al secondo. I dati speri- 
mentali si riferiscono generalmente al numero di stelle prodotte in un em? 
di emulsione al giorno; d’altra parte il rapporto fra il numero di stelle prodotte 
in aria e quello di emulsione », si pud calcolare con la: 


Va ze > Nialia 
Feng ET, 
Ve >} Ne Oje 


ove n, e o, indicano rispettivamente il numero di nuclei per cm? e le sezioni 
d’urto di produzione di stelle dei vari elementi che compongono l’aria e l’emul- 


sione. Assumendo », = (2,1:103/8,64-10%) — 2,43-10-2 (18) e per le o, i valori 
delle sezioni d’urto geometriche (1°), si ottiene 


DE I 0 


Dividiamo ancora per 27 e facciamo Vipotesi che il numero di mesoni uscenti 


(5) M. CECCARELLI: comunicazione privata. 
(1%) J. J. LORD: Phys. Rev., 81, 901 (1951). 
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da stelle in aria sia lo stesso di quelli uscenti da stelle in emulsione; otteniamo 


1 
(19), IT (re to) — 186 108 1,29"107- ES 0,38-10-5 mesoni neutri di energia 


> 0,25 GeV a 50 g'em ? 
per em?-s-sterad . 


Il disaccordo non eccessivo tra il valore cosi ricavato e quello che abbiamo 
:alcolato con la (17), pud molto probabilmente ascriversi alle ipotesi grosso- 
lanamente approssimate, data la attuale insufficienza dei dati ad esse relativi, 
introdotte nell’ultimo ragionamento: in particolare quelle riguardanti la distri- 
buzione angolare dei mesoni r° prodotti, il rapporto delle sezioni d’urto di 
produzione di stelle e del numero medio di mesoni per stella in aria ed in 
emulsione. 

titeniamo percid che il limite di incertezza introdotto dalle precedenti 
ipotesi sia attualmente abbastanza grande da non consentire di attribuire un 
valore significativo alla differenza nei fattori numerici delle (17) e (19) ne di 
poter concludere, da tale disaccordo, sulla probabilita di esistenza di un’altra 
sorgente della componente elettrofotonica non contemplata nel precedente 
quadro. 


Ringrazio i proff. N. DALLAPORTA e G. PUPPI per avermi suggerito Videa 
del presente lavoro e per gli utili consigli e discussioni e per lPinteressamento 
con cui Vhanno seguito. 


SUMMARY 


The contribution from the 7° meson decay into two photons to the electronic 
component of the cosmic radiation, in passing through the atmosphere, is calculated; 
the diffusion equations are used with a generating function for the photons whose 
energy dependance is taken from the experimental data of the Bristol group, and whose 
variation with the altitude is derived from an absorption length of 130 g:em ?. 
The diffusion equations are integrated by a method of successive approximations, 
taking into account both the ionization energy loss and the variation of the cross 
sections for pairs and bremsstrahlung with energy. The results are normalized on 
the integral spectra of the total electronic component taken from the more recent 
experimental data and a ratio 7°/m* = 0.53 of neutral to charged mesons is obtained. 
The absolute values of the intensity obtained in this way are not inconsistent with 
those given by the Bristol group, within the wide limits of the incertainties actually 
existing in the knowledge of some processes used in the comparison. 
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Theoretical Resonance Curves in the Gyro-Interaction 
of Electromagnetic Waves in the Ionosphere. 


F. H. HIBBERD 


Department of Physics, University of Sydney, Australia 


(ricevuto il 16 Gennaio 1953) 


Summary. — In a recent communication Dr. Morzo has considered the 
gyro-interaction of radio waves in the ionosphere when the disturbing 
wave is incident vertically. She concludes that for this case one would 
expect to observe always a single maximum in the index of interaction 
as the frequency of the disturbing wave is varied through the gyro- 
frequency. It is possible to show that with larger approximations in the 
calculations either single or double maxima may occur and that 
there is little difference in the general shape of the resonance curves 
resulting from vertically incident and obliquely incident disturbing 
waves. A qualitative outline is given of the theory of gyro-interaction 
of V. A. BAILEY. 


1. — Introduction. 


In a recent communication M. MoTzo (1) has given theoretical curves for 
the index of interaction of radio waves in the ionosphere when the disturbing 
gyro-wave is incident vertically. These curves purport to show that for ver- 
tical incidence, using data appropriate to Naples, only a single maximum is 
to be expected in the magnitude of the transferred modulation as the fre- 
quency of the disturbing wave is varied about the gyro-frequency. It will 
be recalled that the original curves drawn by BAILEY (?), for a vertical magnetic 
field and a wave incident at 45°, show that if the wanted wave is reflected 
at à level below some 3 to 4 km above the bottom of the E-layer a single 
maximum is to be expected, whereas two maxima, separated by a minimum 

(7) M. Morzo: Nuovo Cimento, 9, 213 (1952). 

(2) V. A. Bartey: Phil. Mag., 26, 425 (1938). 
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which is very close to the gyro-frequency, result if the wanted wave is reflected 
at a higher level. Both single and double maxima have been observed expe- 
rimentally by CUTOLO (#) and by BAILEY (*) and their respective co-workers. 
The conclusion that for a certain direction of incidence of the disturbing wave 
the interaction curve has always a single maximum is unexpected and therefore 
the question has been examined in greater detail. 


2. — The resonance curves for the index of interaction. 


In the absence of specification of the wanted wave one cannot estimate 
the amount of modulation impressed on it. The interaction potentiality of a 
given disturbing wave can, however, still be calculated. It is expressed in 
terms of the index of interaction. This index is the quantity Nvw, where 
N is the electron concentration, » is the electron collision frequency and w 
is the mean work done by the disturbing wave on an electron in the interval 
between collisions. 

On general grounds, providing there is a resonance in the absorption of 
a gyro-wave (*), one would always expect both single and double maxima in 
the index of interaction. Further, the direction of propagation relative to 
the direction of the magnetic field is more important than the question of 
oblique or vertical incidence. Accordingly, with the same data as MoTzo’s, 
for vertical incidence at Naples, the curves representing the index of interaction 
have been recalculated. 

The expressions and data used in the computations are given here again 
for reference. The mean work # done by the disturbing wave on an electron 
per collision with the gas molecules is given by 


where t = e?/m = 2.53 - 108 (e.u.s.)?/gm; 
Z = R.M.S. value of the electric field of the wave at the point of inci- 
dence, estimated by Morzo as 2.2 mV/m for an 800 watt 


transmitter ; 
yp = angle between the electric vector of the wave and the direction 


of the earth’s magnetic field; taken here as 60°; 


(3) M. Curoro, M. CARLEVARO and M. GuerGui: Alta Freq., 15, 111 (1946); M. Cu- 
Toto and R. FERRERO: Nature, 163, 58 (1949); M. CuroLo: Nature, 166, 98 (1950). 

(4) V. A. Baıtey, R. A. Smith, K. LANDECKER, A. J. Hiaas, and F. H. HıBBERD: 
Nature, 169, 911 (1952). 

(*) For the values of N, v and Q corresponding to the lower E-layer, resonance 
in the absorption of the extraordinary component of a gyro-wave always occurs. 
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y — frequency of collisions of electrons with gas molecules; 
6 = w--.Q, where: 
m = angular frequency of the disturbing wave; 
Q — local gyromagnetic angular frequency of the electrons, taken here 
as 7.387 X106 (equivalent to 225 m); 
al mK, ee 
The quantity p’ in the expressions for w is given by 


te? 0 — cot? tN(1 + cos? 0) 
p' — 5 aux E; Fi (17 a ESS 
te? 0 + 2 DE 


where 6 = angle between the direction of propagation of the disturbing wave 
and the earth’s magnetic field; taken here as 30°; 
N = electron concentration. 
Both » and N vary with height and it is assumed that they vary in the 
following manner: 


y = % exp [— ay] 


N = N, exp [by] 


nN 
à 
7 > 
? erg cm” ] 


15 


10 


where y is measured in cm upwards from 
the lower boundary of the E-layer, taken 
as a height of 90 km, and » and N, are 
the values at this base level. The values 
adopted in these expressions are 1) = 1.1: 
10° Fa RO EN — 100 per cm?) db 1053 

The quantity a, measures the atte- 
nuation of the disturbing wave and is 


= | > 
2 BLOC eae given by the relation 


Ww 


Nr 


Fig. 1. — Graphs of Index of Inte- N 

raction Nvw against frequency (in J = { 

terms of e) of disturbing gyro wave. log. a, = — >| k ds = 

The dotted line indicates the gyro- d 

frequency. The curves are for va- . Nop 

rious heights above the lower bound- = — 0.0566 (tg «— cot x) | nee dy, 


ary of the E-layer. R 
where k = absorption coefficient of the wave, 
s = length of the ray-trajectory in the ionosphere, up to the point 
considered ; 
and tg « == cos 0. 
The results of these calculations are represented by means of the curves 
shown in fig. 1. 
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It is thus seen that when the curves given by Morzo are revised they be- 
come very similar to those originally given by BAILEY for oblique incidence. 
It then appears that the transition from a single to a double hump, as the 
height of reflection of the wanted wave increases, will always occur, whatever 
the angle of incidence of the disturbing wave. Further, the asymmetry about 
the gyro-frequency is not very great. 

Although any set of curves such as that shown in fig. 1 is applicable only 
to a given locality and direction of propagation, it is believed that these revised 
curves will be of value in connection with the researches in gyro-interaction 
which are proceeding in Italy. It should be noted, however, that the curves 
depend on the values adopted for b, the constant which specifies the distri- 
bution of the electron density, and for », the electron collision frequency, in 
the lower part of the E-layer. Also there are approximations in the theory 
of BAILEY (?), though current investigations indicate that the corresponding 
errors are not serious. 


3. — Theory of gyro-interaction. 


It is useful now to outline the basic physical ideas of BAILEY’S theory (2), (5) 
of gyro-interaction. From this it will be seen that in general one would always 
expect to find either single or double maxima, according to the height at which 
the wanted wave is reflected. 

When energy from a wave (the disturbing wave) is absorbed in a region 
of the ionosphere, the mean energy of agitation of the electrons is increased 
and this in turn increases the electron collision frequency. This will change 
the absorbing capacity of that region of the ionosphere. When the amplitude 
of the disturbing wave is modulated there will be a resulting modulation of 
the absorbing capacity of the ionosphere. Should a second wave (the wanted 
wave) suffer absorption in the disturbed region this wave will have impressed 
upon it a modulation whose fundamental frequency is the same as that of 
the modulation of the disturbing wave. When the wanted wave is not spe- 
cified one can determine only the effect which a given disturbing wave will 
produce in the ionosphere. This is best described in terms of the amount 
of energy absorbed from the disturbing wave. For this purpose BAILEY has 
chosen the quantity Nvw and has termed it the «index of interaction ». From 
the definitions it will be seen that this quantity represents the mean energy 
absorbed by the electrons per unit volume per second. BAILEY has shown 
that Nvw = 2kU, where k is the absorption coefficient of the disturbing wave 


(5) V. A. BaILEY: Phil. Mug., 28, 929 (1937). 
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and U is the mean value of its energy flux; this is in agreement with the inter- 
pretation of Nvw given here. 

A considerable increase in the amount of transferred modulation may 
result when the carrier frequency of the disturbing wave is varied through 
the gyro-frequency. This is because the extraordinary component of a wave 
suffers greatly increased absorption at the gyro-frequency. The absorption of 
the extraordinary component of a gyro-wave is so much greater than that 
of the ordinary component that only the 
effect of the extraordinary component need 
be considered. At any given height in the 
E-layer, w is determined by the strength 
of the disturbing wave at that height, and 
by the fraction of the wave energy that 
is absorbed. With increasing heights in 
the layer, measured from the lower boun- 
dary, the fraction of the wave energy that 
is absorbed increases rapidly but the wa- 
ve itself becomes weaker. For a given 
disturbing wave frequency, the index of 
interaction N vw attains a maximum value 
at a definite height, as shown by each 
curve in fig. 2. 

The absorption of the extraordinary 
component has a rather sharp maximum at 
the gyro-frequency and so the absorption 
of a wave on the gyro-frequency is much 
Fig. 2. — Graphs of the Index of In- greater than that of a wave whose fre- 
teraction against height inside the quency is a few percent less or greater than 
E-layer. The disturbing wave fre- the oyro-frequency. Thus the index of 
quency corresponding to each eurve . on Fa hee : 
is indicated by the adjacent value interaction attains its maximum value at 
of e. The curves show that the height, à Smaller depth of penetration for a wave 
at which a maximum in the index of whose frequency is equal (or, more 
interaction occurs, increases as the strictly, almost equal) to the giro-fre- 
wave frequency departs from the quency than it does for a wave whose 

sys onfoauency: frequency is slightly removed from the 
eyro-frequency. This is illustrated in 


fig. 2, in wich are drawn the curves corresponding to the gyro-frequency 
(e = 0) and to positive values of e. The curves for negative values of 
e are very similar to those for positive e and are omitted for clarity. At 
levels not far above the lower boundary of the E-layer, the index of 
interaction is a maximum for the gyro-frequeney, whereas at higher le- 
vels the index may be less for the gyro-frequency (because the wave has 
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become so weak) than for frequencies a little less or a little greater than 
the gyro-frequency. Thus, as shown in fig. 1, the curves representing the 
index of interaction have a single hump at levels just inside the lower E-layer 
and this single hump changes over to a double hump at higher levels. In 
fig. 2, the region below the dotted horizontal line corresponds to the single 
hump, the region above the line corresponds to the double hump. Whether 
one obtains experimentally the single or the double maxima in the depth of 
transferred modulation depends on the height at which the wanted wave suffers 
its greatest absorption. CuToLo (5) has actually observed a transition from 
the one type to the other as he varied the carrier frequency, and hence the 
height of reflection, of the wanted wave. 


(5) M. Curoro: Nuovo Cimento, 9, 391 (1952). 


RIASSUNTO 


In un recente articolo della Dott.ssa Morzo é stata considerata la girointerazione 
di onde nella ionosfera nel caso di un’onda disturbatrice con incidenza verticale, e si 
€ dimostrata l’esistenza di un solo massimo dell’indice di interazione allorché si varia 
la frequenza disturbatrice intorno alla frequenza giromagnetica. Si & pensato allora di 
fare ulteriori approssimazioni nei calcoli e pertanto si dimostra come & possibile cosi 
ottenere due curve; una con un solo massimo e una con due massimi. Inoltre si di- 
mostra che l’andamento generale delle curve di risonanza & lo stesso sia per incidenza 
verticale che per incidenza obliqua dell’onda disturbatrice. Si dä un breve riassunto 
delle linee fondamentali della teoria della girointerazione di V. A. BAILEY. 


IL NUOVO CIMENTO VoL: X, N. 4 1° Aprile 1953 


Misure di velocitä delle onde elastiche nei solidi 
a temperatura elevata. 


P. G. BORDONI e M. Nuovo 


Tstituto Nazionale di Ultracustica «O. M. Corbino» - Roma 


(ricevuto il 17 Gennaio 1953) 


Riassunto. — Si & misurata, con un metodo elettroacustico, la velocita 
delle onde elastiche estensionali in cinque metalli puri (Pb, Sn, Bi, Al, Cd). 
In base ai dati sperimentali si & calcolata la derivata logaritmica della 
velocita rispetto alla temperatura, confrontandone poi le variazioni con 
quelle del coefficiente di dilatazione termica relativo agli stessi metalli. 


1. — Seopo del lavoro. 


Le attuali conoscenze sulla struttura dei metalli permettono di attribuire 
le variazioni nella velocita di propagazione delle onde elastiche con la tempe- 
ratura a due gruppi di cause ben distinti tra di loro: 

1) agli effetti di rilassamento, dovuto, per esempio, al movimento delle 
dislocazioni nell’interno dei grani eristallini, agli scambi di calore tra i grani 
stessi, ai mutui scorrimenti del loro contorno, alla eventuale presenza di im- 
purita, e cosi via; 

2) alle variazioni nelle dimensioni medie del reticolo cristallino, prodotte 
dai cambiamenti di temperatura. 

I dati sperinientali di cui attualmente si dispone, riguardano prevalente- 
mente il primo gruppo; le misure relative alle variazioni di velocita in un ampio 
campo di temperature sono invece assai pit scarse, e riesce inoltre difficile 
collegarle tra di loro, essendo state in gran parte eseguite senza preoccuparsi 
di una possibile interpretazione fisica (1). 


(1) Cfr. ad esempio: W. KÖSTER: Zeits. f. Metallkunde, 39, 1 (1948). 
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Si è quindi ritenuto che un’indagine sistematica sulla dipendenza della 
velocita dalla temperatura, specialmente nei pressi del punto di fusione, ese- 
guita con la notevole precisione consentita dai metodi elettroacustici di eui 
attualmente si dispone, ed in condizioni sperimentali particolarmente semplici 
e chiaramente specificate, potesse fornire un utile contributo alla conoscenza 
della struttura fisica dei metalli. 

Nella presente nota sono raccolti i primi risultati ottenuti in tale indagine, 
i quali sembrano avere una portata più generale di quella che ci si potrebbe 
attendere, dato il ristretto gruppo di metalli al quale si riferiscono. Essi mo- 
strano infatti la possibilità di ricavare una legge generale di dipendenza della 
velocita dalla temperatura, analoga a quella relativa alla dilatazione termica 
e valida per temperature non troppo vicine a quella di fusione. Nei pressi di 
quest’ultima, i dati sperimentali suggeriscono invece la possibilità di una trat- 
tazione teorica della dipendenza delle proprieta elastiche dalla temperatura 
basata sulla teoria formale dei fenomeni di rilassamento. 


2. — Metodo di misura. 


I valori della velocita di propagazione delle onde estensionali (*) sono otte- 
nuti mettendo in vibrazione delle sbarrette cilindriche, prive di ogni sensi- 
bile vincolo esterno, e determinandone la frequenza di risonanza fondamen- 
tale »,. La velocita di propagazione c, per una sbarra infinitamente sottile, à 
data dalla formula (?) 


; 2 g2d2 
(1) es JE = Al | a = = 
dove: E = modulo di Young; o = densita del materiale; / = lunghezza della 
sbarretta; d = diametro della sbarretta; o — coefficiente di Poisson. 

Prima di ogni serie di misure le sbarrette vengono sottoposte ad un rin- 
venimento per quanto @ possibile completo, in modo da eliminare gran parte 
delle imperfezioni nel reticolo cristallino, dovute alle lavorazioni precedenti. 

Per ridurre al minimo i vincoli esterni le sbarrette, disposte verticalmente, 
sono sostenute da tre punte di acciaio del diametro di 1 mm, disposte a 120° 
tra di loro, e contenute nel piano nodale passante per il punto di mezzo della 
sbarretta. Tali punte penetrano in fori di diametro leggermente superiore, in 


(*) E sembrato opportuno usare tale espressione per mettere in evidenza il carat- 
tere di estensione semplice della sollecitazione, alla quale, fuori dell’asse della sbarretta 
non: corrispondono spostamenti puramente longitudinali. 

(2) Lorp RayLeiGn: Theory of sound (London, 1937) 2% ediz. vol. I, pag. 252. 
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modo che il contatto con la sbarretta avvenga soltanto su di una porzione della 
parete laterale dei fori stessi, e sia assicurato unicamente dal peso della sbar- 
retta. Si controlla che tale tipo di sostegno non ha alcuna influenza apprez- 
zabile sul valore di »,, in quanto quest’ultima frequenza rimane inalterata 
anche se le punte non appartengono esattamente al piano nodale. 

Durante le misure la sbarretta ed il suo sostegno sono collocati nella muf- 
fola cilindrica di un forno elettrico; la regolazione di temperatura, effettuata 
mediante termometri a mercurio provvisti di un contatto mobile, riduce le 
oscillazioni di meno di un grado; la misura della temperatura € fatta con termo- 
coppie ferro-costantana, tarate in corrispondenza dei punti di fusione di alcuni 
dei metalli su cui sono eseguite le misure (Sn, Bi, Cd, Pb, Zn); il raggiungi- 
mento dell’equilibrio termico tra la sbarretta e l’aria che la circonda & indi- 
cato dalla costanza (*) nel tempo della 1%. 

Le misure a bassa temperatura sono invece eseguite collocando la sbarretta 
ed il suo sostegno, racchiusi in un recipiente a tenuta, in una miscela di alcool 
etilico ed anidride carbonica solida, contenuta in un thermos e seguendo le 


TABELLA I. — Caratteristiche dei materiali. 
| | an te À may 
| poe Naeem. Sl Freq. | Densital Tempe- Tempe-| 
Prove. | prima del Lun- | Dia- | diriso- | 3 ravura | Eaves | 
Metallo vus | rinveni. | SheZza | metro | nanza | 20 °C di di 
| nie mm | mm |a200C \g-em— fusione | Debye | 
| eps OK NOK | 


| | | | 
Piombo /1.G.F. Merck Fusione e | 49,0 | 10,0 | 13813 | 11,35 | 600,4 | 88 


HARZ | tornitura | 
Stagno | C. Erba | Fusione e | 55,0 | 10,0 | 24545 | 7,28 | 504,9 | 160 
| (T.P.) | tornitura 


| Bismuto |I.G.F. Merck Fusione e | 72.0 10,0 |13872 | 9,80 | 544,3 | 100 
| | | | 


| (purissimo) | tornitura 


| Alluminio | Ist. Metalli, Fusione e 
Leggeri | tornitura | 
(99,7%) | 


Cadmio 1.G.F.Merek| Fusione e | 140,0 | 10,0 | 9000 | 8,642) 593,9 | 172 | 
| (in sbarre) | tornitura | | 
| HP) | | 


(?) Tale costanza indica inoltre l’assenza di fenomeni di rinvenimento, durante i 
quali, come è stato precedentemente segnalato, la », aumenta anche se la temperatura 
si mantiene costante. 
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lentissime variazioni di velocita e di temperatura che si hanno durante la fase 
di riscaldamento. 

L’apparecchiatura elettroacustica per Veccitazione e la misura delle vibra- 
zioni, & sostanzialmente quella stessa dettagliatamente descritta in una nota 
precedente (4) e gia adoperata in diverse altre ricerche. Sara quindi sufficiente 
ricordare che essa consente la determinazione di 7) con un errore relativo non 
superiore a 1 ~ 2-10-%, e con valori della deformazione dell’ordine di 10. 


3. — Misure effettuate ed interpretazione dei risultati. 


La Tab. I indica i materiali su cui sono state effettuate le misure, la loro 
provenienza, il processo di lavorazione subito prima del rinvenimento, le di- 
mensioni delle sbarrette e le rispettive frequenze. Per comodita si sono anche 


TABELLA II. — Valori sperimentali della velocità di propagazione di onde estensionali. 


| | 
| Alluminio Ist. Me- 


; Cadmio Merck 
| talli Legg. 99,7% 
| 


| 
| | 

Piombo Merck | Stagno C. Erba | Bismuto Merck 
| Pale 


| RER eng os ge PN 066 
DRS Sigs TO ON PE Sie oK misst | 8K" |} mis 


- I I—-| | je 


238 | 1398 | 289 | 2700 | 201 | 2044 | 245 | 5121 | 298 | 2520 
246 | 1393 | 323 | 2650 | 218 | 2036 | 259 | 5092 | 339 | 2475 
250 | 1389 | 348 | 2603 | 243 | 2022 | 272 | 6071 | 368 | 2441 
285 | 1385 | 372 | 2663 | 255 | 2012 | 289 | 5050 | 409 | 2396:| 
268 (013760 303) 2619") 275 | 2003 | 303 | 5035 | 432 | 2834 


2 
270, 11137217 4264 | 24209) 287. 1.2000..| 318 5011 464 | 2247 
287 || 1361. |} 448) 1,2325 | 317 101982 335 4996 | =) 
304 | 1345 | 474 | 2164 | 318 1979 361 4957 | — |r — 
320 | 1331 | 492 | 2018 | 344 | 1963 | 368 4947 | — | — 
Sy ESTO Poke Be Lassie 530. [4999 ) = | ae 
374 | 1287 | — = 369 1946 391 4917 | — | — 
3992 1.1265.) =4 u OST 1,934,,,,.5418 4888 | — | — 
425 | 1243|. — = 57 390. .1;,1924 | 404 4836 | — | — 
451 | 122 | — AL BORN AS 8265 u RE 
AGE 208.) = = 416 1900 | 481 | 4797 | — | — 
Als 1198| 445 SA) 401, | AT. ee 
487 | 11887] — SM] 4704 1849 | 620 4740 | — | — 
519 | 1153 | — = 497 | 1821 | 552 4685 | — | — 
er  vosGn Ie He = | 620 -|.1788 | :578 | 46295) — hi 
= | | 531 1773 602 4586 | — | — | 
| | | | | 628 LAON EE DES 
| ; PP ON ES a le ace 


(4) Cfr. P. G. BORDONI: Nuovo Cimento, 4, 177 (1947); Rie. Seient., 8, 103 (1948). 


26 - Il Nuovo Cimento. 
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riportati i valori della temperatura di fusione, di quella di Debye e della den- 
sita a temperatura ambiente. 

1 risultati delle misure effettuate sono raccolti nella tab. II. Le curve c(t) 
che da essi possono essere ricavate, non si prestano perd ad istituire un con- 
fronto fisicamente significativo tra il comportamento dei diversi metalli; sembra 
invece preferibile considerare anziché i valori assoluti di Te dic o, come altri 
ha fatto, di de/d7’, i rapporti di tali grandezze rispettivamente alla tempera- 
tura di fusione 7, ed alla velocita limite ec, cui tende ¢ quando la 7 si avvi- 
cina allo zero assoluto. 

Per evitare le notevoli difficolta sperimentali inerenti alla determinazione 
di &, basta osservare che sussiste Videntita 


e 
ir, er 


la quale mostra che la derivata logaritmica della curva sperimentale c(T) ha 
un significato più generale della curva stessa. Infatti il coefficiente A definito 
dalla (2) ha per la velocita un 
ufficio analogo a quello che il 
coefficiente di dilatazione «,= 
= (0 log V/öT) ha per il vo- 
lume: come «, consente di pa- 
ragonare tra loro le dilatazicni 
‘di corpi aventi volumi diffe- 
renti, cosi ß permette il con- 
fronto tra le variazioni die an- 
che per metalli che abbiano va- 
lori assoluti della velocita molto 
diversi tra di loro. Natural- 
0,0 0,2 0,4 0,6 0,8 10 mente, trattandosi di una deri- 
ig. 1. — Variazioni con la ten peratura di ße Vata grafica o numerica, la 
di «. per il Piombo, in funzione del rapporto P(T/T,) è affetta da errori al- 


della temperatura assoluta alla tenperatura di quanto maggiori di quelli re- 
fusione. 


1,507 (°K )"' 


1,25 


lativi alla curva sperimentale 
c(T), ma l’elevata precisione 
delle misure contiene tali errori nei limiti perfettamente accettabili di qual- 
che unita per cento. 

Nella fig. 1 sono riportati, in funzione di 7/7,, i valori di B e, a scopo di 
eonfronto, anche quelli del coefficiente di dilatazione cubica (5). 


(®) Cir. For. L. UFFELMANN: Phil. Mag., 10, 633 (1930); oppure P. HıDNErT, 
W. T. Sweeney: U.S. National Bureau of Standards Journ. of Rs., 9, 703 (1932). 
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Alle basse e medie (%) temperature le due curve hanno un andamento del 
tutto simile; in particolare si nota in entrambe un tratto abbastanza esteso, 
corrispondente a 0,2< T/T,< 0,75 nel quale le variazioni sono sensibilmente 


lineari, e le curve si avvicinano molto a 
delle rette orizzontali. Invece per T/T,> 0,75 
la pendenza della curva B(T/T,) aumenta for- 
temente, e tale comportamento non trova 
alcun. riscontro in quello della «,(7T/T,), la 
cui pendenza si mantiene quasi costante. 

Per lo stagno (fig. 2) i valori di B e di «, 
che si posseggono (7) sono relativi ad un 
campo di temperature assai più ristretto, e 
comunque tutto al disopra del punto di tra- 
sformazione dalla forma tetragonale a quella 
cubica (291 °K). Tali valori sono tuttavia 
sufficienti ad indicare anche per lo stagno 
l’esistenza di un intervallo di temperatura, 
limitato superiormente da 7/7, = 0,75, nel 


Fig. 2. — Variazioni con la tem- 
peratura di Be di «, per lo Sta- 
gno, in funzione del rapporto 
della temperatura assoluta alla 
temperatura di fusione. 


quale le curve «.(T/T,) e B(T/T,) hanno un andamento molto simile e sensi- 


bilmente rettilineo, con 
piccola inclinazione ri- 
spetto all’asse delle 


ascisse. Invece per 


T/T,>0,75 la pen- 
denza della curva 6 


aumenta, anche in 
questo caso, molto pit 
fortemente di «quella 


della curva relativa 


alla dilatazione termi- 
ca. Lo stesso anda- 


mento si riscontra pure 


0,4 0,5 0,6 0,7 0,8 0,9 


Fig. 3. — Variazioni con la temperatura di B per il Bi- 


nel Bismuto (fig. 3) 
per il qualei dati sulla 
dilatazione termica so- 
no tuttavia troppo in- 


10 


smuto, in funzione del rapporto della temperatura asso- certi, per consentire 


luta alla temperatura di fusione. 


un efficace confronto. 


(*) I valori relativi alle basse temperature sono desunti da una ricerca precedente: 
cfr. P. G. BORDONI: Hlastic behaviour of some face-centered metals at very low temperatures 


(in corso di pubblicazione). 


(7) Per i valori di «, cfr. il primo lavoro citato nella nota (5). 
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Pitt complicato appare invece l’andamento delle curve relative all’Alluminio 
ed al Cadmio (figg. 4, 5). In entrambe si nota innanzi tutto la presenza di un 
massimo abbastanza acuto, dovuto probabilmente a qualche fenomeno di ri- 

lassamento. Infatti considerazioni teori- 


09-(K" == > che abbastanza sempliei (5), mostrano 
08 —— ir ame ic che le curve c(») e c(T) possiedono un 
07|— ES ee on mes 


punto di flesso corrispondente rispetti- 
vamente à quei valori della frequenza, 
o della temperatura, per i quali il pe- 
riodo di oscillazione della sbarretta, as- 
sume un valore prossimo al tempo carat- 
teristico di rilassamento. Le derivate, 
sempliei o logaritmiche, di tali curve 
Fig. 4. — Variazioni con la tempera- rispetto a » oppure a 7, debbono quindi 
tura di Be di «, per l’Alluminio, in  Presentare, in corrispondenza al flesso, 


funzione del rapporto della temperatura dei massimi, come quelli che si osser- 
assoluta alla temperatura di fusione. vano nelle figg. 4 e 5. 


0.2 0,3 0,4 0,5 0,6 0,7 


Alle temperature in cui l’effetto di 
rilassamento non & sensibile (7 < 0,5 7,, per l’Alluminio; 7 < 0,65 T,, per il 
Cadmio) le curve «, e 6 hanno sensibilmente 
lo stesso andamento (?) in entrambi i me- 
talli. 

In definitiva il confronto tra le diverse 
curve che caratterizzano le variazioni rela- 
tive di velocita e le corrispondenti curve 
che caratterizzano invece le dilatazioni ter- 
miche, rivela, nonostante la presenza di ac- 
cidentali fenomeni di rilassamento, l’esistenza 
di una sostanziale similitudine nel compor- 
tamento dei diversi metalli riguardo alle va- 


Fig. 5. — Variazioni con la tem- 


riazioni di velocita ed alle dilatazioni ter- peratura di ß e di «, per il Ca- 
miche, almeno per temperature minori di dmio, in funzione del rapporto 
0,7 T,. Dal punto di vista teorico tale simi- della temperatura assoluta alla 


litudine appare perfettamente giustificata dal temperatura di fusione. 


fatto che le variazioni di volume e di velo- 
cita sono entrambe dovute ad una stessa causa, e cioè alla non-linearitä delle 


(8) Cfr. ad esempio: C. 7ENER: Elasticity and anelasticity of solids (Chicago, 1948); 
oppure: P. G. BORDONI: Nuovo Cimento, 7, 144 (1950). 

(°) Per i valori di « relativi all’Alluminio cfr.: Tables annuelles internationales de 
constantes et données numériques (Paris (1931-1934), 11, 1, n. 4, 2. Per il Cadmio cfr.: 
International Critical Tables (New York, 1933), 1, 459. 
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relazioni che legano — in media — le azioni interatomiche alle dimensioni 
del reticolo. 

I risultati sperimentali ottenuti suggeriscono quindi la possibilita di rica- 
vare dalla equazione di stato dei solidi, analogamente a quanto gia si & fatto 
per «, (1°) una relazione teorica che dia f in funzione della temperatura. Non 
sembra perd che una tale relazione possa giustificare anche lo scostamento 
della curva B(T/T,) dalla «.(T/T,) che si verifica nei pressi del punto di fusione. 
Infatti 7; non compare direttamente nell’equazione di stato, e d’altra parte 
le ipotesi semplificative poste a base di quest’ultima, divengono sempre meno 
valide man mano che cresce la temperatura. 

Sembra piuttosto che il forte aumento della pendenza che si riscontra nella 
curva B(T/T,) quando ci si avvicina al punto di fusione indichi l’esistenza di 
qualche fenomeno di rilassamento, il cui massimo, alle frequenze considerate 
pud anche corrispondere a temperature maggiori di T,. 


4. — Conelusioni. 


Le misure di velocita delle onde elastiche eseguite in funzione della tem- 
peratura ed a frequenze comprese tra 5 kHz e 25 kHz su Pb, Sn, Bi, Al, 
Cd, mostrano che: 


1) Le curve che danno — dlog e/öT in funzione di T/T, hanno sostan - 
zialmente lo stesso andamento per tutti i metalli considerati, nonostante la 
presenza, nei provini di Al e di Cd, di fenomeni di rilassamento il cui effetto 
& d’altronde sensibile soltanto in un ristretto intervallo di temperatura. 


2) Tali curve possiedono un tratto intermedio, abbastanza esteso, che 
appare sensibilmente rettilineo e poco inclinato sull’asse delle ascisse. 


3) Per ognuno dei metalli considerati e per T<0,1T,, la curva che 
da — Ologe/öT ha inoltre lo stesso andamento di quella che caratte- 
rizza 0 log V/oT. 


4) Quando ei si avvieina alla temperatura di fusione la pendenza della 
curva — log JT aumenta fortemente, mentre nulla di simile si verifica 
nella clog V/eT. 


(22) Cfr. ad esempio: J. €. Starter: Introduction to chemical physies (New York, 
1939), p. 219. 
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SUMMARY 


The velocity ¢ of extensional waves has been determined as a function of temper- 
ature 7 for five pure metals (Pb, Sn, Bi, Al and Cd). The vibration frequency had 
its lowest value for Cadium (9000 cps at 20°C) and the highest one for Aluminium 
(25 483 eps at 20 °C); the strain was of the order of 10-8. From the experimental data 
the coefficient 6 = — 0 log c/2T has been computed, and the curves B(T/T,) traced T, 
being the melting point temperature. From the inspection of the above curves, and their 
comparison with the corresponding ones for the thermal expansion coefficient 
a, = Olog V/0T suggest the following conclusions: 1) For the five metals investi- 
gated the 6(7/T.) curves look substantially alike, notwithstanding some relaxation phe- 
nomena, which occur in the Aluminium and Cadmium samples. 2) All these curves in 
an intermediate temperature range are straight lines and in the same range ß 
keep almost constant. 3) For every metal among the five considered, and for 7 < 0,77, 
the 8 and «,-curves are rather similar to one another. 4) When 7 gets near the melting 
point value the $ curve becomes more and more steep, while nothing of the kind happens 
in the «,-curve. 
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Diffrazione dei raggi X nei liquidi. 


. CENNAMO 


Istituto di Fisica Sperimentale dell Universita - Napoli 


(ricevuto il 19 Gennaio 1953) 


Riassunto. —- 


Si descrive il metodo sperimentale seguito per la misura 
dell’intensitä dei raggi X diffratti dai liquidi e si riportano i risultati rela- 
tivi al benzolo ed al cicloesanolo. 


Allorquando un fascetto monocromatico di raggi X attraversa un liquido, 
l’intensità dei raggi da questo diffratti presenta uno o più massimi in corri- 
spondenza di particolari valori degli angoli di diffusione che dipendono sia 
dal liquido che dalla lunghezza d’onda della radiazione impiegata: interponendo 
una lastra fotografica normalmente al percorso del fascetto che attraversa il 
liquido si ottengono anelli di diffrazione più o meno sfumati. 

Questo comportamento mostra che nella sostanza allo stato liquido resta 
indubbiamente una regolarità ed infatti la distanza tra le molecole di una 
sostanza allo stato liquido calcolata per diffrazione dei raggi X, e la distanza 
tra le molecole della stessa sostanza allo stato solido determinata col metodo 
delle polveri sono tanto vicine che in prima approssimazione si possono 
considerare le stesse. 

Diverse teorie sono state proposte per spiegare la diffrazione dei raggi X 
prodotta dai liquidi (KEESOM e DE SMEDT (1), C. V. RAMAN e K. R. RAMA- 
NATHAN (2), DEBYE (3), P. ZERNIKE e J. A. PRINS (*), H. MENKE (5), G. W. 
STEWART (5), ece.), ma qualunque possa essere l’interpretazione € evidente che 


(4) W. H. Kersom e J. DE SMEDT: Proc. Acad. Amsterdam, 27, 112 (1923). 
(A) Ge en “Raman e K. R. RAMANATHAN: Nature, 3, 185 (1923). 

(3) P. DeBye: Phys. Zeits., 28, 135 (1927). 

(4) P. ZERNIKE e J. A. Prins: Zeits. f. Phys., 41, 184 (1927). 

6) H. Menke: Phys. Rev., 33, 593 (1932). 

(6) G. W. Stewart: Phys. Rev., 32, 153 (1928). 
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la presenza di anelli di diffrazione nei liquidi sta ad indicare la presenza di 
piani reticolari e quindi di raggruppamenti regolari dovuti alla permanenza 
di forze intermolecolari più deboli, ma della stessa natura di quelle esistenti 
allo stato solido. In aleune note precedentemente pubblicate (7) noi abbiamo 
riportato l’andamento delle intensita dei raggi X diffratti da alcune sostanze 
allo stato liquido sia al variare della temperatura sia al variare della concen- 
trazione in aleuni miscugli: ci proponiamo 
ora di descrivere il metodo seguito in queste 
misure. 

Tl dispositivo sperimentale utilizzato à 
illustrato in fig. 1: la scatola A, a pareti di 
piombo dello spessore di 7 mm, rigidamente 
fissata al tubo a raggi X a mezzo del so- 
stegno S, presenta nella sua parete anteriore 
un collimatore costituito da un tubo cilin- 
drico cavo di piombo sulle cui basi € prati- 
cato, in corrispondenza dell’asse, un foro del 
diametro di 1 mm. Un terzo foro, in asse coi 
precedenti, € praticato sulla parete della sca- 
tola ed è, dalla parte rivolta verso la va- 
schetta contenente il liquido da analizzare, 
opportunamente svasato. Il fascetto di 
ragei X, cosi delimitato, dopo aver attra- 
versato il liquido posto nella vaschetta V, 
giunge sulla lastra fotografica ZL; lungo il percorso perd incontra un di- 
schetto di piombo D del diametro di 4 mm e dello spessore di 3 mm con 
al centro un foro di 2 mm, chiuso da una laminetta di piombo di spessore 
opportuno: la intensita del fascetto incidente & cosi ridotta in modo che l’an- 
nerimento da essa prodotto risulti uguale all’ineirea a 0,6-0,7 e cioè perfetta- 
mente misurabile. 

Al disotto del cilindro © sulla parete anteriore della scatola una fenditura 
orizzontale F, lunga 15 mm ed alta 5 mm si prolunga nell’interno in un tubo E 
a pareti molto spesse della stessa sezione. 

Lungo E, il cui asse é opportunamente inclinato rispetto al tubo a raggi X, 
passa quindi un fascetto di raggi X, ridotto in intensità da una laminetta di 
piombo di spessore opportuno e da un filtro di nichel dello stesso spessore di 
quello utilizzato per rendere monocromatico il fascetto che attraversa il colli- 
matore. 


(7) F. Cennamo: Rend. Acc. Lincei, 9, 339 (1950); 10. 310 (1951), 10, 475 (1951); 
12, 294 (1952). 
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Una spessa lamina di piombo G, nella quale sono incise sette fenditure, 
posta all’estremo del tubo E, permette l’uscita di sette fascetti di uguale inten- 
sita: l’uguaglianza delle intensita dei fascetti emergenti dalle fenditure è stata 
sempre assicurata dal fatto che i rispettivi annerimenti prodotti sulla lastra 
fotografica, misurati al microfotometro, risul- 
tano uguali. 


Ss 

. A . I 

Disponendo invece lungo queste fenditure 100 
delle laminette di piombo di spessore diverso, 

i fascetti emergenti dalla lamina @ risulte- so 


ranno di diversa intensità e quindi in corri- 
spondenza sulla lastra fotografica produrranno 0 2 4 B DR 
delle tracce con annerimenti diversi. Noi ci 


siamo avvalsi di queste marche per costruire 

una curva di taratura la quale consentisse di 

trasformare la distribuzione degli annerimenti prodotti dai raggi diffratti dal 
liquido sulla lastra fotografica in distribuzione d’intensitä, ed allo scopo 
abbiamo proceduto nel seguente modo. 

Com’é noto, i raggi X, per quanto riguarda la loro azione fotografica, solo 
apparentemente sono simili alle radiazioni visibili ed a quelle ultraviolette: 
l’esame microscopico di uno strato sensibile esposto all’azione dei raggi X 
mostra. a differenza della luce visibile, che i raggi X producono attraverso 
l'intero spessore dello strato sensibile un’uguale distribuzione di granuli di 
argento depositato: variando l’intensità del fascetto incidente varia il numero 
di granuli di argento depositati per unita di volume, ferma restando la uni- 
formità di distribuzione per l’intero spessore dell’emulsione. Pertanto una 
lastra fotografica, esposta all’azione dei raggi X, esaminata per trasparenza 
si comporta come una lamina a spessore costante ma con coefficiente di assor- 
bimento variabile da punto a punto a seconda della intensita del fascetto di 
raggi X che ha colpito ciascun punto della lastra. 

Indicando con D, la deviazione segnata dal galvanometro del microfoto- 
metro per un punto non colpito dall’azione dei raggi X, e con D la deviazione 
misurata invece in corrispondenza di un punto colpito dall’azione delle radia- 


x 


zioni, l’annerimento S & definito dall’espressione seguente: 
B D, 
(1) S = log io Di: 


_ L’annerimento S cosi definito varia in funzione del tempo di esposizione t, 
delVintensita I del fascetto incidente e della lunghezza d’onda 2 della radia- 
zione, per cui si puo scrivere: 


(2) S = fiL,tA). 
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A parita di J e t, Vannerimento & una funzione piuttosto complessa di 2 
come si pud osservare nella fig. 2 riportata dal lavoro di R. BERTHOLD ed R. 
GLOCKER (8): S decresce con legge lineare molto lentamente per frequenze 
crescenti sino a v = 3,23-1018, aumenta bruscamente in corrispondenza di tale 
frequenza, decresce con la stessa legge nell’intervallo di frequenze 3,23-6,00-1018, 
subisce in corrispondenza di tale frequenza un ulteriore aumento brusco ed 


I- 560 
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Pig. 3: 


infine decresce con l’aumentare della frequenza sempre con legge lineare molto 
prossima a quella relativa ai due intervalli precedenti. 

Le brusche variazioni sono dovute all’assorbimento selettivo del bromo e 
dell’argento dell’emulsione e si manifestano in corrispondenza del passaggio 
da lunghezze d’onde superiori a lunghezze d’onda inferiori a quelle corrispon- 
denti ai rispettivi massimi di assorbimento di queste due sostanze. 

Poiche, come si & fatto osservare, la diminuzione di S colle lunghezze d’onda 
€, per ciascuno dei tre intervalli considerati, molto piccola, l’annerimento si 
puö, per lo meno in prima approssimazione, in ciascuno di questi tre inter- 
valli considerare indipendente da A, ed inoltre, determinato per uno dei tre 
intervalli spettrali l’andamento della funzione S = f(/,t) si potrà ricavare 
quello relativo agli altri due intervalli moltiplicando le intensita corrispon- 
denti ai diversi annerimenti per un coefficiente numerico facilmente ricavabile 
dalla fig. 2. 

A tali criteri ci siamo ispirati nella determinazione delle curve di taratura 
delle lastre a mezzo delle marche precedentemente descritte. 

Lo spettro di un ordinario tubo a raggi X con anticatodo di rame contiene, 
oltre alle radiazioni CuK, e CuK, caratteristiche anche delle radiazioni con A 
più piccole, ed il filtro di nichel, che con lo spessore da noi adoperato rende 


(8) R. BERTHOLD e R. GLOCKER: Zeits. f. Phys., 31. 261 (1925). 
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sufficientemente monocromatico il fascio agli effetti della diffrazione dei liquidi, 
eliminando la CuK, lascia passare, ed in buona parte, anche queste radiazioni. 

In fig. 3 si & riportata la distribuzione spettrale dei raggi X emessi da un 
anticatodo di Cu per una tensione applicata agli elettrodi di 44 kV (°); nella 
fig. 4 invece é riportata la stessa distribuzione spettrale filtrata perd attra- 
verso una lamina di nichel dello spessore di 10-2 em da noi utilizzata. Il sue- 
cessivo filtro di piombo, usato per la costruzione delle marche, riduce del tutto 
la intensita della radiazione CuK, e molto poco quella delle radiazioni dure e 
cid in dipendenza del fatto che la curva di assorbimento del piombo presenta 
valori molto elevati per À — 1,4 A ed un massimo molto meno accentuato 
per À — 0,7 À. 

L’annerimento delle marche di controllo si pud quindi considerare prodotto 
da radiazioni le cui lunghezze d’onda sono comprese nell’intervallo spettrale : 
À = 0,5-0,2 À. 

Com’e noto, l’intensità del fondo continuo emesso da un anticatodo varia 
proporzionalmente al quadrato della tensione applicata al tubo, mentre Vin- 
tensita della radiazione monocromatica caratteristica varia all’incirca propor- 
zionalmente a (V — V,)#2, in cui V, è la tensione di eceitazione della radia- 
zione monocromatica: mantenendo quindi costante sia la tensione applicata 
che la corrente che passa nell’ampolla si puö ritenere che resti costante la 
distribuzione spettrale emessa dall’anticatodo e quindi quella relativa alle pic- 
cole À trasmessa attraverso i due filtri. 

Poiche nell’intervallo spettrale innanzi citato S si pud considerare indi- 
pendente da A, potremo, in queste condizioni, ritenere l’annerimento solo fun- 
zione di J e di t e precisamente si pud porre: 


(2) Selle tn, 


essendo p un coefficiente numerico che tende ad 1 per piccole intensita. 
Perö per determinare il valore dell’intensitä 7, del fascio emergente dalla 
i-esima fenditura non & evidentemente possibile avvalersi deila relazione: 


tf; = I, exp [|= un], 


nella quale 7, rappresenta l’intensitä del fascio incidente su ciascuna delle sette 
fenditure, u il coefficiente di assorbimento del piombo, x, lo spessore di piombo 
attraversato, sia perché, come si è detto, il fascio non è monocromatico e 
quindi essendo u = f(a, A) i calcoli sarebbero piuttosto laboriosi e sia perchè 
la parte assorbita dallo schermo: 


Lai = I,(1 — exp [— url), 


(9) J. Bouman: Selected Topies in X-Ray Crystallography (Amsterdam, 1951), p. 116. 
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è, com’è noto, in parte riemessa sotto forma di radiazioni secondarie (radia~ 
zioni di fluorescenza, radiazione diffusa senza cambiamento di lunghezza d’onda, 
effetto Compton) e queste radiazioni secondarie apportano anche un loro con- 
tributo, sia pure modesto, all’annerimento della lastra fotografica che risulta 
pertanto maggiore di quello che si avrebbe, a parita 
dell’intensita del raggio trasmesso, in loro assenza. 
In altre parole, su di un punto A della lastra foto- 
erafica L (fig. 5), oltre al raggio con intensita J, tra- 
smesso attraverso lo spessore x, della laminetta di 
piombo, giungeranno anche le radiazioni parassite 
secondarie riemesse a causa dell’assorbimento attra- 
verso lo schermo e provenienti da tuttii punti di questo. 
Abbiamo pertanto preferito procedere alla taratura 
delle marche per via sperimentale ed allo scopo sono 
state eseguite, sulla stessa lastra, fotografie dei sette 
fascetti emergenti dallo schermo @ con tempi di posa 
diversi avendo, cura che durante l’esposizione le condizioni di funzionamento 
del tubo restassero rigorosamente costanti in modo da poter considerare tali sia 
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la composizione che l’intensità di ciascun fascetto. 
Riportando (fig. 6) sull’asse delle ordinate i logaritmi degli annerimenti 
prodotti da ciascuna fenditura e su quelle delle ascisse i logaritmi dei corri- 
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spondenti tempi di posa si sono ottenute, in accordo con la relazione: 


log S = plogt + costante 


delle rette. 


I coefficienti angolari di queste rette danno i valori di p corrispondenti a 


ciascuno degli spessori di piombo 
attraversati. 

In fig. 7 & riportato l’anda- 
mento di p in funzione di Aw: 
con l’aumentare di Ave quindi 
col diminuire dell’intensita del 
fascetto incidente sulla lastra, p 
tende al valore unitario. 

Determinati i valori di p e 
quindi delle espressioni 1?, Si ri- 
portano (fig. 8) sull’asse delle or- 
dinate gli annerimenti S e su 


quelle delle ascisse i corrispondenti valori t? e per ciascuna mar 
una retta il cui coefficiente angolare esprime, in base alla 
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500 


intensità del fascetto trasmesso attraverso la marca in esame. 
Nella fig. 9 sono riportati i risultati relativi a diverse marche d’intensitä 
ottenuti con lastre diverse, sempre nelle stesse condizioni di funzionamento 


del tubo, ma con tempi di posa 


ne 


20 x 10cm 


ca si ottiene 


valore della 


diversi: si ottengono, riportando sull’asse delle 
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ascisse gli spessori di piombo attraversati e sulle ordinate i logaritmi naturali 
dei rapporti S/t? = J, delle rette tra loro parallele, il cui coefficiente angolare, 
col segno cambiato, rappresenta il coefficiente angolare di cui si à innanzi 
parlato. Tale coefficiente & risultato uguale a K= 252 em. 

Se si tien conto del fatto che le sostanze attraversate sono il piombo ed il 
nichel, che le radiazioni trasmesse sono comprese tra A = 0,28 À e À = 0,50 A 


Ax 


Fig. 9. Fig. 10. 


e che in quest’intervallo per il piombo é, in media w,— 340 em”! e per il nichel 
si ha w, — 90 em, il valore di K ottenuto & abbastanza vicino a quello 
K'= 280 em che si ottiene dalla relazione: 


Sea 


Die Ta 


AC 


? 


nella quale x, ed x, rappresentano rispettivamente gli effettivi spessori di 
piombo e di nichel attraversati. Il valore sperimentale di A, inferiore a quello 
cosi ottenuto, pud giustificarsi pensando all’influenza delle radiazioni secon- 
darie di cui si € innanzi parlato. 

Poiché, anche nelle stesse condizioni di funzionamento del tubo, l’anneri- 
mento delle marche pud variare da lastra a lastra e per variazioni nel tempo 
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di posa e per le variate condizioni di sviluppo, si deve, al fine di rendere con- 
frontabili i risultati ottenuti con lastre diverse, nell’assegnare il valore del- 
Vintensita relativa alla marca più intensa, determinarne il valore J, per inter- 
polazione della curva annerimenti-intensitä relativa ad una lastra della stessa 


fig.11b 


serie e poi stabilire quello da assegnare alle altre marche di controllo dalla 
relazione: 


I, = I, exp [— ke. 


Nella fig. 10 @ riportata la curva di taratura delle lastre Ferrania tipo U, 
relative ad una serie di fotografie. 

Nelle fig. 11 a) e b) sono riportate le microfotometrie degli anelli di diffra- 
zione ottenuti per il benzene e per il cicloesanolo: la vaschetta nella quale era 
posto il liquido in esame era a pareti di mica, per cui sulla lastra fotografica 
oltre all’anello di diffrazione prodotto dal liquido compaiono in parte a questo 
sovrapposti i Laue prodotti dalla mica che perd, come del resto appare dalle 
fig. 11 possono, nell’eseguire la microfotometria dell’alone, facilmente evitarsi, 
scegliendo una direzione opportuna. 

Nella fig. 12 sono riportate le distribuzioni dintensita delle stesse sostanze 
in funzione degli angoli di diffrazione espressi in radianti; i valori ottenuti per 
le intensità nelle varie direzioni sono stati moltiplicati per l’espressione J/cos? 0. 
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Questo fattore di correzione € stato introdotto allo scopo di poter confrontare 
i valori ottenuti per le intensitä relativamente a punti diversi della lastra foto- 
grafica i quali, essendo quest’ultima piana, capitavano a distanze dal punto 
nel quale era posto il liquido colpito dal fascetto di raggi X diverse fra loro. 

Infine nella fig. 13 sono riportate le bande di diffrazione caratteristiche delle 
due sostanze, le quali sono state ottenute sottraendo il continuo che si estende 
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Fig. 12. Fig. 13. 


da 0,2 a 0,6 radianti. I massimi di diffrazione sono di poco spostati tra loro, 
pero la larghezza delle due bande di diffrazione, misurata in corrispondenza 
delV’intensita I’= %I,„ è notevolmente diversa: per il cicloesanolo si ha, 
come del resto appare dal confronto delle microfotometrie, una banda di dif- 
frazione molto piu stretta. 

Le due distribuzioni d’intensita relative a queste due bande di diffrazione 
verranno utilizzate in un prossimo lavoro in corso di pubblicazione. 


Lavoro eseguito nell’Istituto di Fisica Sperimentale dell’ Universita di Napoli. 
Si ringrazia vivamente il Direttore dell’Istituto, prof. ANTONIO CARRELLI, per 
i mezzi messi a disposizione. 


SUMMARY (*) 


The experimental method is reported emploied in determining the intensity of 


X rays diffracted by liquids and results are given referring to benzene and cyclo- 
hexanol. 


(*) Editor’s translation. 
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Cascades in the Gamma-Ray Spectrum of ’„Po. 


F. DEMICHELIS and R. MALVANO 


Istituto di Fisiea Sperimentale del Politecnico - Torino 


(ricevuto il 25 Gennaio 1953) 


Summary. — Coincidence absorption measurements have been conducted 
on the y-ray spectrum of #44Po. We have proved the existence of a 
strong y-ray cascade of 1,120-1,761 MeV respectively. On the basis of 
our results the two level schemes suggested by E. RUTHERFORD and J. 
SURUGUE have been compared. 


1. — Introduction. 


The study of the decay spectrum of natural radioactive nuclides, initiated 
long before the discovery of artificial radioactivity, has not yet given a know- 
ledge of the Band y spectra of natural nuclides as complete as for nearly all 
the artificial ones. 

The incompleteness of the data about those spectra is not caused by the 
lack of systematic studies on the subject but by the great spectral complexity 
of the nuclides belonging to the natural radioactive series. 

We wished to extend the experimental knowledge of the case of y-ray 
cascades, which are useful in the studies of nuclear levels. Taking also into 
account certain already well know intuitions (1) (?), we have settled our at- 
tention on the case of *4Bi(RaC), of the “Ra series, which presents a very 
complicated decay spectrum. 

It decays (half-life T— 19,7 min) 99,96% into “Po, thus emitting 6 part- 


icles, and 0,04% into "Tl, emitting «-particles (). 


() E. RUTHERFORD, W. B. Lewis and B. V. BowpEN: Proc. Roy. Soc., A 142, 
347 (1933). 

(2) J. SurvuGuE: Journ. de Phys. et le Radium, (8), 7, 145 (1946). 
(2) F. Raserri: Elements of Nuclear Physics (London, 1937), p. an. 
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*4Po, making electromagnetic transitions from the excited to the ground 
state, emits many y-rays of different wave lengths and intensities. 

Our aim in this research is to study the y transitions of ”4Po and in 
particular to verify experimentally the existence of possible y-ray cascades. 

In order to make the subject clearer we set forth more elaborate data 
about spectroscopy of the y-rays and of long range « particles emitted by 


“Po, both of fundamental importance in our problem. 


2. — y-ray spectrum of Po. 


The information usually given concerns the energies and the number of 
quanta emitted by Ra, and its equilibrium decay products; in this case for 
each « particle emitted by “Ra we have in the average a disintegration of 
every nuclide of the “Ra series. 

The data of Table I, obtained from the results given in various papers (*) 
have been acquired experimentally in different ways, i.e. measuring, by means 
of a ßB spectrograph: 

— the energy of the Compton electrons produced in the interaction of 
y-rays with matter and emitted at a certain angle with the direction of the 
incident y-ray; 

— the energy of the internal conversion positrons which originate in 
the interaction of y-rays, having energy > 1 MeV, with the Coulomb field 
near the nucleus; 

— the energy of the internal conversion electrons, which are pulled off 
from the electronic cloud quite close to the nucleus due to the interaction 
of the y-rays. 

The above results have been partly confirmed and improved by M. MLA- 
DENOVIC and A. HEDGRAM (°) in very recent measurements; incidentally we 
may observe that with the procedure and the improvement of the measur- 
ements, the y-ray spectrum of ®;Ra and its equilibrium decay products become 
more and more complicated, because new energy rays are discovered, fortu- 
nately of very small intensity. 

We can state, therefore, that the essential feature of this spectrum is now 
well defined. 


(4) C. D. Ezzrs and G. H. Aston: Proc. Roy. Soc., A 129, 180 (1930); C. D. Eruis: 
Proc. Roy. Soc., A 143, 350 (1935); G. D. LaATysnev: Rev. Mod. Phys., 19, 132 (1947); 
W. Borne: Zeits. f. Phys., 96, 607 (1935). 

() M. MrapeNovic and A. HEDGRAM: Proceedings of the Conference on ß and + Radio- 
activity (Amsterdam, 1952); not yed published. 
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The same cannot be said about the intensity of the y radiations. 


407 


AS 


intensity we mean the average quanta per « particle of “Ra. 

The inaccuracy about the y-rays intensities is due partly to purely expe- 
rimental reasons, which could be avoided by using very strong radioactive 
sources, and partly to uncertainties of the theories which tend to explain the 


interaction phenomena. 


TABLE I. — Gamma-ray spectrum of radium and its equilibrium decay products. 


| 
Average quanta 


| Energy per photon 


| man per alpha ray of 26Ra (MeV) | 
(a) (b) (a) (b) 
226 Ra —> ?%Em 0,01 | — 0,184 = | 
= Te | 
PP pi 0,11 = 0,241 = 
| 0,26 = 0,294 = 
0,45 == | 0,350 - | 
Bi — 24Po 0,03 = | 0,426 = 
— = 0,498 = 
0,66 = 0,607 0,606 
0,06 = 0,766 — 
0,07 = 0,933 = 
0,21 0,20 1,120 1,120 
0,06 0,05 1,238 1,234 
: f 1,290 
0,06 “5 1,379 11.370 
= = = 1,414 (*) 
| — 0,06 | = 1,520 
| 1,620 
0,29 | 1,690 
| 0,26 | 1,761 | 1,761 
= 0,04 == 1,820 
| 2,090 
| 0,07 Nu 2,198 eee 
| = 0,04 = 2,420 
| (c) | (c) 
| | 0,202 
219Po — 2%Pb 7:10 | 0,798 
1,068 


Roy. Soc., A 143, 350 (1934). 


(ec) W. Borne: Zeits. f. Phys., 96, 607 (1935). 


place through direct interaction. 


(b) G. D. LATYSHEV: Rev. Mod. Phys., 19, 132 (1947). 


(a) ©. D. Eruis and G. H. Aston: Proc. Roy. Soc., A 129, 180 (1930); ©. D. ELLIS: Proc. 


| : (*) In the spectrum of the internal conversion electrons there is a line corresponding to | 
a 1,414 MeV photon that however is not found in the Y-ray spectrum. This is a well-established | 
case of non-radiative nuclear transition where the transference of energy to the electrons takes 
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We can thus understand why, in Table I, there are such discordant values 
for the intensities according to the different experimental methods and the 
different experimenters. 

On the other hand it must be emphasized that the precise knowledge of 
the various y-rays lines is very important in order to get the internal con- 
version (°) and the internal pair formation coefficients (*), and thus to establish 
the type of transition (dipole, quadrupole...,) involved in the y-ray emission. 


3. — *4Po long range « particles spectrum. 

Po, emitting « particles, decays into *$Pb. The « particles spectrum 
of Po has been accurately determined by E. RUTHERFORD (!) by range 
measurements as well as with a magnetic spectrograph method. 

The disintegration energy, i.e. the kinetic energy of the « particle plus 
the recoil energy of the nucleus, the energy excess of the long range « particles 
above the energy of the normal ones, and the number of « particles per di- 
sintegration (1) are given in Table II. 


TABLE 11. 
Energy | Energy excess Number 
| of disintegration | of excited level of alpha-particles 
| (corrected for recoil) above normal level per disintegration 
(MeV) | (MeV) x 10-8 
| | 
7,829 | 0 105 
8,437 | 0,608 | 0,43 
| 9,112 | 1,283 (0,45) 
| 9,241 | 1,412 | 22 
| 9,493 | 1,663 | 0,38 
| 9,673 | 1,844 | 1,35 
| 9,844 2,015 | 0,35 
| 9,968 | 2,138 1,06 
| 10,097 2,268 0,36 
10,269 2,439 1,67 
10,342 2,513 0,38 
10,526 | 2,697 1,12 
10,709 | 2,880 | 0,23 


(6) H. R. HULME: Proc. Roy. Soc., A 138, 643 (1932); H. M. TAyLor and N. F. 
Morr: Proc. Roy. Soc., A 138, 665 (1932); H. R. FowLer: Proc. Roy. Soc., A 129, 
1 (1930). 

(7) J. C. JAEGER and H. R. HULME: Proc. Roy. Soc., A 148, 708 (1935). 
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Comparing Table I with Table II it can immediately be deduced that the 
observed long range « particles 
are emitted by “Po formed in 3 
an excited state by ß decay from 
its parent *Bi nuclide. 

From this standpoint the 
energy of these omogeneous 
groups of « particles represents 
the energy of the excited nuclear 
levels. In an extremely short in- 
terval (< 10-11 s) the excited nu- 
cleus decays into the ground state 
and the energy excess of the « 
particles is emitted in the form 
of electromagnetic radiation, 
.either as a single photon or as 
successive photons, in which case 
the excited nucleus passes through Fig. 1. - The solid line levels correspond to 

measured long range «particles; those drawn 


one or more transition energy in dotted and dashed lines have been assumed 
levels to exist according to Rutherford’s scheme. 


2,697 


2,198 
2,142 


1,761 
1, 667 


0,766 
0,608 


0,426 


(0) 


Obviously the number of « par- 
ticles per disintegration is not directly connected with the intensity of the y-rays 
emitted by “Po. On the contrary 
mess vaasıs nero the above two processes compet 
ev Fr) Sree eee e with each other; and while the 
emission probability of long range 
x particles increases with the 
energy of the « particles, and de- 
pends from the angular momen- 
tum of the excited state, the 
probability of emission of the as- 
sociated y-ray is connected with 
the absolute magnitude of dipole 
or quadrupole electrie or ma- 
enetic moment of the transition: 
j.e. it depends on the distribution 
of charges inside the nucleus 


0 
Fig. 2. — The solid line levels correspond to itself (*). 
measured long range « particles; those drawn Knowing the energies and the 
in dotted and dashed lines have been assumed ; ee 3 
to exist according to Surugue’s scheme. intensities of the y-rays and the 


(8) G. Gamow and C. L. CRITCHFIELD: Theory of Atomic Nucleus and Nuclear 
Energy-Sources (Oxford, 1949). 
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energies and the number of « particles per disintegration of “Po, it is possible 
to build up schemes (figg. 1 and 2) (1) (?) of energy levels: these energy level 
schemes account fairly well for almost all experimental data we have 
about the « particles and y-rays energy. 

Building up these schemes it was necessary to consider not only direct 
transitions but more complicated ones between levels where the nucleus looses 
its excess energy in successive stages; i.e. it was necessary to admit, without 
experimental proof, the existence of y-rays cascades, and one or more trans- 
ition levels not connected with measurable long range « particles. 

According to the purpose of our research (§ 1), we tried to verify experi- 
mentally the existence of such cascades (°). In order to test which one of the above 
schemes were the most acceptable we 
measured by an absorption method (2°) 
the energy of the corresponding y-rays. 

This was made possible thanks 
to the new very efficient technique 
of y-ray detection. 


4. — Experimental apparatus. 


The experimental apparatus used 
in our research is suitable to deter- 
mine the absorption curves of two 
coincident y-rays giving rise, na- 
mely, to a cascade of °4Po. 

F, and F, (fig. 3) are two scin- 
tillation detectors consisting in 931A 
photomultipliers, placed so as 
to «see» the scintillations eventually produced inside Nal(Tl) crystals (*). 

The voltage pulses at the output of the photomultipliers F, and F, are 
conveniently amplified through A, and A, and then brought to the grids 
(input) of two asymmetrical univibrators U, and U, acting as pulse shapers. 

After a process of differentiation, the pulses are applied to the grids of 
a twin triode 7,, which indicates the possible coincidence of two pulses at 
its cathode output. 

The numerator N records the total number of coincidences; ©, and C, are 


(°) F. DEMICHELIS and R. Matyano: Nuovo Cimento, 9, 1106 (1952). 

(1) A. C. G. MITCcHELL: Rev. Mod. Phys., 20, 296 (1948). 

(*) Provided by the Physikalisches Synthetische Einkristalle Laboratorium - Kiel 
(Germany). 
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two counting rate meter (C.R.M.), which indicate the total counting rate of 
pulses at the output of the pulse shapers. 

The d.c. voltage for the photomultipliers is given by the power supply A:T 
and measured by means of the electronic voltmeter V with an accuracy of 
0,5% (or better). 

The photomultiplier voltage is automatically regulated by R, and R, so 
that, for. the same geometrical and physical conditions of the source and the 
detectors, the counting rate of the pulses is constant altough other concurring 
conditions (temperature, voltages on the single dynodes,...) may noticeably 
change (11). 

The total voltage at the ends of the photomultipliers is placed at a voltage 
as high as 940 volt in order that the resolving time of the coincidence appa- 
ratus be indipendent from the pulse heights of the photomultipliers ('?). 

The absorbers a, and a, consist of a set of lead disks of various thickness 
and the same diameter of about 2 cm; they are placed between the source 
and the scintillation counters and entirely shield the one from the others. 

The solid angles formed by the absorbers with the source and with the 
counters are very wide; but this fact cannot give rise to any inconvenience 
because accidental coincidence with Compton photons produced inside the 
absorbers are very improbable, as one can see by simple computation. 

The source S (radium bromide in equilibrium with its decay products) has 
an activity I of about 35 uc and is contained inside a small glass enclosure. 

From Table I we deduce that “Ra, “Em, “Pb, y-rays have too low 
energies; while the *SPb y-rays have too small an intensity to sensibly in- 
fluence our experimental results, which therefore exclusively concern the *{iPo. 


5. — Measurements and experimental results. 


Repeated coincidence measurements have been made with various thick- 
nesss of absorbers. 

A first series of measurements (symmetrical arrangement) have been done 
by placing an equal thickness of absorber in front of each detector and by 
counting the number of the coincidences for each thickness of lead from 
z=%=0to x — 12 mm. 

The counting rate of pulses from each detector being kept constant by 
a suitable stabilizer, as before said, if we introduce between source and de- 
tector an absorber, the detector counting rate would remain constant and the 


(1) R. AscoLı: Nuovo Cimento, 9, 615 (1952). 
(2) F. DemicuELis and R. MAILVANO: Nuovo Cimento, 9, 1227 (1952). 
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photomultiplier voltage, and consequently the deteetor sensitivity, would 


change. 


In order to avoid this inconvenience we act properly on the stabilizer so 
that we read always the same photomultiplier voltage on the electronie volt- 


meter. 


After all to be sure that the detector sensitivity is always constant, we test 


Fig. 4. 


at intervals, the constancy of the rate of 
pulses given by the C.R.M.: C, and C, 
with no absorber in front of the de- 
tectors. 

A second series of measurements 
(asymmetrical arrangement) have been 
done placing a constant thickness 
(2 mm) of lead in front of one detector 
and a variable thickness from 2 = = 2 
to æ — 16 mm in front of the other. 

The stabilizers À, and À, are regu- 
lated in the same way as before men- 
tioned. 

Points of fig. 4(a) and 4(b) repre- 
sent the experimental data respectively 
for the symmetrical and asymmetrical 
arrangements. In both cases the 
abscissae represent the values of 2 — x, 
and the ordinates represent the ratio 
between the effective coincidence count- 
ing rate N, for various absorber thick- 
nesses and the effective coincidence 
counting rate N, for minimum absorber 


thickness. On the other hand the effective coincidence counting rate is given 
by the total coincidence counting rate, minus the accidental coincidence count- 


ing rate N,. 


This number is given by the well know formula 


NENNE 


where N, and N, are the counting rate of pulses read on each of the C.R.M. 
and t is the resolving time of the coincidence circuit. 

We have measured the value of t by means of two different sources of 
radium, each placed so as to radiate only one detector and, in order to be sure 
that the resolving time t be indipendent of the y-ray energy, we have mea- 
sured this value with various absorber thicknesses from 0 to 16 mm of lead. 
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The value thus determined was constantly (1,28 + 0,01):10” s during all 
our measurements. 

As both N, and N, are affected by a probable error which is less than 1%, 
we may determine the accidental coincidence counting rate with an error less 
than 2%. Taking into account the probable error in the determination of 
the total coincidence counting rate, the pro- 
bable errors which affects the experimental 
points are represented in figg. 4(a) and (b). 

Another series of measurements concern. 
the determination of the effective coinci- 
dence counting rates against the resolving 
time Tr of the coincidence apparatus. In 
front of each detector, 2 mm lead absor- 
bers were constantly placed in order to 
eliminate completely the spurious coinci- 
dences caused by Compton photons. 

Fig. 5 represents the diagram of the ef- 
fective coincidence counting rate N and the 
diagram of the ratio A between the acci- 
dental coincidence counting rate against 
the resolving time 7 (*). In the diagram 
are also represented the probable errors which affect the experimental points. 


6. — Discussion of the experimental results. 


It is well to discuss how to deduce the y-ray energy of the cascades from 
the experimental data and how to calculate their total intensity. 

As to what concerns the energy of the cascades we can easily see that when 
a radioactive nuclide decays to the ground state emitting two y-rays in cascade, 
the decreasing of the effective coincidence counting rate versus the absorber 
thickness is given by the following formula 


N,= N, exp [— (m + u) — %)] 5 


in the symmetrical arrangement and by the formula 


exp [— pu(@ — %)] + exp [— WT o) | 
’ 


| 
Ï 
2 
a 


N,=N,, 


in the asymmetrical arrangement. 


. (*) The ratio A is conveniently corrected for spurious countings of the photo- 
multiplier [see formulae (4), (5), (6)]- 
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N,is the coincidence counting rate for lead thickness # and N. is the 
coincidence counting rate for lead thickness a; mu, and w# are respectively 
the absorption coefficients of the two rays of the cascade. 

On the contrary the preceeding relations become respectively 


(1) N, = Ns, > pioios, EXP [— (ts, + m) — Ho] , 
> L exp = u isrete ilar exp [— 14,(& — %) | 
@) MN. > pile, Ep, 


if the nuclide decays giving rise to more than one cascade. 
The summation index 7 extends to all the cascades; the coefficient p; is 
the intensity of the i-th cascade and Qi, 
1 9 : 10 X-x mm 15 is the efficiency of the counter, for the 
| 
| 


energy of j-th y-ray belonging to the 2-th 
cascade. 

However a practically straight dia- 
gram, having as abscissae 2 — x, and as 
ordinates N,/N, represents these for- 
mulae. The slope of represented straight 
lines is directly dependent on the exi- 
sting cascades. 

The solid line of fig. 6(a) fits best the 
experimental points of the fig. 4(a) when 
the absorber thickness exceeds 1 mm. 

In effect for lead thicknessses below 
1 mm we have a very great absorption 
owing to the existence of many spurious 
coincidences of backscattered quanta in- 
side the crystal detectors. 

The dotted and dashed line is obta- 
ined from formula (1) admitting the 
existence of the 5 cascades of Ruther- 
ford’s scheme (fig. 1); the intensities p,; 
(Table I) are those of LATYSHEV (*) the 
values of mi, are obtained from the table of the absorption coefficients of 
y-rays in lead; for the values of 0; see the Appendix. 

In the contrary assuming the existence of only the following 3 cascades 


*) We have used the data of Larysury, being almost certainly more precise. 
( , & yi I 
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of Rutherford’s scheme: 


1) 0,933 MeV, 1,761 MeV; 


? 
2) 1,376 MeV, 0,760 MeV; 
3) 1,120 MeV, 1,761 MeV; 


we have the broken line of fig. 6(a). 

The dotted line is finally obtained by assuming the existence of the cascades 
of fig. 2 scheme (SURUGUE). 

In fig. 6(b) the different lines have a similar meaning as those in fig. 6(a), 
but they are derived from equation (2) valid for the asymmetrical arrangement. 

Examining both figg. 6(a) and (b) we must conclude that the complete 
scheme suggested by RUTHERFORD seems absolutely inacceptable; on the 
contrary the partial scheme consisting only on the three cascades 1), 2) and 3) 
seems to be in accordance with our experiments. But in this case the two 
fairly strong 1,241 and 2,198 MeV y-rays of the two cascades 


4) 1,241 MeV, 0,426 MeV, 


5) 0,499 MeV, 2,198 MeV, 


would not have any place, because their energy do not correspond directly 
to the energy excess of any long range « particles; except if we assume that 
the hypothetis 1,241 and 2,198 levels have an extremely short life. 

The scheme suggested in the work of SURUGUE could be experimentally 
justified because the dotted lines run very near our experimental lines: how- 
ever we want to point out that this scheme allows the 1,761 and 1,120 MeV 
rays to have the same intensities, while this fact does net agree with the 
experimental data of Table I. 

We are dealing now with the total intensity of the cascades. 

The accidental coincidence counting rate N. may be written in the fol- 
lowing way 


(3) N SEN N 
— 2712,21? > 4:0: > 4505 + Malo + > 1:01 Qin, i > qje1 Qn, | ; 


where I is the source disintegration rate, Q, and 9, the solid angles subtended 
by the detectors, n, and n, the pulse rates of the detectors without source, 
diy 9; the intensities of each y-ray and 0,, 0; the detector efficiencies of the 
receivers for the different energies. 
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From (3) we get 
(4) 27N,N,— 2nn — 27 D 40:1Qın, — 27 ¥ 4:01 Om = 
i 3 
== DANONE >) 110: >! 4505 - 
a 2 


All terms on the left of (4) are experimentally known; on the other hand 
the effective coincidence counting rate is given by 


(5) N= ION > P:i0:i0; + 


Thus the ratio À between (4) and (5) becomes 


tI > Gis > 40; 
(6) A=? 
> P:0:0; 


The experimental behaviour of A versus 7 is represented in fig. 5. The 
slope of the asymptote & to the line indicates the values of the ratio 4. 
D 0:0: 9,0; is experimentally known (see the Appendix); hence we can 


derive Ÿ p,0:0; whose value results 0,50 + 0,05. 


o; and 0, being known we have in conclusion for 


2 PiQiQs 
the complete scheme of RUTHERFORD . . . . f 0,56 
the partial scheme of RUTHERFORD . . . . . . 0,53 
theschemeloMSURUGUE Sos PRE ee gg OEE 


We can deduce from the above results the existence of a strong cascade of 
1,761 and 1,120 MeV y-rays respectively and we can consider more plausible 
the partial scheme suggested by Rutherford both for the energies and for the inten- 
sities; anyway by means only of absorption measurements we cannot exclude the 
validity of the scheme suggested by Surugue (*). 


(*) Note added in the proofs. — We regret that the conclusions deduced by A. H. 
WAPSTRA (Academisch Proefschrift (Amsterdan, 1953)) on the basis of 6 —y coinci- 
dence measurements are not in agreement with our results. 
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APPENDIX 


To interpret the experimental results we had to know the dependence of 
the y-ray detection efficiency (*?) from the energy of the incident y-rays. 

In order to resolve this problem we compared the experimental absorption 
curve (line S of fig. 7) of the y-rays 
emitted by our radium source with 
the theoretical absorption curve, 
drawn on the assumption of cons- 
tant efficiency (line 7). 


Line R represents the ratio S/T between the value of the ordinates of 
lines S and 7. Assuming, in the first approximation, that the energy depen- 
dence of the detector efficiency be of the simple type 


E 
(7) 0 = À |do + di nent 


from the experimentally known ratio S/T and the intensities of each y-ray we 
can determine the coefficients «, and ay. 


(3) J. V. Dunwortu: Rev. Se. Instr., 11, 1677 (1940). 
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We obtained in this way the solid line of fig. 8. The coefficient A of the 
equation (7) depends on the geometrical conditions of the receivers. 

The dotted line of fig. 8 represents the theoretical variation of the energy 
absorbed by an Nal(Tl) crystal, versus the energy of the incident photons. 

From the solid line of fig. 8 we have deduced the values of o/A necessary 
to interprete our results (1). 


(1%) F. NoRLING: Phys. Rev. 58, 277, (1940). 


RIASSUNTO 


Sono state fatte misurazioni di assorbimento in coincidenza sullo spettro y del 24Po. 
Abbiamo sperimentalmente confermato l’esistenza di una intensa cascata di raggi y 
di 1,120 ed 1,761 MeV. Sulla base dei risultati ottenuti si è proceduto ad un confronto- 
tra gli schemi di livelli energetici suggeriti rispettivamente da E. RUTHERFORD e 
J. SURUGUE. 
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Application of second quantization methods 
to the classical statistical mechanics (II). 
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Centre de Physique Nucléaire de l'Université Libre de Bruxelles 


(ricevuto il 28 Gennaio 1953) 


Summary. — The classical mechanics of indistinguishable particles 
discussed in I is further developed. The mechanical foundations of thermo- 
dynamics are discussed with the classical wave theory and the « quantized » 
fields in phase space. It is shown that the formalism eliminates automa- 
tically the Gibbs paradox. The classical statistics is treated with the 
Gibbs-von Neumann method and also by a generalized form of the Boltz- 
mann method which allows to take into account the interactions. The 
introduction of finite cells in phase space leads to formulas similar to 
those of the quantal statistics for free particles, but no more so for inter- 
acting particles. The statistical treatment of the « quantized » field in 
phase space leads immediately to the canonical grand ensemble. It is 
shown that the theory of the « quantized » fields in phase space can be 
derived from the theory of a «non-quantized » field in phase space, in 
which there are only continuous distributions of matter, by a procedure 
of «quantization » altogether similar to that of the quantum theory of 
fields, but not involving the Planck constant, or any other universal 
constant. The «non-quantized » field in phase space corresponds to. a 
theory less accurate than the classical mechanics and gives an approx- 
imation of the same kind as the molecular chaos hypothesis of the kinetic 
gas theory, it contains both the dynamics and the heat theory of a con- 
tinuous medium. A new derivation of the Boltzmann equation involving 
a special kind of time average is given and a similar equation is established 
for the two particle distribution function. 


1. - Introduction. 


In the preceding part of this paper (1) we have shown that it is possible to 
develop in the classical mechanics a wave formalism in phase space which pre- 
sents many of the features of the quantum wave mechanics. The possibility 


(1) M. SchöngEre: Nuovo Cimento, 9, 1139 (1952). This paper will be indicated by I. 
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of that treatment results from the fact that the theory of the Hamilton 
equations of motion for a dynamical system 


dq, CH Op. eel. (oes 
di op Ch nog SE TOR 


(1) 


is equivalent to that of the first order partial differential equation 


. of Eee © oH = 
2 il le IN = = 
(2) et E 2, = Op, Op. 0q,)’ 


which becomes the Liouville equation of statistical mechanics when f(q, p) is 
a probability density in the phase space of the dynamical system in consi- 
deration. The equation (2) being of the first order, the square of the absolute 
value of any solution O(q, p) is also a solution. Therefore we may introduce 
a wave function ©(g, p), in general complex, such that the probability den- 
sity be the square of its absolute value: 


(3) f(a, p) =| OG, P)|?- 


The consideration of the wave function gives some essential new possibilities, 
because it is not restricted to have only real and positive values, as the pro- 
bability density. Thus, when our dynamical system is formed by n identical 
particles and we denote by Tr, the set of coordinates and momenta of the 
j-th particle, we can impose to the wave function the symmetry condition 


(4) DO von One) 
D indicating any permutation of the indices. This leads to a classical form 


of the Bose statistics. A classical analogue of the Fermi statistics is obtained 
by taking, instead of (4), the following condition 


(5) DP .O(t,, 9) Tn) SE od, 5007) Ep) ’ 


D, and M_ denoting even and odd permutations, respectively. When no 
symmetry conditions are imposed to the wave function, we get the Boltzmann 
statistics. From a pure mathematical point of view, we may say that the 
classical dynamics of a system of identical indistinguishable particles is equi- 
valent to the theory of the symmetrical, or of the anti-symmetrical, solutions 
of the corresponding linear partial differential equation of the first order (2). 
It is convenient to consider sometimes this more abstract correspondence 
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rather than that based on the interpretation of the solution of (2) as a pro- 
bability amplitude, as will be seen. 

The «second quantization » formalism gives a unified treatment of the 
evolution equation (2) and the symmetry conditions (4) or (5). In the case 
of the Bose statistics, the absorption and emission operators w(t), y*(r) are 
characterized by the following commutation rules: 


[ [ot v*(c = de 7‘); 


| Wa, po = [p*(z), Pe )]- = 0. 


(6) 

+ denotes a point in the phase space 2 of a single particle and the bracket 
[A, B]- of two operators A and B their commutator: 

(7) [A, B]- = [4, B] = AB— BA. 


In the case of the Fermi statistics, the commutation rules for the absorption 
and emission operators are: 


(8) [ya ya) = dem), Wo, ya) = Pro), vr): = 0, 
(9) [AD Bl, = AB + BA. 
The differential operator L is replaced by the operator Te 


(10) B= Risin + Fe 


pot? 
(11) Rein = | vr), ylr)ac , 
2 
à od a ed 3 
2 — — * *(-! ed Sao SEN ER 7 ! ! ; 
(12) Root |» (t)p QE op He ap! w(t! p(t) Ar At 
Q 


@(x— x!) denoting the potential energy of a couple of particles and m the 
mass of the constituent particles. The partial differential equation (2) is re- 
placed by the following equation for the wave functional x: 


(13) DR. 


We showed in the preceding paper that the « second quantization » formalism 
is entirely equivalent to the wave theory of a n-particle system, when the 


28 - Il Nuovo Cimento. 
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wave functional y is an eigenfunctional of the operator for the number of 
particles N,, 


(14) Nop = | vtteww(e)ae, 


corresponding to the eigenvalue n 
(15) Non. 


The general solutions of (13) do not correspond to determined numbers of 
particles and may be considered as describing assemblies of systems of identical 
particles, each system of the assembly having a different number of particles. 
There is no difficulty whatsoever in extending the «second quantization» for- 
malism to the case of indistinguishable composite constituenss. As a matter 
of fact, the «second quantization » formalism can even be developed when f 
is a function of any number of points 7 in an abstract space ©, in which a 
measure can be introduced suitably, and has an evolution equation of the form 


(16) Aa = Kf, 


K being a linear operator. This general «second quantization » formalism 
will be discussed in a forthcoming paper. 

The main objective of the present paper is the application of the classical 
wave mechanics in phase space and of the corresponding «second quantization » 
formalism to the statistical thermodynamics. As a consequence of the simi- 
larity between the classical and quantal wave mechanics, the von Neumann 
density operators can be used also in the classical wave mechanics. In section 3 
the main results of the quantal form of the Gibbs theory due to von NEU- 
MANN (?), DIRAC (?) and KLEIN (*) are recalled and put under a form adequate 
to our purpose. The canonical and microcanonical ensemble are discussed 
in sections 4 and 5. It is shown that the thermodynamical behaviour of the 
systems described by the symmetrical and anti-symmetrical solutions of the 
Liouville equation is the same in equilibrium. It is however possible to obtain 
a difference in their thermodynamical behaviour in equilibrium by associating 


(2) J. von NEUMANN: Got. Nach., 245 (1927). 

(5) P. A. M. Dirac: Proc. Cam. Phil. Soc., 25, 62 (1929); 26, 376 (1930); 27, 
240 (1930). 

(4) O. KLEIN: Zeits. f. Phys., 72, 767 (1931). 
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to the von Neumann density matrix theory a suitable division of the phase 
space of a single particle into finite cells, i.e. by a suitable modification of the 
density operators describing canonical and microcanonical ensembles. Such 
a division is not quite consistent with a purely classical treatment, but has 
the advantage of giving a simple approximation to the quantal theory. It is 
shown that in the case of systems of non interacting particles the results are 
the same as in the quantal theory, with a suitable choice of the cell size, but 
no more so for interacting particles. 

The density operator method is applied to the statistical treatment of the 
« quantized » field in phase space in section 9. The present analysis differs 
in an essential point from the quantal statistics of the fields developed by 
SCHEIDEGGER and McKay (°), by taking into account the special kind of non 
ergodicity of the field arising from the existence of the integral of the motion 
N,,, which was not considered by those authors. The existence of the 
integral of the motion N,,, requires a modification of the definition 
of the density operator corresponding to the canonical ensemble of 
fields. The modification of the definition of a canonical ensemble of fields 
associated with the integral N,, has the remarkable consequence of leading 
automatically to canonical grand ensembles of particles. The canonical 
grand ensemble thus appears as of a more fundamental nature than the 
canonical ensemble. Another important difference with respect to the work 
of SCHEIDEGGER and McKay lies in the fact that a potential interaction 
between the particles is included in our treatment. 

The main result of the introduction of the indistinguishability in the 
classical statistical thermodynamics is the automatic elimination of the Gibbs 
paradox. It is shown that the Gibbs paradox subsists when the classical 
wave mechanics is used as the basis of the statistical thermodynamics 
without taking into account the indistinguishability of the particles, as 
happens in the quantal statistical mechanics. This result is of philosophical 
interest because it shows that a limitation of the individuality of the particles 
is required, even in a consistent classical picture, in order to avoid the Gibbs 
paradox. 

The classical wave mechanics affords the basis for a very general extension 
of the Boltzmann treatment of statistical mechanies based on the division 
of the phase space into cells. A similar technique can also be applied to the 
« quantized » fields in phase space, as shown in section 8. It is shown in section 7 
that the application of the generalized Boltzmann treatment to the « quantized » 
field in phase space leads also to a kind of canonical grand ensemble of particles. 


(5) A. E. SCHEIDEGGER and C. D. McKay: Phys. Rev., 83, 125 (1951). 


424 M. SCHÖNBERG 


The generalized Boltzmann treatment can also be applied to systems of part- 
icles in interaction and gives an alternative to the density matrix treatment. 

There seems to be a considerable difference between the conceptions of 
statistical mechanics used in I and in the present paper. In I the wave 
functions in phase space ©, or the wave functionals y, were used as the tools 
of the statistical mechanics, whereas now the density operators are employed. 
The difference results from the existence of two essentially different appli- 
cations of the Gibbsian ensembles: 

(a) The use of a probability distribution in the phase space as a means 
for the analysis of the motions of systems whose initial states are incompletely 
specified. 

(b) The use’ of ensembles of similar systems to obtain a mechanical 
foundation for thermodynamics. 

The essential difference between (a) and (b) is not clearly seen when pro- 
bability distributions are employed to describe the ensembles. But it is no 
more so in the classical wave mechanics. In I we had mainly in view the 
application (a), whereas now we deal mostly with (b). Strictly speaking, the 
difference between (a) and (b) may still be overlooked to some extent in the 
unrealistic case of distinguishable identical particles, but the difference must 
be always taken into account in the theory of indistinguishable particles, as 
shown in section 6. It is shown in sections 4 and 5 that the density operators 
for the fundamental canonical and microcanonical ensembles are actually formed 
with mixtures of delta-like wave functions which describe completely specified 
states, i.e. completely known values of the coordinates and momenta of the 
particles. The reason why a classical wave formalism separates the two 
problems (a) and (b) can be understood by taking into account the assumption 
of the randomness of phases underlying the use of the von Neumann density 
operators (see for the discussion of this important point in the quantum me- 
chanics the § 84 in The Principles of Statistical Mechanics of TOLMAN, Oxford, 
1946). There is an important point to be noticed when systems of indistin- 
guishable particles are considered: 

The wave function O(,; 7) describes the amount of information regarding the 
positions and velocities of the particles at the initial time t,, but it does also define 
their statistical nature (bosons or fermions), whereas the probability density 
| to; t)|* describes only the amount of information regarding positions and velo- 
cities. Even when the coordinates and momenta are given exactly at the initial 
time lo, it is still necessary to consider a wave function, in order to specify its 
behaviour under permutations of the particles. 

The necessity of the wave functions, even when the initial positions and 
velocities of the particles are assumed to be exactly known, is a consequence 
of the central part played by the groups of permutations in the theory of 
indistinguishable particles. The way in which the wave function comes in in 
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the description of the motion is discussed in section 2. The imposition of 
a symmetry condition leads to an identification of a number of different 
trajectories of a system in its phase space, and also to the exclusion of certain 
trajectories in the case of fermions. The situation is clarified by the com- 
parison of the usual (incomplete in the case of indistinguishable particles) 
classical description of the motion by the Hamilton equations and the complete 
description given by the adequate solution of the Liouville equation with their 
quantal analogues: the Heisenberg and Schrödinger representations plus 
symmetry condition for the wave function. The classical Hamilton equations 
correspond to the g-number Hamilton equations of the quantal Heisenberg 
representation, which do not include a specification of the statistics of the 
particles. A 

It is obvious that the solution of the Hamilton equations allows to get 
the value of the classical wave function at any time in terms of the wave 
function at the initial time. But the knowledge of the general symmetrical 
or anti-symmetrical solution of the Liouville equation is not sufficient to get 
the general solution of the Hamilton equations. From this point of view the 
classical mechanics of systems of indistinguishable particles appears as simpler 
than the usual classical mechanics of distinguishable particles. This simpli- 
fication has an important bearing in the discussion of the Einstein and Zermelo 
cycles, as well as in the formulation of the ergodic problem, as shown in 
section 2. 

The introduction of the classical wave functions clarifies the meaning of 
the unitary transformations in Hilbert space associated with the motion of 
a classical system, which were introduced by Koopman (°). The Koopman 
Hilbert-space can be taken as that of the classical wave functions. The wave 
functions of the systems of indistinguishable particles correspond to closed 
sub-spaces of the Koopman Hilbert-space. 

The importance of the introduction of the indistinguishability of the part- 
icles is shown in a particularly striking form by the existence of the « quantized » 
fields in phase space associated with indistinguishable particles. The most 
powerful tool of the classical statistical theory, the grand canonical ensemble, 
finds its most natural justification in the theory of the « quantized » fields. 
Thus it is seen that the field formalism, although apparently more complicated, 
leads often to the simplest solutions, because of its more powerful methods. 
Moreover the field formalism throws light on very difficult problems, such as 
the thermodynamical irreversibility. The work done in recent years in the 
quantum theory of fields has shown that the S matrix of fields in interaction 
is different from the unity matrix, so that there is irreversibility. The S matrix 


(5) B. O. Koopman: Proc. Nat. Acad. US'A,: 17. 315 (1931): 
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formalism can obviously be applied to the theory of the y-field in phase space 
and its S matrix is not the unity matrix, as a consequence of the interaction 
between the particles. Thereby the behaviour of the y-field is irreversible. 
The bearing of this result on the thermodynamical irreversibility will be di- 
seussed in a forthcoming paper. We shall now limit ourselves to the remark 
that the field formalism shows that even the statistical theory of a system 
of a finite number of indistinguishable particles leads to the statistics of a 
system with an infinite number of degrees of freedom. 

The fact that the classical theory of indistinguishable particles leads to 
a field formalism similar to that of the quantum theory of fields — the clas- 
sical particles being « quanta » of the fields in phase space — shows that there 
is the possibility of developing a kind of super-classical theory in which there 
are only continuous distributions. This super-classical theory can be obtained 
by ignoring the non commutability of the y and y* and keeping the same 
equations of motion for the e-numbers y, y* thus obtained. The equations 
for y. and y* can be derived from a variational principle, put under hamil- 
tonian form and « quantized». This « quantization » is quite similar to the 
true quantization of a field in ordinary space, the c-number Poisson brackets 
being replaced by commutators with a suitable numerical factor. Neither 
the Planck constant nor any other physical constant is involved in that 
« quantization », which leads to the formalism of the « quantized » field in 
phase space, i.e. to the exact classical theory of indistinguishable particles. 
The super-classical field formalism gives a kind of Hartree approximation to 
the classical wave mechanics, in which the self-forces are not subtracted 
(see the discussion of the quantal analogue of the super-classical formalism 
in the Appendix of HEISENBERG’s book, Die Physikalischen Prinzipien der 
Quantentheorie, Leipzig, 1930). 

The super-classical theory is discussed in sections 12 and 18. It may be 
used in some cases to get a convenient treatment of the motion of a continuous 
medium. The existence of heat is immediately shown by the fact that the 
kinetic energy of the macroscopic motion is less than that of the microscopic 
motion. A definition of the entropy similar to that of Gibbs can be obviously 
introduced but leads to the same difficulty as in the other cases, the total 
entropy being time independent. This point is quite interesting, because it 
shows clearly that an incomplete specification of the motion is not sufficient 
to get an increasing entropy, as often assumed without any satisfactory 
justification. 

The equation for e,= p*y, derived from the super-classical theory co- 
incides with that of the Vlassov formalism (7) for a system of a large number 


(7) A. A. Vrassov: Journ. Exper. Theor. Phys., 18, 840 (1948). 
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of particles of the same kind with two-body forces. The Vlassov method of 
introducing the temperature in the case of stationary motions leads to the 
same value given by the general dynamical definition of local temperature. 

In order to clarify the reasons of the unsatisfactory behaviour of the 
entropy in the super-classical theory we give in section 13 a new derivation 
of the Boltzmann equation starting from the integro-differential equation 
relating the one and two particle probability distributions obtained from the 
Liouville equation (that equation plays a central part in the work of Yvon (°), 
BORN and GREEN (*) and KIRKWOOD (!%1!). Our derivation of the Boltz- 
mann equation has some analogy with that of KIRKWOOD (1), but is essentially 
different insofar the time averaging is done in such a way that the three 
particle collisions do not come in. This results from the fact that our time 
averaging is not taken along the actual notion of a molecule or pair of mole- 
cules, but rather along a virtual motion of a two molecule system. The time 
averaging along a virtual motion appears as the most natural one, once the 
integro-differential equation relating the one and two particle distribution 
functions is written in a convenient form, thus it is seen that the complication 
of the Kirkwood derivation results from the choice of the method of time 
averaging. Our proof shows clearly that the time variation of the entropy 
is due to the time averaging and not to the incomplete specification of the 
motion of the particles or to the molecular chaos assumption. The continuity 
equation of the super-classical theory is seen to coincide with the form taken 
by the integro-differential equation relating the one and two particle proba- 
bility densities when the hypothesis of molecular chaos is introduced. 

The method of time averaging along a virtual motion can also be applied 
to the integro-differential equation relating the two and three particle distri- 
bution functions and gives a new approach to the theory of transport phe- 
nomena in cases other than the diluted gas. 


2. — The Schrodinger and Heisenberg representations in the classical and quantal 
wave theories. 


Let H, be the hamiltonian of a system of n identical particles of mass m: 


n 


2 
= u Pi 1 
a H, = Sd =— += D(x,— xy). 
(17) om 132 2, (x; v) 


(8) J. Yvon: La théorie statistique des flwides (Paris, 1935). 

©) M. Born and H. $. GREEN: Proc. Roy. Soc., À 188, 10 (1947). 
(2) J. G. Kırkwoon: Journ. Chem. Phys., 14, 180 (1946). 

(1) J. G. Kirxwoop: Journ. Chem. Phys., 15, 72 (1947). 
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We shall denote the Liouville operator Z of our system by L, and its wave 
function by ©,. (The notations in this paper will be the same as in I, as far 
as possible). The unitary operator of the motion is U,(t): 


(18) U,,(t) = exp [— #L,] , 0, (t) = U,(t)0,(0). 


The passage to the Heisenberg representation is a contact transformation 
determined by U,,(t): 


n 


(19) AMG) = UA HAM Use), O7 = 7:9, = 9,00); 

A® and A(t) are the operators which correspond to the same time inde- 
pendent physical A quantity in the Schrödinger and Heisenberg represen- 
tations, respectively. The equation of motion in the Heisenberg represen- 
tation is: 


d 


20 ] 
(20) is 


AD (¢) = AS le 
j 1 


Hence we have the Hamilton equations 


er Re a, 
dr pa ? diva 0x2 ? 

because 

| , x . CH, 
(22) [xr, L,10 = ix\(H,, 0), — WH, 2,0) = — 7 (BA, %1),0 = 1 2 (OFF 

L 
and 
(23) Ip, 0 =— Hu, pa = — 1 0, 
CX, 


the Poisson bracket being defined in the usual way: 


; n 0 = > Ss 
(24) (0,%), = > en a) 
VEN 


Ox, Op, epı Ox, 


The operators for the coordinates and momenta are commutable in the Schrö- 
dinger representation of the classical wave formalism, so that the corresponding 
operators at the same time in the Heisenberg representation will also be com- 
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mutable. Moreover, since H, is a function of the # and p, it does commute 
with these operators. As à consequence of (21) we have: 


d CHO (t 
(25) Ë x(t), 220) ler = ) pr) =e 
L 
In a similar way we get: 
d d id 
(26) pro, s | pi. pr) z 5 x(t), ze) = 0. 


Therefore the operators for the coordinates and momenta at two times f 
and t-+ dt do commute, and by successive infinitesimal time displacements 
we see that the same holds for any times ¢ and f. Since the operators ‘at 
the time t= 0 do coincide in the two representations, it follows that the 
v(t) and p(t) are diagonal for all the values of t. The Hamilton equa- 
tions (21) may therefore be replaced by c-number equations, which are essen- 
tially the Hamilton equations of the ordinary form of the classical theory. 

It results from the preceding considerations that, for the functions of the 
coordinates and momenta, the Heisenberg representation attached to the wave 
formalism in phase space is essentially the ordinary form of classical mechanics. 
It is obvious that, in the Heisenberg representation, the symmetry conditions 
(4) or (5) do appear under the form 


(27a) DO” (x, p) = OF (a, P) (Bose particles) 
or 
(27b) DO (a, p) = + %@, P)- (Fermi particles) 


When the coordinates and momenta are given at the time t = 0, we must 
take for bosons 


(28a) OT, Tr) = (05 Tao +.) Ta) = 
es > or 7,0) -- Ott, 0)) = 6™(z, t(0)) , 


n! permut. k 
and for fermions 


(28b) OP (Ts very Tr) = 9,0; Try +++) Tn) = 
Ô(T—T(0)) Ô(T — T0) oan OC ae 7.(0)) | 
Hi Ô(T: — 71(0)) (Ts — T(0)) .... Ô(T: — Ta(0)) = ÿmr, x(0)) 
mn! CODE wees ee a Nar iets 


öl. — 7(0))  O(ta— TD) OF — T(0)) 
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The introduction of the wave functions shows at once a characteristic feature 
of the indistinguishability: The 0” are not, in general, eigenfunctions of the 
operators for the coordinates and momenta Rn and 9)”. The 0” are never- 
theless eigenfunctions of all the operators of the form f(t‘”,. 9), RS) 
being any function of the 7, invariant for permutations of the 7. This results 
from the fact that the maximal theoretically possible information is to know 
that there is a particle at each of the points 7,(0), without distinguishing which 
of the partieles is at each of the points. A set of 6n independent invariant 
functions is the following: 


ik 1...” 
= [ St) = SS Se, Elo) a aD .... 562 — FREE) N ao) 
| j M CET ks 
29) 17 2 
( | eae 
ie = (&) ae (mit (x) — mya)ma) Yo) 
x Ir, ne = = "pop u YO == PO .... po, 
2 5,75, 


The index « denotes a component of a vector. The eigenvalue of the operator 
fe”, .... T°) which corresponds to the eigenfunctions 0° is f(7,(0), ..., T„(0)). 
Although there are 6n independent functions &, 7, it is not sufficient to 
give the eigenvalues &”(0), n°°(0) of the corresponding operators to determine 
the n points 7,(0). Indeed, when the &”(0) are given, the x}(0) are determined, 
but for the order, as roots of the equation 


(30) ut — ET HE. + DE) = 0, 


and a similar conclusion holds for the w?(0) when the &?(0) are given, but no 
information is thus obtained on the way in which the different coordinates 
must be coupled. This shows how awkward it is to work with permutation 
invariants in the phase space wave formalism. 

It is easily seen that the solution of the Hamilton equations allows to solve 
immediately the Liouville equation for the wave function. Indeed, the ge- 
neral solution of the Hamilton equations gives us the coordinates and mo- 
menta of the particles, treated as distinguishable, as functions of the time ¢ 
and of the position and momentum vectors a,, b, at an arbitrarily chosen 
initial time f,: 


| == Ft; a, .., an; bi, On) 


(31) 
| Pilt) = Gilt — ty; ay, -.., Gn; Bi, ..., b,). 


The equation of motion for any function R(t; x, p) is: 


(32) La OB soe ae 
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Hence, when R is a solution of the Liouville equation, we have: 
(33) R(t; F(t—t; a,b), G(t—t,;a,b)) = Rit,; a,b). 


By applying (33) to the wave function ©, we get the solution of the wave 
equation : 


(34) 9,(; x, p) = O,(to; Fu —t; x, p), Gt —t; x, p)). 
Equation (34) shows that U,„(t) is a substitution operator 
(35) U,,(t)O(x, p) = O(F(—t; x, p), G(—t; x, p)), 


which replaces the x, by the F,\(— t; x, p) and the p, by the G,(—t; x, p). 

The equation (33) is not dependent on the symmetry conditions imposed 
to the wave function. Let us examine now the effect of the symmetry con- 
ditions. The equations 


(36) 


describe the motion of the particles, treated as distinguishable, which cor- 
responds to a permutation of the initial coordinates and momenta. If we 
assume that ©, is symmetrical at the time f,, we get: 


Nee alles Ox. mn Ong Bd) = 
n 1 n 1 n 


= Old; Ay, wey Any b, so. b,) == 9,(; x(t), p(t)) . 


Hence the trajectories in the phase space of the system corresponding to 
{x(t), p(t)} and to {x(t), … k 9 Pr} do describe the same state of 
motion of the system of bosons. In the case of fermions we get 


(38) Ont; xx, es PO 
the sign corresponding to the parity of the permutation {Tet Tar soap pS 
conclusion is the same as in the case of the bosons. Hence, when the particles 
are indistinguishable the set of trajectories corresponding to the different 
{at xs Pe x} is lumped into a single state of motion. It is obvious 


by reasons of symmetry that: 


15e. 


(39) ln), Pia, = Pl). 
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Hence the symmetry (or anti-symmetry) of the wave function at one instant 
of time t, ensures the symmetry (or anti-symmetry) for all the values of t, as 
a consequence of (37) (or of (38)). 

The preceding results show clearly that the consideration of the wave 
function is necessary to give a complete dynamical description of a 
system of indistinguishable classical particles, even when there is no uncer- 
tainty on the initial values of the coordinates and momenta of the particles, 
as was emphasized in the Introduction. The situation is similar to that of 
the quantal wave mechanics. The quantal equations of motion of the oper- 
ators for the coordinates and momenta are also Hamilton equations, they 
have exactly the same form of the corresponding Classical equations (21), the 
only difference consisting in the non commutability of the operators for the 
coordinates and momenta. As happens in our classical theory, the Hamilton 
equations of the quantal Heisenberg representation do not specify the stati- 
stical nature of the identical particles. This is done by the choice of the time 
independent wave function of the quantal Heisenberg representation. 

There is one aspect of the dynamical behaviour in whose analysis the in- 
troduction of the indistinguishability is of fundamental importance: the 
quasi-periodicity of the motion on a finite surface of constant energy in phase 
space. The duration of a Poincaré cycle must now be defined in terms of the 
states of motion, instead of trajectories. Since a state of motion includes in 
general n! trajectories in the case of a system of n indistinguishable particles, 
the duration of the true cycle will be shorter than in the theory of distinguis- 
hable particles. Thus the Einstein point of view is seen to be satisfactory 
(see the interesting remarks of BORN: Natural Philosophy of Cause and Chance, 
Oxford, 1949, appendix VI). The necessity of defining the representative 
point in phase space « without regard of the individuality of the particles » 
(BORN) has a rigorous justification in the present theory. 

In the case of systems of indistinguishable particles, the ergodic problem 
is no more the same as in the ordinary classical theory, because now we must 
consider time averages in a state of motion and not along a trajectory of the 
system in phase space. The time averaging in a state of motion canbe split 
into two steps: (a) time averaging along each of the trajeetories pertaining 
to the state of motion; (b) an ordinary average of the time averages along 
the different trajectories. It is obvious that the ergodic condition for systems. 
of indistinguishable particles is less stringent than in the ordinary theory, 
especially in the case of systems with a large number of particles. We shall 
discuss elsewhere in more detail this important point. 

The above results show that the theory of the motion of systems of 
indistinguishable particles is simpler than that of systems of the same structure 
whose particles are distinguishable. From a very general point of view this 
is quite obvious, since in the latter case we must determine all the solutions 
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of the Liouville equation, whereas in the former case only a part of these 
solutions are required. It is well known that a n particle system has only 6n 
independent constants of the motion, hence the corresponding Liouville equation 
has only 6n independent solutions. These can not be all symmetrical in the 
variables of the n particles, because this would imply that all the solutions 
are symmetrical. It is easily seen that the number of the independent anti- 
symmetrical constants of the motion is equal to that of the symmetrical 
constants and thereby less than 6n. The preceding remarks lead to an im- 
portant conclusion: 

The knowledge of all the symmetrical (or anti-symmetrical) wave functions is 
not sufficient to obtain the general solution of the Hamilton equations, wh ich depends 
on the knowledge of 6n independent constants of the motion. 


3. — The von Neumann density matrix. 


The analogue of the Gibbs statistical mechanics in the quantal theory was 
given by voN NEUMANN (?) and DIRAC (8). The description of an incompletely 
known state is given by the density operator R of von Neumann, which 
satisfies the quantal analogue of the Lionville equation for the classical distri- 
bution function 
ae 


ny 
UHR, 
don Pl 


(40) 
H denoting the quantal hamiltonian operator. If it is only known that at the 
initial time f, the system was in one of the orthogonal states of wave funct- 


ions ga, with the respective probabilities q,, the initial density operator must 
be taken as 


(41) Re) = Zu. 
with 
(42) yaa 


the P, being the projection operators corresponding to the states ga. There 
is no loss of generality in assuming that the y, are a complete set of orthogonal 
and normalized functions, since this can always be done by taking some of 
the q, equal to 0. In the representation whose basic vectors are the y,, the 
matrix elements of R(t,) have the simple form: 


(43) <a’ | Rt) | a") = qa'do',a" - 
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In the particular case of a completely specified quantal state W, R is a 
projection operator 


(44) R = Py; 


P,, being the projection operator of the state W. Let us consider the parti- 
eular case of a system of n identical spinless particles and the representation 
in which their coordinates are diagonal. The matrix elements of R in this 
representation are: 


(45) CR A …, %, | R(t |}, [RE HY = DAC AT ICE 


a 


(46) p(t) = exp - . (tt Fa: 


The probability of finding the n particles in the elements of volume dx, at 
the: time t is € 0... Xn P(E) | x4, -..) Xn dx..dx,. 

We shall now consider indistinguishable particles. In order to get the 
Fermi statistics we would have to take particles with half-integer spin. The 
introduction of spin variables would represent an unnecessary complication 
for the goal we have in view. Moreover, since we are dealing wtih a non 
relativistic theory, there is no inconsistency in imposing the Fermi statistics 
to spinless particles. We shall thereby use the simpler formalism of spinless 
particles for the two statistics. 

In the von Neumann quantal statistical mechanics, an entropy is attributed 
to the incompletely specified state described by R: 


(47) SP=— k Tr (Rlog R). 


Tr denotes the trace of the operator. The definition (47) is also applicable 
to a completely specified state and gives to any of such states the entropy 0 
as well known. The trace of an operator can be computed in any represen- 
tation, by using that corresponding to (43) we get: 


(48) Sato) a GAL ge Ms Bip RE 


We shall now assume that the , are eigenfunctions of H, corresponding to 
energy levels E,. In the case of degeneracy we shall have the same value E, 
for g, levels. The expectation value of the energy in the condition described 
by R(t) is: 


(49) (H>,,= Tr (R(t) H Tr? GB 
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In the present case R is time independent, thereby we may drop the time 
variable. The most important ensemble of this type is the canonical ensemble: 


exp - Be 
(50) da De ; = > exp |— 5 ; 
(51) ot | ie 
The parameter Y has the following value: 
(52) WY — — 6 log Z 


It is well known that the canonical ensemble describes the behaviour of 
a system in thermal equilibrium with a bath at constant temperature 7 


(53) = KT, 


k being the Boltzmann constant. The entropy of the system is Scan 


(54) Scan so =k log Z + "can ’ 


(55) <Hycan = > 32 Ey exp |— 


E, 

KT 

Y is therefore the Helmholtz free energy of the system F: 
(56) Vz — hs Z = <A pon TS can =L: 


The formalism of the von Neumann density operators becomes more com- 
plicated when the wave function of the system fulfills symmetry conditions. 
The allowed wave functions of the system form a closed sub-space of the 
Hilbert space of the system without symmetry conditions. Let us denote 


+ 
by P” the projection operator on the closed Hilbert sub-space corresponding 


to the symmetrical wave functions of a system of n particles and by P” the 
projection operator on the closed sub-space of the anti-symmetrical wave 
functions. In the construction of the ensemble the allowed P% must satisfy 
the conditions: 


+ + + + + + + + 
(57a) > Pw = pw ee = 4 P® Fnpn— a (bosons) , 


(676) pw = pn (Eve = ÿ,P® Popn— 2 (fermions). 
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The density operators of the canonical ensembles in the cases of indistinguis- 
hability are 


SH, 
on Ru or 0 | 
with 
4 E, 
(59) PO log Z®, ZY = Yg,exp 1-5 
Eq 


Ya and J being the weights of the energy level #, for bosons and fermions, 
respectively. The difference of the weights renders the free energies different in 


the various quantal statistics. 


4. — Canonical ensembles in the classical wave formalism. 


The introduction of the von Neumann operators in the phase space wave 
mechanics and in the theory of the « quantized » fields in phase space can 
obviously be done by applying the rules of the preceding section. From the 
formal point of view the only differences arise from the necessity of consi- 
dering continuous spectra of eigenvalues. There are nevertheless some very 
important differences resulting from the fact that the classical wave functions 
in phase space describe both completely specified states and incompletely 
specified ones, such as those corresponding to approximately known initial 
values of the coordinates and momenta of the particles. In other words, the 
most significant differences arise from the fact that the pure states of the 
wave mechanics in phase space are not entirely similar to those of the quantal 
wave mechanics. 

The above difficulty indicates that it is necessary to distinguish between 
two different applications of the Gibbs methods: 

a) The use of a probability distribution in phase space as a means for 
the analysis of the motion of systems whose states are incompletely specified. 

b) The use of ensembles of similar systems in order to obtain a mecha- 
nical foundation for thermodynamics. 

It is clear that the problem of the analysis of the motion of a mechanical 
system, whose initial state is incompletely known, is essentially different from 
that of the foundations of thermodynamics. This difference was not em- 
phasized in Gibbs’ book, because of his definition of statistical mechanics as 
the branch of mechanics dealing with the motions of assemblies of non inter- 
acting similar systems. In the Elementary principles in statistical mechanics, 
Gress considered different applications of the methods of statistical mechanics, 
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in the above sense, some of them entirely unrelated to thermodynamics. Later 
some authors have assumed that the problem of the mechanical foundation 
of thermodynamies was essentially a particular case of that of the analysis 
of the motion of a mechanical system whose condition is incompletely spe- 
cified, because the thermodynamical specifieation of the state of a body does 
not include a complete mechanical specification of the condition of its micro- 
scopie parts. This point of view is misleading, since the relevant fact is that 
the thermodynamical laws are not dependent on the complete microscopical 
specification, so that such a specification is not required, even if it would 
be known. 

The distinction of the two problems (a) and (b) in our theory comes out 
in a clear cut way as a consequence of the nature of the wave formalism in 
phase space, in which both completely and incompletely specified conditions 
are described by wave functions and appear as pure states, whereas an en- 
semble in the sense (b) is described by a von Neumann density operator. 
The density operator RX), of a canonical ensemble of n-particle systems, 
in the case of distinguishable particles, is 


(60) Reg = exp - — =] 
Ü op 


as a consequence of (51), or rather of the analogous formula in our classical 
wave theory, with 


fi TELA CE tn 
(61) Ww = — 0 log exp — — | Zr - Tn) durer. 
hence 
Ur al ho ee 


dr) em): 


Te ee Hole ihe} 2] 
ex D u LE Cael Ô 


In the present case the density operator R®, is built with the eigenfunctions 
of the operators x and p (commutable) having continuous spectra. A com- 
plete set of orthogonal and normalized eigenfunctions of the operators for 
the coordinates and momenta of the n-particles is given by the 6(t,— 171) ... 


On in): 
| a Ô(T1 — 71) es O(Tr—- a) oo sm — 7.) Ag ee =) 
(63) 


| joer 6(t%1— er ee OG — Go) = pe” 6(%1— Ta) + Ô(En— Tri 


29 - Il Nuovo Cimento. 
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The normalization of the eigenfunctions corresponds to the well known rule: 


(64) [dr did m). Dome 


7 


= 0(%1— Ti) tee O(T,— ae 
Equation (62) can be written as follows: 


[or 


NS 
can | Ts ve.) Tr? = 


(65) Gage 2) ce 


= Es .dr„exp 7 bli u). Aero 


The right hand side of (65) shows that the mixed state described by RX is formed 

with the eigenstates of the coordinates and momenta 6(t,—1,) ... 6(t,—T,), 

Y„— H(t 

the respective proportions being exp = p (7) 
probability of finding the n particles in the elements of volume d,7 of Q. 

Let us consider now the cases of indistinguishability. The corresponding 


canonical density operators are given by a formula analogous to (58), with 
ap. 


dt,...dt,, hence this is the 


phase space projection operators P”: 


(66) Re, = Pin exp 


It is easily seen that 


(67) a EG ee een 

hence: 

(68) PRE LR MNT OL | 2 = Dr, 2") 
We have: 

(69) Tal tee Tal LRU, Ken 


A ae PAG 2 i] 
= | az, ye AT, EXP pee Ole, 2,0% (7% 


Equation (69) shows that the RD, are built with the wave functions 6°?(z, r) 
of states of completely known coordinates and momenta, which are analogous 


to the pure states of the quantum theory. Since permutations of the 7 in 
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the ö%(r, 7) do not lead to different states, we must use the following normal- 
ization: 


rm £ = 
(70) fer 


Q 


VO — HAUT) 
0 


| dr, ... dr, — 1 . 


In order to get the thermodynamical interpretation we may use the von Neu- 
mann definition of the entropy. The same reasoning of section 3 shows that 
y® is the Helmholtz free energy: 


(71) FD = WY — KT log 


al exp | N GER ccc az, 6 


2 


This formula differs from (61) only by the divisor n! under the logarithm. 
It is well known that such a divisor must be introduced to render the free 
energy an extensive quantity. In the usual form of the classical theory the 
divisor n! is introduced a posteriori, but in our theory this is not necessary 
because the formalism leads automatically to the correct value, as a conse- 
quence of the indistinguishability of the particles. 
The von Neumann definition (47) gives for the entropy of the system of 7 
indistinguishable particles the formula 
k W® — H,(t)| YP—- A, (t 
(72) se=— [exp = LT sea (7) 
2 


A dt... Ota 


account being taken of the normalization (70). In the case of distinguishable 
particles the von Neumann definition gives the ordinary expression of the 
entropy: 


i J 
(73) N —= | exp 


Q 


It is easily seen that in the case of a gas of non interacting particles of mass m 
enclosed in a volume V: 


(74) Jexe|- Ep | Cie n= amer. 


a 


CrmkT)"? Vie 


(75) yo — __ kT log = 


’ 
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and by applying the Stirling formula we get: 


7 


Van 5 
(76) SY = nk |log— +5 log (2mmiT) +3 


With this expression of the entropy there is no Gibbs paradox. 

The results obtained with the von Neumann density operators in our clas- 
sical wave theory are not so satisfactory as in the quantal statistics, because 
the entropy still gains an additive constant term by a change of units, as in 
the ordinary Gibbs theory. This is due to the fact that there is no division 
of the phase space 2 into finite cells. The introduction of finite cells leads 
also to a different thermodynamical behaviour of the systems of classical 
bosons and fermions, as can be expected. 

Let us divide the phase space © into non overlapping cells w, of volume w 
and introduce the functions f(r; w,) which are equal to 1 when 7 lies inside w, 
and to 0 when r lies outside w,. We shall assume that the energy of the 
system varies little when the n points t, move inside the cells of any of the 
Sets (D; ++) Oy - In order to build a density matrix we shall pick up a point 7, 
in each of the w, and take 


ae i à Ye (x 
VER 2 (ny) Lig os 4 N Re BI 
CPE Ty |e lee tes en Ar > exp 5 1 An/|. 
A 


rn / 2 
(77) DS (CT on) DENT any gat Coe ee 


(78) Dit, Tr} Way ++) on) = (n!œ") Ust... SM (ES Oa, )--f (Tr, 0) 


permut. k 


r being the number of the distinct A, ($1, s.,...,s,) their multiplicities and 


| Kt; où) Mt; @,)....Mtı; @a,) 
(79) De, Uns C,3e..s Wa ) = (n ! Orme | It; ©.) fe; O3.) +. à eee ©) 
n . . | 
| 
| obs @,) 1 ; ©)»: : i Ges 3) | 


The 7% „ are determined by the normalization of the probabilities: 


1 ee u: 
(80) en > Gal, ’ An) exp 5 es un 2 Tan) == il ’ 
LOS Vote, 6) 
(81a) OA ee Ads LS, s,! 
al LOL N 
(81d) Ware 
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n 19: (a, An) is simply the number of times that DE eng; Os ces O2) 


appears in the right hand side of (77). It follows from (80) that: 


1 des 
— D Ga (Ar y+) An) EXP on Tan) 
Ai). 


Lye 


(82) exp |- Free 


The sum in the right hand side of (82) is larger in the case of bosons than in 
the case of fermions. When the particles do not interact, 7, can be easily 


evaluated with the saddle-point method (see SCHRÖDINGER: Statistical Thermo- 
dynamics, chapter VII, Cambridge, 1948): 


2 
(83) D, = O|nlog ls + > log (are aE fu Cbg | 2 
7 2 


£, is determined by the equation : 


(84) W = > 


Thus in the case of a gas of non interacting particles we get the well known 
quantal values of the free energy, with the only difference of the weight 
introduced by the spin, provided the volume of the cells w is taken as h®°. 

The division of the phase space Q into finite cells corresponds to the Boltz- 
mann approach to statistical mechanics. The above results show that in the 
case of free particles the Boltzmann finite cell method applied to our classical 
theory of indistinguishable particles gives results similar to those of the quantal 
theory, with the choice w — h?. In more complicated cases the introduction 
of finite cells will not, in general, lead to the quantal results, even when 
O— he. 

By using the density operators R® we can also get formulas for the 
average values of the occupation numbers of the different cells w, which 
are similar to those of the quantal statistics of systems of non interacting 
particles. For instance, in the case of a system of n free bosons the 
probability of finding the set of occupation numbers n, , N, + with Sn =n 

2 


(n) Peon oye 
yO ,— dm pz [2m 
is exp = |. Hence the average value of n, is 


ies — > m pil2m 
7 


/ \ n> [ > ; gd 
(85) LS a Ny exp |= ee : 
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>’ denoting that the condition In, =n is satisfied. It follows from (85) 
4 
and (82) that 


op 1 2 
2 = exp| 4] 1 


In a similar way we get: 


(87) (m= 


The values of the ¢, depend obviously on the choice of o. 


5. — The microcanonical ensembles in the classical wave mechanics. 


The microcanonical distribution is the most fundamental tool of the classical 
statistical mechanics because it describes the thermodynamical behaviour of 
an isolated system. The average values taken in a microcanonical distribution 
are also important because of the theorem of ORNSTEIN (1?) and KRoù (1%) 
which shows that the average value of a quantity in the microcanonical di- 
stribution will coincide with its time average whenever the time average be 
the same for all the trajectories lying on a surface of constant energy in phase 
space (see the very interesting discussion of the microcanonical distribution 
given by P. HERTZ in Statistische Mechanik, vol. I, part 2, of the WEBER-GANS 
Repertorium der Physik, Leipzig, 1916). 

The matrix elements of the density operator R/,. of the microcanonical 
ensemble corresponding to an energy interval H— E + AE are: 


7e. / | (n) I! " ree 
(88) Tyye zal Renio tiara Te 
E+AE 


== Cm [ dt ... dr, O(t; = ia) ... öT, = Tn) 6(t;—7,) ... == Tn) . 
E 


+ 
In the cases of indistinguishability RY, is replaced by Rihe: 


E+AE 
(89) RA or far, de, dr", 7) 2(r/, 7) . 
E 


(2) L. S. ORNSTEIN: K. Akad. Wetens. Amsterdam, 818 (1911). 
(5) J. Krod: Bul. Acad. Sci. Cracovie, 563 (1913). 
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The C are determined by the relations: 


E+dE 


(90) ca fa ar, — 1, 
EB 
E+AE 
1 2 2 
(91) 7 or | a, ade: 10: 
E 


The von Neumann definition of the entropy gives: 


E+dE 
(92) Sm. == — KO log 0% | Gad ee, 
{ 
E+AE 
(93) SRE C® log or | dr, ...dt, =— k log C2. 
E 


It is easily seen that in the case of non interacting particles: 


E+AE 
= a 3 3 AEH 
(94) log fa ee) os V + log (27zmE)? — > (to sa —- i) + log a 


£ 


Hence we have: 


Val Die 3 AB 
(95) Se uk los (2rmE)?? + log — + 2 a | | k log E 


bo 
Le 
t 


The temperature can be determined by the well known thermodynamical 
formula: 


(96) (3) = = 


The last term in the right hand side of (95) may be neglected and we get: 


(97) E = ra I 
Ve, 5 
(98) De = nk log (2rmkT)®? Ze log n SF 9 


This is the same expression obtained with the canonical ensemble. 
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In order to get a different thermodynamical behaviour in equilibrium for 
the classical bosons and fermions it is necessary to introduce finite cells in 2, 
as in the case of the canonical ensembles. It is worthwhile to remark that 
the density operators of the microcanonical distributions are built with 
6(t,— Ta). O(Tr— Tn) or 6%(t, T) which describe states analogous to the pure 
states of quantum theory, as happens also in the case of the canonical en- 
sembles. The finite cells can be introduced as follows: 


+ 
ie / (n) Im Re (n) y R 
Ti Tal Log ici ee, Un) = CSE D G(H,(T,, 5 oreg a E) 


Ann 
(99) > Dm, aie Tas Op, + ©.) Dit OO} Tu (OZ EEE) ©) ’ 

ee Hage BAR 
(100) G(æ —E)= | 

| 0 ana EN Ame 


The entropies of the ensembles (99) are: 


k rs 

= — ai Ce > AGE .. An)G( (Ta, » erty T,,) — E) log 1070 (He Tan) — E)] 
k 

= — OP, log Pa I Ga (day ces Ela )— B) = — blog CP 
2 ERW 


In the derivation of (101) we used the normalization condition: 


(102) SOD Gane NE nee) D) le 


The use of the ensembles R,, „.. is obviously related to the ordinary elementary 
method of the energy states, in which a macroscopical state of energy B is 
defined by the numbers of particles in the different one particle quantal states. 
This method of the energy states appears rather unsatisfactory and has been 
criticized (see P. JORDAN: Statistische Mechanik, Vorwort, Braunschweig, 1944). 
Our treatment shows clearly that the so-called energy state is simply a mixture 
described by a suitable density operator. From this point of view the special 
enumeration procedures for bosons and fermions appear entirely justified, 
since they follow immediately from the structure of the density operators. 
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6. — Inadequacy of the elassical wave functions for the definition of the entropy. 


The results of the two preceding sections allow us to understand why the 
wave functions are inadequate for the definition of the entropy. Let us con- 
sider a normalized wave function 0,(t, .…, Tn) 


(103) | On (T9 -+-) Tn) | Ta «ee dtr = 1, 


Q 


describing a condition of a n-particle system in which the amount of inform- 
ation is less than the maximal possible. In the case of distinguishable particles 
we may use @,, to describe an ensemble whose entropy is given by the Gibbs 
formula: 


(104) S,=-—-k | On try Zn)” log | Otis +++) Tn) |? TI IT 
Q 


For instance, a canonical ensemble can be obtained by taking: 


HT HG et) 
mad : À 


(105) | On (Try +++) Tn) P= x 


The situation is no more the same when the particles are indistinguishable. 
Thus, in the case of fermions, it is not possible to find any wave function such 
that 


RESTE ou 
(106) ka, Tr) P= exp | en | | 
because the left hand side vanishes when two of the points 7, do coincide in 2. 
In the same way it is seen that it is not possible to describe a microcanonical 
ensemble for fermions by a wave function 0)”. 
In the case of bosons, it is possible to take 
EIS H(t, Cor) Tn) 


(107) OO (Ts DICH) am) = exp = 20 ’ 


but the normalization (103) does not lead to the correct value of 7”. If we 
try to avoid this difficulty, by taking quite arbitrarily 

1 (+) {2 
(108) MO (rs...) desde; 


n! 
Le] 
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instead of the normalization (103), we will get trouble from (104), since 


(109) —k [ We a OL? a 


Le] 


RE 7 7 > 71 dt, ... dr, , 


=k J exp oe mean nn 55 Ca) 


is not the correct value of the entropy. It would be necessary to divide the 
right hand side of (104) by n! to get the correct expression for the entropy 
and this is also not satisfactory. 

The possibility of introducing ensembles in two essentially different ways 
is an advantage of the present formalism, because it allows to distinguish the 
two totally different applications of the ensembles denoted by (a) and (b) in 
the discussion at the beginning of section 4 Whereas (a) may be considered 
as a rather subjective matter, (b) is obviously an entirely objective problem 
which can not depend on the limitations of our information. 

The above considerations lead to an important question: What is the nature 
of a condition of a mechanical system described by a density operator? It seems 
satisfactory to assume that a density operator, which is not equivalent to a wave 
function, describes the condition of a mechanical system submitted to the action 
of very small and irregular external influences not included in its hamiltonian. 
The importance of such small irregular external actions was already shown 
by BOREL (11), a long time ago, and has again been emphasized more recently 
in the Statistical Thermodynamics of FOWLER and GUGGENHEIM (§ 102 of the 
Introduction, Cambridge, 1952). Boren pointed out that the effect of the 
small irregular perturbations, which do always exists, will very rapidly put 
a closed system into a condition corresponding to a microcanonical ensemble. 

As a matter of fact, it is possible to combine the two kinds on ensembles 
and to consider density operators built with wave functions describing incom- 


pletely specified states, such as the RiP and Ric, discussed in sections 4 
and 5. It is not altogether clear whether the consideration of these more 
complicated ensembles, which combine the methods of Boltzmann and Gibbs, 
is strictly necessary in a purely classical theory. These more complicated 
ensembles may be of practical convenience to get approximations to the 


quantal theory, because of the finite cells introduced by ‘quantization. 


(4) E. BOREL: Introduction géométrique à quelques théories physiques (Paris, 1914); 
Ene. Sei. Math., tome IV-2, suppl. 2 (Paris, 1915). 
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7. — Statistical treatment of the fields in phase space along Boltzmann lines. 


The «quantized » fields in phase space give the most satisfactory and 
complete theory of the systems of indistinguishable particles from the purely 
mechanical point of view. It will be seen that the same is true also from a 
thermodynamical point of view. The von Neumann density operators allow 
us to describe the microcanonical and canonical ensembles in the wave me- 
chanics of the phase space, we shall see that the field formalism allows us to 
get the grand ensembles, i.e. the most general and powerful tool of the sta- 
tistical mechanics. 

It is remarkable that the theory of the classical « quantized » fields in 
phase space leads very naturally to the consideration of finite cells in the 
phase space ©. Let us introduce the operators N (wm), already discussed in I: 


(110) No) = | y*(r)y(r) dr. 


oO 


It was shown in I that the N(») have the eigenvalues 0, 1, 2, ... (co) and that 
(111) [N(w), N(w')] = 0, 


for any pair of volumes w and w’. The group of all the N(w) is a complete 
set of commutable operators for the field. Indeed, the wave functional of 
the field y may be expanded in a series of the following type (Fock expansion) 


foo} 


1 
(112) % = x + > iat 


| OA Ta) y, Ta) Ar, den, 
o 

4, being the wave functional of the vacuum and 

(113) (Try =, Tn) = Y* (tH) WERK » 


as shown in section 5 of I. When y, is normalized, KART Tu) Oty can AT 
gives the probability of finding n particles in the field and located in the ele- 
ments of volume dr, containing the corresponding 7,. The wave functionals 
are therefore determined by the probability amplitudes for the phase space 
distributions of any number of particles and this shows that the group of all 
the N(w) is a complete set of commutable operators. 

In the representation in which all the N (») are diagonal, the wave functional 
of the field is a functional of a domain function N'(w) taking the values 
0,1, 2,... (co) for any volume w: 7[N’(o)], \z[N (o)]|? being proportional 
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to the probability of finding N'(w) particles in the volume wm. Therefore the 
« quantized » field formalism leads naturally to the consideration of finite vo- 
lumes of m and their occupation numbers. The occupation numbers N’ take 
however all the positive integer values both for fermions and bosons. The 
operator for the energy of the field is #: 


s if 
HU = Be = vH) y(t) dr + sl p*(t)p*(t' )O(x — x')y(t')p(t)drdr’ . 


2 Q 
The first term in the right hand side of (114) corresponds to the kinetic energy 


and the second one to the potential energy. We shall denote them by Hin 
and #,,,,° When @(0) is finite, we may write: 


When the y(t) and y*(r) are expanded in series of orthogonal and normalized 
functions 


(116) p(t) = > a;u,(r), w(t) => aux (T), 
7 à 
the u,(t) being the functions of a complete set, we have 


(117) [a,; ay |= Onars [a,, api la, a], 0, 


the signs — and + corresponding to bosons and fermions respectively. We 
have also: 


(118) Ne Na Nos ce 
3 


The N, have the eigenvalues 0, 1, 2... (co) in the case of bosons, but in the 
case of fermions the only eigenvalues are 0 and 1. By introducing the ex- 
pansions (116) into (114) we get 


(19) H=— D apay dd |puayy+ 
BUSA 


IL 
= > ’ De Ai jn Azıv dj € 13 2 P(x, ee x) la" Ie Se 
LIA! MER, A 
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with 


(120) a) p2 ay = fup(tou(r)p?àr 


(121) 24, 2° | Dax) 2", A) = ETAT (Te) B(x, — + )AT dr. 


2 


In the representation in which the N, are diagonal, the wave functional of the 
field is a function of the eigenvalues N, of the N,. The functions f(r; @,) 
introduced in section 4 are orthogonal and normalized but do not constitute 
a complete set. Nevertheless, if the cells w, are small in all the dimensions, 
the set of the u,(r) = f(t; w,) will be approximately complete and we may 
use the formalism corresponding to the expansions into the functions of a 
complete orthonormal set to get a description of incompletely specified states 
of the field. In this treatment we have: 


; iL ff 
(122) CK | p? | A N Pi On zu pi = — | P° dr ; 


cat, P(x, — x2) | A", AT = Dur a Org On" anv 


(123) x te 3 
| D;ran == —— | ar, | D(x, a Xp) Ot, , 
OO" 
op! az 
1 (w) (wm) (w) il 
(124) Ho = 5 > Ny pi += > at N zu By 1zr = > NO . 


In the above formulas we did not assume that the volumes of all the cells 
are equal. 


(@), (w), 


We shall now assume that each incompletely specified state 7,(N,; Na», ---) 
has an entropy S[x,| 


(@) (w) (@) 


(125) Sl x0] = KDE Lea)? log) SUN), (PE x N) |? = 1). 
À N, A £ 


N, 


Let us consider the y determined by the conditions of extremal entropy with 
given average number of partieles N and given average value of the energy BP: 


(@), 
(126) = (5%) ’ Ds Pr > > Ya (N | 2N; 
(w) 


N; 


a N) (@) pi il (w) (@) It 
(127) 13 Ve = a D | Xo(N je | My (Ee IM Bu) + Na Mar Para |. 
ae 
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By introducing the Lagrange multipliers ka, kß and ky we get the equation 


(128) ù (SL Lo] ane kal Xo Nee) tri KB Loos Bog es) . ; ky (Xoo ah =0 ’ 


hence 
(@) 
19) u) = 
(wy Ih LA (w) Pi il (@) (@) Il 
== exp —= (y + aE Ny —p{ SN, (Be — 5 Da a) +3 > Ny NO ae 
A! al 2m 4 ary 


The formula (129) is obtained directly by giving a variation to | 7, |? in (128), 
without splitting it into the factors y* and y,. By variating separately x”, 
and x, we get: 


(130a) Gol LOG | Kol 2 el ON ee Ce 
(1306) ER Leo Ll en Bee 


These equations are satisfied as a consequence of (129) by taking y, «, P and y 
real. We shall denote the 7 thus obtained by 7: 


vCan * 


(131) Xeon = exp 
@) (pi, 1 fl (0) (0) 
Fe Ni (PE is 2 734 Ny ND ul] = 
It follows from the first equation (125) that: 


can | 


(oy 
(132) Sean =k 21 PR y+ + Le: > N; SF 


(@ ; Dirar 1 (@), (@), 
= aba ak en - Ny Na Du. | & 
2 a 7! 


FL 2m 7 


By taking into account the second equation (125), (126) and (127) we get: 


(133) Sn) = ky’ + aN + BE) y=7y +1. 


Therefore we have 


(134) DE RL Be 


: 
=k 
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T, u and P being the temperature, chemical potential and pressure, respectively. 
The incompletely specified state described by yx, corresponds thereby to a 
canonical grand ensemble. 

In the particular case of non interacting particles 


(@) ey 2 
- ae 
(135) Io (WY? = € gh exp [— Ny(« + Bex)] ; (ex = =| . 
so that the average number of particles in the cell @, is: 


4 (@) N 


can can 


(136) N - 5h 


(©) (w) (wm) 
— ex > N, exp - N;(« + pes) LE > exp |— Nyr(x + Per)| - 
A N 
Ar A 


7: aA 
Ny 


In the case of fermions the N’ take only the values 0 and 1, hence: 


(137) No == et exp a0 Ber] IT [1 + exp[— «— Bex]] = 
ME 


ar BR 1 


ae a | (wo) N' ee Er er 5 
exp [— «— Be] $ D Kean (N)! exp [x + Pa] +1 


© (@) (0) À 
> N;exp te - (2% - Bes 


(138) LN, SOX — N3= =. TER > ik > | ,,(&) le = 
ve NAVAL CSD | 77 © (@) 7 Xean + ie ae 
> exp |— N;(x + Bo ; 
N;=0 
0 Ce] 
= Lol $ exp [— Mie + Ben]! = 
0m a = 
a 
0 Tee } 1 
= — o 4 —- SS ee  , 
OB ts ORB LO Peal een fac Bey] 


8. — The Boltzmann method in the classical wave mechanics. 


In sections 4 and 5 the Boltzmann division of the phase space into finite 
cells was combined with the von Neumann method of the density operators. 
We shall now apply to the wave mechanics in phase space the method of 
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expansion into series of orthogonal functions which do not constitute a com- 
plete set. The wave function @* will be replaced by an approximate one 
O™ represented by a series 


n,@ 


1 
159, OS ee D YEA ape) Ae ee ee ee eee 


nw N ! OEE PROD Te An, 


the D, being defined by (78) and (79). The normalization condition for 
O® is: 


ñn,o 


(140) RC RCE ne Je Gdn tA, PEER 
D RE ere Er 


o n 

Ar and A are respectively symmetrical and anti-symmetrical with 
respect to the indices À, À,…, An. We shall now consider the following de- 
finition of the entropy: 


1 does n 


k we: 
(141) N ss SW Ges an ae 
Ale 


The average value of the energy in the state 0° is: 


Nn,0 


(142) HI Oe.) = [ | On (Ti: wees Tn) | Bu, SOC, a) dr, eine dr, — 


a 


< (+) © 
(43) HE, = [ (Data, toy Tan, pracy Oy MPH la) Cry der 


a 


The condition of à system in thermal equilibrium can be defined by the 
oe" whose coefficients are: 


(144) A CD É Soe | (Ai AE. <A, for fermions). 


It is easily seen that the entropy S)*)°"of the state O7" is: 


[0] 


(145) SE =; HO) Na. 
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Hence, by taking as usual 0 = kT, ne will be the Helmholtz free energy. 
The Y= are determined by the normalization condition (140): 


Nn,0 


(146) pee =: Pos; 2, Ge (Areas An )exp 


These expressions of the free energies are equivalent to those obtained in 


(n) 
oO 


section 4 for the density operators R 
The above method can also be applied to systems of distinguishable particles 
by using the following expansion for Or 


(147) Goltz. Tn) = D, Ad. D Tay Un; Day 5e OA) 9 
Are A 
(147a) > 9+(Aı, RG) A5) ln ua, ls 
AIR 
1 
(148) DIS cs Tns DY, 3. (072) in Alan 3,) .. En; Op): 
Vo" 


The average value of the energy in the state © 


n,0 


(149) BIC, ol = cae . 2. Lan (A. Am) Ag, SOS pu He 


grees 


(1494) Hy, ...,4, = [De Tiny WA,» «++ on) Hat Ta) dr 


Q 
We shall take as entropy: 


(150) Se = > 0 (ER) An) Ady ccs 18 log ! in) fone 


EN 


é ye o H 
(151) AU, — exp| Bact 
ne 20 
The free energy WY, is: 
(152) Ph OR 0 log | > Gs (Anse An) CXS al 5 
AudzparsAn 


30 - Il Nuovo Cimento. 
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9. — Statistieal treatment of the phase space fields by density operators. 


The definition (51) of the canonical density operator must be modified in 
the application to the fields, because of the existence of the integral of mo- 
tion N,,. The canonical density operator is now characterized by rendering 
the entropy extremal with the conditions of given average values of the energy 
and the number of particles, whereas the expression (51) is obtained by ren- 
dering the entropy extremal with the only condition of given average energy, 
as well known (1°). By applying the method of the Lagrange multipliers and 
taking into account the von Neumann definition of the entropy (48), the 
expression (49) of <H) and the following equation 


(153) EINE 


we get the variational equation 


(154) d(S— kBKN > — kBxH)) = — k I ög.[1+ log da + 6:N, + BE] = 0, 


kp, and kß, being the Lagrange multipliers. Thus we get 


(155) qa = exp[—f,N,— Ba — 1], 


and by taking into account the normalization } q, = 1 we may write 


a 


exp [— BN, — pHa] 


156 Ve an = ex —BN'— BE, 

( ) q > exp = BLN u BE x ply Pa Pa ] ’ 
with 

(157) CES exp | AN BE: 


Finally we get, instead of (51) 


(158) Roan = EXP [7 — BN op — Bs COM 


(5) Kôpr Husımı: Proc. Phys. Math. Soc. Japan, 22, 264 (1940). 
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since the energy operator of the field is denoted by #. The entropy of R,,, is 


(159) Sean =— &  (9— AN — BB.) exp [9 — AN — BE] = 
=e k(n m BLN ENT ra Bx FO» can)» 
hence: 
u ii TAY 
eae! Lie apne calla = EP ene) NETTER 


The canonical ensemble for the field corresponds to a canonical grand ensemble 
for the particles, as shown by (156). Thus the statistics of the field leads 
automatically to the introduction of the grand ensemble. This fact shows 
that the grand ensembles are of a more fundamental nature than the petits 
ensembles in the theory of indistinguishable particles. The quantal formula 
(158) can be taken over without modification into the theory of the « quant- 
ized » fields. 
It follows from the formulas (60) and (138) of I that 


do) exp [— BN op] — > exp [— B,n]P; = 
n=0 
ls > exp [— fin] (Pa, tn co) LC area 
n=1 . 
a 


P, being the projection operator on the linear manifold of the eigenfunctionals 
of N,, corresponding to the eigenvalue n and P(t,, ...,T,) the projection ope- 
rator on the direction of the vector y(t,,...,7,) in x-space: 


af 
(162) Ey en) 


== ee Cay) er Ta) à 


À 
(163) Pit ++) Tn) = 75 y*ltı) + P*(Tn)Pop(Tn) «++ Y(t) - 


The operators Ay and 3,4 being commutable, we have: 


(164) exp EE Pa] = exp = BeHyin] exp = B:P pot | © 


Let us consider now a numerical function F(r,,....,7,) and the corresponding 
operator Fon: 


(165) Fop = a JF . Tr) (Ta) ce PAT) PTT) «++ P(t) AT, ... Ate. 


a 
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Is was shown in I that 


| Ile CN corte 


(166) N Inne, 
: | (2 ae e M Yee) for nÆn, 

and 

(167) Pets, (") | POS POS de de 


By taking into account (166) and (167) we get 


(168) tante > "| lt, PE) den. 


so that 


nt (eo) 2 ln} 
(169) exp[—» Fp] = SP: + > | exp |— v(") Fu, ee Pa ee 
‘=r y 


n=r, 


y being a number. We get from (169) 


(170) exp [— Po Hin) = Po + > | exp Be Pa SP; A, cdd 
n=1, L AM y = 1 
(171) exp ES Pa Hot | es P; = Jen SF 


TAGE Or ooo Ok. 


+ S | exp 
2 


Q 


AT eet 


so that: 
a 
(172).  exp[— BA] = P, + | exp, pi] P(r) du + 


A 


FAST | exp [— BH, (Ti Ta) Pl, 2 Tr) AT,» den. 


n=2. 
2 


By combining (161) and (172) we get finally the spectral decomposition 
of Roan? 


(173) Roan = Poe? +> | exp [4 —- Bit — BoH alts» «+; Tn) |P (tay 2 Te) dr... dm. 


2 
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Equation (173) shows that we must take: 


(174) CHENE ss exp [— fn] [exp [— Bob (T1 + Ta) | AT... Ta - 
n=1 , 


Q 


Thus we get the same 7 for both fermions and bosons, as should be expected 
from the results obtained in section 4. A difference in the thermodynamical 
behaviour can be obtained by the introduction of finite cells. This can be 
done by a suitable modification of the density operator Ron: 


10. — The super-classieal field. 


We shall now prove that the theory of the « quantized » fields in the phase 
space © can be obtained by a kind of « quantization » of e-number fields in Q, 
Let us consider the variational principle 


ty 
(175) 2 = 6{ Ade, 


to 


x 
(176) A=i ] yet) 3 polls 7) dr — | we (t; 7) Lan (t) pelts 7) AT — 


ö a 
J ; 3 / / ! I 
4 | we (ts Thpé(t; T )Lpot(T, T)yelt; T w(t; t) AT dr 
ip © sa OD (a x) 0 .0D(x—x') © 
al 4ki ary a Rae] OGG) Jar 5 % d ; 
( Gh) Lxin(?) ah oe Ly UT T ) v dx Op rv ax’ op’ ’ 


for the c-numbers y(t, t) and ye (t, 7). The Euler equations are: 


a p 2 
(178a) at w(t; T) = — = ag Pi; a) 42 
@@(x—x') 0  0P(x—x) 0 | 
al *(t: 7! Br SEE Mel - ! 
1 futur] Ox aye Ox! api a SE ? 
0 p © 
= *(t: ey RL: te 
(1785) An We (t; T) m dx Wo (t; T) 7 


(ed — x), 0 Pa — =), 2 a En PA RE 
fil az (Op |e TE RARES 
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We can develop a hamiltonian formalism in the ordinary way. The eonju- 
gated momentum of y.(t; r) is JZ, (t;r) and the hamiltonian X.: 


© 


7) 
(J 79) IT, (t, T) = ae = ips (657), 


: 9 ( 0 wie | 


= w(t; 7) AT — Alt) = 


= | we (t; tT) Lxin(t)pe(ts T) AT + 
2 


oun 
ate | pe (ts tT) pe (t3 T')Loot(t, yet; t' p(t; t) dr dr” . 


2 
The equations (178) can be obtained as Hamilton equations: 


GATE OR, Ô k ER, 
me) RATES ee 


The field described by y. will be called the super-classical field. Our « quan- 
tized » field in phase space can be obtained from the super-classical field by 
replacing the y. and 71, by operators y and 71, satisfying the commutation 
rules: 


[ Cyté; x), IL, (6; v')]4 = iölt— 1’) 
[yd; 7), pl; t) ds = LT; T), IT, (t; Eyl — 


(182) 


Thus we get the commutation rules of the «quantized » field in the phase 
space 2 (Heisenberg representation). The Hamilton equations 


. oy OA Olly Ki 
(183) à Ot = [y, A] Mir. aap: =[I,, R) = —i— 


give for y and y* equations of the same form as (178). From now on we shall 
use the same notation for the corresponding quantities of the super-classical 
and quantized fields, except when it will be necessary to distinguish both cases. 

From the fact that our fields are in a six dimensional space or a seven 
dimensional phase-space-time follow some interesting results. Thus, instead 
of a four vector of density and current, we have a seven vector. Indeed we, 
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get from equations (178) or their analogue for the « quantized » field 


A 


g 5 P : 
34 i L p* Æ 
(184) ap (YY) + div, (v nm”) 
— div, (ys | | p*(t; T') oe p(t; t') dr’ v) — 0, 


so that the time component of the current seven-vector is y*y, the space 
components are y*(p/m)y and the momentum components — p*{] pr; T')° 
2 


-(0D/0x) p(t; t’) dtr \y. By integrating both sides of (184) over the momentum 
space, we get the ordinary continuity equation 


+ © + 


(185) 4 | wap) = aiv. je? pap) — 0, 


m 
- © 


whereas the integration over the configuration space gives the continuity 
equation in momentum space: 


(186) “(| ptpax) + ai (— ay ae eC eye 
A pry x) 7 OLY, | —— dxy* | Y Neha À de py) = > 


y L 2 


By integrating over the entire phase space we get: 


(187) ler Oz 


2 


This equation in the case of the « quantized» field shows that N,, is an 
integral of motion, and in the case of the super-classical field that the total 
mass is constant. 

The « quantized » field in phase space gives a classically exact description 
of the motion of systems of indistinguishable particles. The description cor- 
responding to the super-classical field is less accurate than the classical one, 
but it may be useful in the discussion of the behaviour of matter in bulk when 
the fluctuation effects arising from the molecular structure become negligible. 
In the super-classical approximation there are no particles but only a con- 
tinuous distribution of matter with density mo,(t, x): 


4 co 


(188) 0.1, x) = | | pelt; v)|*dp - 


-@ 
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The microscopie kinetic and potential energies of the super-classical field are: 
(189) FO Saud | JET) DE 
(19% kin an, | We 5 7); P P, 

2 


1a? | 
(190) = à] o.(t, x)o-(t, x’) B(x — x') dx dx — 


y 


bol ke 


=3/ | p(t; T) |?! p(t; t') |? D(x — x’) dr dr’ 


2 


It is easily seen that, as a consequence of (184), we have: 


(191) a I (zen, == FES +) = 0. 


The potential energy 4,,, can not be rigorously expressed as the space integral 
of a density of potential energy, because of the existence of actions at a di- 
stance. It is nevertheless possible to get a coarse-grained potential energy, 
which is the space integral of a coarse-grained density of potential energy, 
when the forces have short range, by dividing the space into cells V_ large 
enough to take: 


(192) Ma 2 A 
a 
(193) KU (Vo) = = | Olt, æo.lt, x')O(x — x") dx dx’ . 
Fe 


It results from (184) that 


3 
(194) „tet; ta} + div Je ot; Tps | En 
© 


with 


(195) ol; t) = p*(t; Ty; 7) , 
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for both the super-classical and « quantized » fields. Hence we get: 


+o + © 


R 5 
(196) Al o(t; Typ, dp + aiv. |? o(t; t)pxdp + 


| o(t; t') D(x — x’) dr’ y(t; 7) dp = 0. 


In the case of the super-classical field there is no difficulty in introducing a 
macroscopic velocity, of the matter. W.(t, x): 


+o 


1 
(197) al 0.(t; T)pdp = a(t; x)W.(t; x). 


—-o 


This is not possible for the « quantized» field, because the operator 


+0 


pre; t)y(t; t)dp does not admit an inverse. With W. we get in the usual 


way the total derivative d/dt: 


d a 6 
1 => : 
Asst dt ct 0x 


It is easily seen by, taking into account (185) and (197), that 


+o 
A 


(199) Al oct; JE dp = p.(t; x) LA 


ot 


er — W. div {o.(t; x)W.}, 


-o 


and since 


; N ee (Be Pen, do LV ot 
ee [Een Bay = | (F w.) Vea) outs D dp eWex0.lt; %) > 


— © — © 


we get from (196) the equation of motion 


AW Oleg. a Be 
CO a — ells #) ap, | ot: x!) D(x — x") dx’, 
r 
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the kinetie stress tensor being: 


+ © 


2 a ((®- ee soit: t) dp - 
(202) sp” — m | (Pe — Wa.) (2 — Weg) outs 2) dp 


— © 


The equation of motion (201) corresponds to the absence of external forces. 
There is of course no difficulty in introducing such forces. 

The consideration of the macroscopic equations af motion leads to a macro- 
scopic definition of the kinetic energy different from the microscopic kinetic 
energy given by (189): 


(203) HE i o.(t; x)W? dx . 


The explanation of this difference is immediate by taking into account that — 
a part of the total microscopic kinetic energy corresponds to heat. We shall 
define the local temperature T,(t,x) by the condition that the total micro- 
scopic kinetic energy of a macroscopic element of matter in its rest-system 
of reference be equal to */,k of the number of its particles multiplied by 
the local temperature: 


(204) = | 0.(t; T) (2 Æ A dp= : ko.(t; x)T(t, x). 


It is easily seen that: 


+O +0 
m 


ä p = ¥ Cy m 
5 fete T) (P = w) dp = | et T) on: dp — 5 o(t, x)W2. 


oo — © 


(205) 


Hence the density of microscopic kinetic energy is the sum of the density of 
macroscopic kinetic energy plus */, ko,(t; x)T.(t, x). 
We shall now take as density of entropy s,(t, x): 


+ © 
(206) (LA = =k | ot; 7) log o.(t; t) dp . 


=. 


This entropy has the usual unsatisfactory behaviour of the Gibbs entropy: 
the total entropy of the super-classical field is time independent. Indeed, let 
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us compute the time derivative of s.(, x) by means of (184): 


a ie a 
(207) En SE, x) = — x | (los dir) + Ua e.(t; T) dp = 


+0 
= k div, {o-(t; x)W.} + k [ 108 0-(t; T) div. | oc(t; DE | dp — 
- x | tos 0.(t; T) aiv, | 0.(t; | 0.(6;%.) dr) dp — 
u ia) 
EAN : 5 
— EAN, | o.(t; T) log o.(t; tT) dp . 


[ee] 


Tt follows from (207) that for a volume V’ of boundary 2’ we have 


+ 0 


of { 
(208) | SL, X) dx — [az [(n) o.(t: T) log o.(t, t) dp , 


y’ x 


n denoting the outer normal at the points of the boundary 2”. In the part- 
icular case of V'— V the integral in the right hand side vanishes and the 
total entropy remains constant. It is easily seen that equation (208) can be 
written as follows 


+0 


Leal es ere sae [az [r ı? je 


7 


w.\) o.(t; T) log o.(t; tT) dp — 


“ a 
y’ yy 


— | (n-W.)s,(t, x) dd’. 


yy 


Both terms in the right hand side of (209) are fluxes through the boundary 

and this shows again that there is no creation of entropy, with the definition 

(206). Since the second term corresponds to the flow of matter, the first must 

correspond to the flow of heat. It is important to notice that the first term 

in the right hand side of (209) behaves satisfactorily under a change of units, 

because the contribution of the additive constant due to log o, is nil. 
Equation (184) for the super-classical field can be written as: 


RP Ne Ry 
(210) (5 FAP A =) alt: T) — 0, 
à 5 
(211) F (t, x) -— | Q(t; t') P(x — x')dt'= — = p(x) 
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Equation (210) is the same derived from the Vlassov formalism (?) for a system 
with a large number of similar particles. In the case of a stationary motion 
of the super-classical field we may apply Vlassov’s method of solution by taking: 


2 
(212) ot) = 0s(X)0.(E) , (e = pe) : 


It is easily seen that 
(213) o(é) = ae”? , 


x and p being constants related by the normalization condition: 


(214) [ o(e)dp — 1. 
0.(x) is of the form 
(215) O-(X) = exo) exp [— Bip(x) — p(xo)}], 


the potential g(x) satisfying a non linear integral equation: 
(216) p(x) = 0x0) | D(x — x') exp [— B{y(x) — p(xo)}] dx’. 
5 


We shall not discuss here the solutions of the Vlassov integral equation (216). 
The parameter ß is related to the temperature, according to Vlassov, by the 
equation: 


1 
217 = ——. 
se P= EF 
We shall now prove that this follows from the general definition (204). In the 
case of a stationary solution (212) the macroscopic velocity vanishes 


(218) W(t, x)= 0, 


and (204) becomes simply: 


+ 0 + © 
(219) : Ee) fotos dp =— « a exp [— Pe)dp = = : 


— © 


It is easily seen that the stationary solutions of the type (212) are cha- 
racterized by the condition of rendering the entropy extremal with a given 
energy and a given mass. The entropy corresponding to the definition (206) 
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is therefore convenient for the discussion of the equilibrium, notwithstanding 
its unsatisfactory behaviour for non stationary motions. 

In the case of short range forces it is convenient to express the second term 
in the right hand side of (201) in terms of a potential stress tensor defined by 
the equation: 


0 NE ” ap Sr: ! 
(220) EB — 9.(t; x) | 0x; &)- est PES 
OL A OL x 


Equation (201) becomes: 


Won 


6 
291 vof, es eal 
Wen leer aarp 


(kin) (pot) 
(Dra. Fre 


13. — The Boltzmann equation. 
Equation (184) can be written as follows: 


(222) fe ee elle: T)— pis 9.(t: ot. 


| at Fm Ox J | o.(t; T')D(x — x’) dr’. 


Ox 

Le) 
Equation (222) has a remarkable resemblance with the fundamental equation 
of the kinetic gas theory, the Boltzmann integro-differential equation. The 
difference consists in the replacement of the right hand side of (222) by the 
Boltzmann collision integral. This change is of course of capital importance, 
as shown by the fact that (222) leads to the time independence of the entropy, 
whereas the Boltzmann equation leads to the H-theorem. 

We shall now obtain the Boltzmann equation in à Way which allows to 
see clearly the different assumptions on which it depends. It is often assumed 
that the entropy increase allowed by the Boltzmann equation results from 
the incomplete specification of the motion or else that it is due to the hypo- 
thesis of molecular chaos. Neither the incomplete specification nor the mole- 
cular chaos lead to an increase of the entropy, as we shall see. Let fr(f; T15 ---, Tw) 
be the distribution function in the phase space of a system of N identical 
particles. fy satisfies the Liouville equation and is normalized to N!: 

fx a ee en 
(223) + Hay 0, | f(t} Tay oe) Tw) dti ee Oty = NE. 


Q 


The r-particle distribution functions 


1 \ 
(224) Pl; u, Tr) = RE | fr(bs ty Tr) Appa Ai 5 
2 
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satisfy the equation of Yvon (8), BORN and GREEN (?) and KIRKWOOD (1): 


(225) à JONES Tan secte) A its ut ee u 


= At 441 | > D(X p44 a= x), fy P(t; Tips u) 
U r+t 


2 = 
The hypothesis of molecular chaos consists in assuming that: 


N—1 N—-1 
(226) (U5 Try Te) = nr CG GP; ne) = y JG IG; 2) 


In this case we get from (225) (*): 


Bory (2 4 Pe. © \ 9g) a wana 
(227) Es \ra N Is jit; 1) | sg, | Hts Pl Sag: 
2 


For large N this equation is equivalent to (222). The factor (N — 1)/N arises 
from the exclusion of the self-forces, as is clearly seen for N = 1. 
The Boltzmann definition of the entropy of the N particles system 


(228) SAT fic 1) log f(t; t)dr, 
5 


associated to the molecular chaos hypothesis leads to a time independent 
entropy, when the variation of the one particle distribution function f is 
computed with (227). This results from the reasoning of section 12, since (227) 
goes over into (222) by replacing the potential ® by ®N/(N—1). Therefore 
the time variation of the entropy is not a consequence of the incomplete 
specification of the motion associated with the description by a probability 
distribution, and the situation is not changed in this connection by the in- 
troduction of the molecular chaos hypothesis. The essential step is the re- 
placement: of the third term in the left hand side of (227) by a suitable time 
average which transforms (227) into the Boltzmann equation. We shall now 
define more precisely this time averaging. 


(*) It is necessary to assume that the spatial distribution is not homogeneous, 
otherwise the third term in the left hand side of (227) will be nil. This indicates 
that (226) is not quite satisfactory an assumption. In section 14 a better method 
will be given. 
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It follows from the definitions (177) that equation (227) can be trans- 
formed into the following 


A Neeson |? ; 
(229) at f(t; %) = — IN: | D, {f(t; %) ft; Te) } dr, + (U N— 2) Lxin(Tı) f(t5 7), 
a 
the operator D, being: 
(230) DE, te) = Rb a) Dun) Loe (ts T2)} E(t, Te) « 


Let us consider a function F(t,(t), t.(t)) of the two points tft) and T,(t) 
in Q, which correspond to a two-particle system with the interaction 
@(x,— x), the two particles having the same mass m. It is clear that 
D,F(1(t), t2(t)) is simply the time derivative of F(x (t), t2(t)) taken along the 
motion of the particles. We shall assume that the forces have a range small 
enough to be treated as a first order infinitesimal quantity in connection with 
the space variation of F(t, 7). The time of collision corresponding to a 
collision finishing when the two particles are at the points 7, and 7, of 2 will 
be denoted by 6(t,, 72) and treated as a first order infinitesimal quantity in 
connection with the space variation of F(t,, t). In other words we shall take: 


(231) F(t, (¢— OG) tte 6)) = F(x,(t), Pilt— 0); x2(t), p:(t—9)) = 


Of 0 | = ler, 7 = & 
~ m |? u(t) ° 2x) + Pel = als (x, (8), Pı(t 9); x2(t), pat 0)) ; 


(232) 0 = Art), t2(t)) . 


It follows from (231) that, for a slow spatial variation of F, we can take: 


(233) | dr; (t) IE = ‚low D,F(n(t'), t2(t')) ar \= i Dygn(ta(t)) [Foo 1,(t)) draft) — 


Q 2 


-1; { F(x,(t), pilt— 0); xt), palt— 0)) — Frl), t2(t))} deat) = 


Q 


1 Ce ä 
PAG)" en) Fr), T(t) dt,(t) — 
2 


m 
— L [apa je x (8), Pat — 9); xt), Pe(t— 6)) = Frl), T2(t))} | pa(t) — pit) | dqe(t) ’ 


(234) dg.(t) = dad’ (t) da? (t) , xs axis // p(t) — pt) - 
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We shall now replace in the right hand u of (229) Difflt; t1)f(t; t2)} by its 


time average (1/6) | D,{f(t: m(t')) f(t; 72(t'))} dt’, making use of the equation 
1-0 

(233). Thus we get the Boltzmann equation for the modified distribution 

function f(t; 7,): 


R G Pi Q 
23 + —|j(t — 
ee Is m oO | (tm) 
we = = 
u. | dp |: ft; x1, pt) f(t; x2, pt) — fit; u) ft; ve) — p2| dq, 
2 F 
(236) Pi = pilt— 0"), po = p,(t— 0"), = O(ty, Ta) - 


The above derivation of the Boltzmann equation is less elementary but 
simpler than the ordinary one, since it does not require the discussion of the 
direct and inverse collisions. This advantage becomes clearer when collisions 
between particles with internal degrees of freedom are considered. Another 
important feature consists in the fact that triple and higher order collisions 
do not come in in our method of time averaging. This corresponds to the usual 
approximation in the derivation of the Boltzmann equation, in which only 
binary collisions are taken into account. 

The Boltzmann equation leads to the H theorem, which shows that the 
entropy may increase but never decrease. Thereby it is clear that the increase 
of the entropy defined by (228) results from the time averaging. It has already 
been shown by KIRKWOooD (!!) that the derivation of the Boltzmann equation 
requires a time averaging, so that the Boltzmann equation is satisfied by a 
time-coarse-grained distribution function. The method of Kirkwood differs 
essentially from ours, because we did not average the distribution function 
over a time interval along the real motion, but averaged instead the effect 
of the interaction along a virtual motion of a two particle system. Our method 
of time averaging is simpler and in a certain sense quite naturally suggested 
by the form of the equation of evolution of the distribution function. It must 
nevertheless be recognized that the need of a time averaging does still require 
a deeper justification. The current point of view is based on the consideration 
of the time necessary to make an observation (see KIRKWOOD (!°)), but it is 
hardly conceivable that the increase of entropy should depend essentially on 
the fact that our observations require a time large with respect to the duration 
of the molecular collisions. It does not seems unreasonable to think that the 
definition of the entropy is not altogether satisfactory for non equilibrium 
states and requires some modification. This is suggested by the « quantized » 
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field formalism. Indeed, from the point of view of the «quantized » field 
formalism the analogue of the Boltzmann definition would be the introduction 
of some kind of entropy operator with a structure corresponding to that of 
the entropy of the super-classical field: 


(237) S,=—k| 0.(t; 7) log o.(t;r)dr. 
0 


2 


It does not seem possible to define an operator log o, but it is possible to 
divide © into cells and to take as entropy operator 8: 


op’ 


op 


(238) Ss = —k > N, log N,. 
À 


The entropy of a state of the field would then be the average value of the 
entropy operator (coarse-grained entropy): 


(239) ete re ({ Ixzledu=1). 
lt -space u -space 


It is easily seen that 


BADEN 


ODL 


=) > >| ERA) 
2 


IN (0, Tyan Un) LO N (0; Ty, u, Tn) DTT, 


the ©, being the coefficients of the Fock expansion (112) and N (0,5 Ti, Tn) 
the number of points 1,,..., T, inside the cell w,: The definition (239) is 
obviously not equivalent to those discussed previously. Formula (125) sug- 
gests that, if an entropy is to be attached to a pure state x, it should be 
taken as S[ x]: 

(241) Sly] =—k | |x|*log! x|*du- 


space 
We shall discuss this point further in a forthcoming paper. 


14. — It is interesting to notice that the method of time averaging of 
section 13 can be applied directly to the exact equation 


0 0 © 0 
(242) (= = i x=) f(t; T) Sr fn (t5 Tay Te)” 2x, Dix, — X2) AT: , 
2 


Le 


31 - Il Nuovo Cimento. 


470 M. SCHÖNBERG 


without introducing the molecular chaos hypothesis. Thus we get a gene- 


ralized Boltzmann equation involving modified distribution functions fand 7%: 


(243) r HB ie) = 


et m Cx, 


+0 
1 = x 
Sa Jap. | a (t; Xi, pt; x2, Pi) — =; T1 Tz) }| Pa — p: dqp - 


— © V 


It is preferable to perform first the time averaging and then to intrudoce 
in (243) the molecular chaos hypothesis in terms of the modified distribution 
functions f. 

Let us consider now the equation (225) for r= 2. It is easily seen that 


this equation can be written as follows: 


2 i 
(244) Fr fy (t; Ti, Te) = | (Hat, T2, Ts); rl; T15 Te, Ts))s dt; — 


SON = 3)(B,(, Ts), fe; Ti, T)) o= Te | D? fats Tas To Ta) AT, — 
a 


za 3)(H;(tı, Tz), alt; Ti T2))2 - 


The operator D is the three particle analogue of D;: 
(245) DFE T2, Ts) = (F(t, Ta, Ta), Hat, Te, Ts))3 : 


It seems plausible that the same kind of time averaging which led to the 

Boltzmann equation will also give satisfactory results when applied to (244). 

We shall replace in (244) DP f(t, tm, Ts) by its time average during the col- 
t 


lision of one of the particles with the couple of the other two (1/6) | DIE 
t-9 
(nt), T2(t'), t9(t')) dt’, assuming that the collision finishes at the time ¢ with 
the particles at the phase space points 7,, 7,, T; and denoting by 0 = 0(t%, T2, Ts) 
the time of collision. Since f@(z,, tT, T3) is symmetrical with respect to the 
variables of the three particles, we may assume that particle 3 collides with 
the couple (1, 2), the distance between 3 and the other two being larger than 
|x,—x,|. In order to simplify the results we shall neglect the difference 
between the forces exerted by particle 3 on the other two, so that the variation 
of p, = Pi— Pz will be taken as due only to the forces between particles 1 
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and 2. We shall also assume that the variation of p, is small. With the above 
assumptions the positions of all the three particles and p, may be taken as 
varying little during the collision and it is easily seen that: 


(246)  Flnlt—0), 129), rft—0)) = 


~ 


(x x,(t— 0), x,(t— 6), x3(t— 0); p-(t— 9), Pr(t— 9), P3(t— 0)) = 


(Ps = Pr + Pr) . 


This equation replaces now (231), F is assumed to be a slowly varying function. 
The analogue of (233) is 


t 


(247) [an IE i DPF (rt), rt), T(t’) dt’ 


a t- 0 


m — 2 = F(r;(t), T2(t), Ta(t)) dr.) — 


= dp;(t) X,(t), x(t), x(t); Pr), Prt— 9), pat — 9)) — 


— © 


— F(n(t), t2(t); Soy | Ps(t) — po(t) |Aqa(t) , 
(248) dg;(t) = da(t) da@(t), a axis // pa(t) — p,(t) . 


It follows from (247) that the analogue of (234) for the two particle distribution 


function is: 


aes à 
(249) a (t; T1, 2) + ( ill T1 T2), Halt, T2))2 = 


-,| dp | (t5 X15 Pi} X25 Pis Xss Ps )— 196; Ty T3 T)}|Pa— Pı— P2| qs 5 
= y 


al ' 
(250) pt= 5 {mt 9) + pb}, Pi = {pl 9) — pift)} P3= palt—0"), 


0'= OT, Te, Ts) - 
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We can also get the analogue of the Boltzmann equation by using the 
three particle analogue of the hypothesis of molecular chaos, the superposition 
principle of KirKwoop and BoGGs (15): 


Pts toy Ta; Tey MOUS Ty Ta) 


ft; ta)f (ts Te) f(t; Ts) ; 


(252) WE Tr, To Ts) = 


I thank Prof. PRIGOGINE for many useful discussions which helped me to 
clarify several points of the present paper. 


(6) J. G. Kırkwoop and E. M. Bocas: Journ. Chem. Phys., 10, 394 (1942). 


RIASSUNTO (*) 


Il lavoro sviluppa ulteriormente la meccanica classica delle particelle indistingui- 
bili trattate in I. Si discutono i fondamenti meecanieci della termodinamica in termini 
della teoria ondulatoria classica e dei campi « quantizzati» nello spazio delle fasi. Si 
dimostra che questo formalismo elimina automaticamente il paradosso di Gibbs. Si 
tratta la statistica classica col metodo di Gibbs e von Neumann nonché con una forma 
generalizzata del metodo di Boltzmann che consente di tener conto delle interazioni. 
L’introduzione di celle finite nello spazio delle fasi conduce a formule simili a quelle 
della statistica quantistica per particelle libere ma non cosi per particelle interagenti. 
Il trattarnento statistico del campo « quantizzato » nello spazio delle fasi conduce im- 
mediatamente al grand ensemble canonico. Si dimostra che la teoria dei canıpi « quan- 
tizzati » nello spazio delle fasi si puö derivare dalla teoria di un campo «non quantiz- 
zato » nello spazio delle fasi, in eui si hanno solo distribuzioni di materia continue, per 
mezzo di un processo di « quantizzazioue » del tutto simile a quello della teoria quanti- 
stica dei campi, ma tuttavia escludente la costante di Planck, o qualsiasi altra costante 
universale. Il campo «non quantizzato » nello spazio delle fasi corrisponde a una teoria 
meno esatta di quella della meccanica classica e da una approssimazione dello stesso 
tipo di quella data dall’ipotesi del caos molecolare della teoria cinetica dei gas e com- 
prende sia la dinamica sia la teoria del calore di un mezzo continuo. Si dà una nuova 
derivazione dell’equazione di Boltzmann che richiede un tipo speciale di media tem- 
porale e si stabilisce una funzione simile per la funzione di distribuzione di due par- 
ticelle. 


(*) Traduzione a cura della Redazione. 
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NUOVO CIMENTO Vor. X, N. 4 1° Aprile 1953 


Sull’energia di legame dell’He’ e del Li‘. 


G. MORPURGO 


Istituto di Fisica dell’ Universita - Roma 
Istituto Nazionale di Fisica Nucleare - Sezione di Roma (*) 


(ricevuto il 3 Febbraio 1953) 


Riassunto. — Facendo uso di funzioni d’onda descriventi un deutone ed 
una particella « ((3) e (4)) l’energia di legame del Li® viene calcolata al 
variare della distanza tra il baricentro del deutone e quello della parti- 
cella «. La massima energia di legame, quando si impongano le usuali 
condizioni di saturazione, si ottiene per una distanza infinita (figg. 1 e 3) 


cid essende dovuto al fatto che l’energia potenziale cresce — in valore 
assoluto — troppo poco al diminuire della distanza tra il deutone e la 


particella «. Di questa lenta crescita sono presumibilmente responsabili 
le eondizioni di saturazione. Analoghi calcoli vengono svolti per l’Heÿ. 
In accordo col momento magnetico nel Lif lo stato S (3) dä luogo per ogni 
distanza ad una maggiore energia di legame che lo stato P (4); nell’ He® 
capita l’opposto. 


ntroduzione. 


Il presente lavoro contiene una discussione del problema dell’energia di 


legame dell’ He’ e del Lif. Esso si riallaccia ai lavori di INGLIS (1), MARGENAU (?), 
TYRRELL (?), HUMBLET (+). 


E noto che se in un’espressione della forza tra due nucleoni della 


) 
) H. Mareenav and K. CARROLL: Phys. Rev., 54, 705 (1938). 
3) W. TyRRELL Jr.: Phys. Rev., 56, 250 (1939). 

) J. Humprer: Physica, 14, 285 (1948). 
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forma: 
(1) Jr.) =— Af () {w-+bP,+hP,P,+mP,}, 
(2) wtbt+th+m=1 


si aggiustano i parametri in modo da essere in accordo (°): 

a) con Venergia di tripletto e di singoletto del deutone (cid implica due 
relazioni indipendenti tra A, ae w + m); 

b) con l’energia di legame dell’H? (terza relazione tra A,a e w + m) 
allora Venergia di legame della particella « risulta essere alquanto troppo 
grande (°). 

Ma non appena dalla particella « si passa alla considerazione dei nuclei 
ad essa immediatamente superiori le energie di legame calcolate con la sopra- 
detta espressione della forza risultano estremamente inferiori al vero. Cid & 
in particolare per il Li‘, il primo nucleo stabile superiore alla particella «. 
La ragione di cid è dovuta in parte al fatto che i calcoli diventando più com- 
plicati, essi devono essere svolti in maniera meno perfezionata; ma in parte 
è dovuta anche al fatto che il Li® é il primo nucleo nel quale a causa del prin- 
cipio di Pauli, le restrizioni imposte sulla forza dalle condizioni di saturazione 
dei nuclei pesanti si fanno sentire (nel caso dei nuclei fino alla particella « 
compare nell’espressione dell’energia soltanto la combinazione: w + m). 

Lo studio del Li® si presta percid allo studio delle proprieta di saturazione 
della forza. 


1. — Diseussione dei lavori precedenti ed impostazione del problema. 


Per l’Heÿ sono stati fatti due tipi di calcoli dell’energia di legame (*) (4): 


a) calcoli nei quali si parte da un’approssimazione zero nella quale i 
cinque nucleoni vengono, tutti e cinque, riguardati come muoventisi nello 
stesso campo di forze elastiche ed in cui la differenza tra la forza fittizia ela- 
stica e la somma delle interazioni tra i nucleoni viene trattata poi perturba- 
tivamente ; 


(5) I parametri cosi determinati risultano essere per una f(r,,/a) gaussiana, come 
assumeremo nel seguito: 
A = (35,6 MeV; a = 2,25 1078 0m; 2 a 2 m =10;8.. 


(5) Vedi, ad esempio, L. ROSENFELD: Nuclear Forces (Amsterdam, 1948), p. 305, 
per la discussione di questo fatto rilevato la prima volta da RARITA e PRESENT (Phys. 
Rev., 51, 788 (1937)). 
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b) calcoli variazionali nei quali la funzione variazionale à caratterizzata 
da due parametri l’uno associato, per cosi dire, alle dimensioni della parti- 
cella x, Valtro, che indicheremo con 1/f, caratterizzante la distanza del neutrone 
esterno dalla particella «. 

Questo secondo tipo di calcolo è il più interessante e conduce alla conclu- 
sione che la massima energia di legame si ha per / = 0, ossia per una distanza 
infinita tra il neutrone e la particella « (7). 

Nel caso del Li® sono stati fatti soltanto i calcoli del tipo a); ossia i sei 
nucleoni vengono riguardati in approssimazione zero come muoventisi indi- 
pendentemente in uno stesso campo di forze elastiche, e la differenza tra la 
forza fittizia elastica e la somma delle interazioni viene trattata poi perturba- 
tivamente. Il risultato di questi calcoli & che Venergia di legame ottenuta è 
sempre assai inferiore a quella effettiva (31,94 MeV). Nel caso del Li® tuttavia 
non sono stati fatti calcoli in cui si possa esaminare l’andamento dell’energia 
di legame al variare della distanza tra i nucleoni esterni e la particella «. 

Scopo di questo lavoro è di colmare questa lacuna e presentare questi cal- 
coli. Pit preeisamente, come risultera dalla forma delle funzioni variazionali 
usate, che scriveremo nel num. seguente, il Li® verrà riguardato come com- 
posto di un deutone e di una particella x. E l’energia di legame sara espressa 
in funzione di un parametro rappresentante le dimensioni della particella «, 
un parametro rappresentante le dimensioni del deutone ed un parametro rap- 
presentante la distanza del baricentro del deutone da quello della particella «. 

Si noti che una descrizione del Lis del tipo ora detto & consistente con il 
fatto che il Li® possiede un livello eceitato a 2,12 MeV (*) e con il fatto che il 
momento magnetico del Lif, é, grosso modo, eguale a quello del deutone (°). 


2. — Funzioni variazionali usate e cenni sui calcoli. 


Le funzioni variazionali con le quali per il Li® e per I’Hes si & fatto il cal- 
colo sono le seguenti: 


— Li’: 
3) BD, = 3-1 Pull )a(1)w(2)B(2)o(3 jae(3)0(4)B (4) u( 5) (5 J0(6 )ox(6) [56), 


(4) D, = > (— 1)’Pu(1)a(1)u(2)(2)0(3)a(3)v(4)h(4)u(5 )x(5) (6) (6) |56), 


(7) Nei calcoli in questione TyRRELL ha fatto uso di un’interazione gaussiana ed 
HungBLer di un’interazione yukawiana; 1 risultati sono qualitativamente gli stessi. 

(8) W. M. Harris: Phys. Rev., 84, 1249 (1951); F. AJZENBERG e T. LAURITSEN: 
Rev. Mod. Phys., 24, 321 (1952). 

(2) II momento magnetico del deutone è 0,86 magnetoni nucleari, quello del Lif 
0,82 (vedi, ad esempio, J. BLATT € V. Waisskopr: Theoretical Nuclear Physics (New- 
York, 1952)). 
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dove u, v ed «x, ß sono rispettivamente le funzioni di spin isotopico e di spin, 


> indica la solita somma sulle permutazioni, e si & posto: 
P 


(5) 156), = exp — or + ris + ia + re + Von + 184) — 
Tg ou. Us SR ER Sr LA - 
a les ae 
€ % J 
07 
(6) 156), = Pr, tr) (nr Pre r,)]z|56): - 


Tanto (3) quanto. (4) sono invarianti rispetto a traslazioni (1°) e corrispon- 
dono ad uno spin pari ad 1; (3) è uno stato S, (4) uno stato P. Il significato 
dei parametri o, ß e 6 nelle (5) (6) è, dopo la discussione precedente, ovvio. 


— He: 
(8) Fy = 3% (— 1)?Pu(1)ac(1)u(2)6(2)r(3)a(3)o(4)6(4 (5 )x(5)] 5). , 
(8) [i gee > (— 1)?Pu(1)a(J )u(2)3(2)v(3)a(3)o(4)B(4)u(5)o-(5)[5), 
dove 
(9) 15), = EXP (re + Tis + Pia + Vas ra) + 
a LPs Sie Uy ie ae 
ED fr 1 ) ' 
(10) 15), = Ar, (nn +r, + rs + ra):15). 


Senza entrare nei dettagli dei calcoli diamo i risultati di essi che possono essere 
utili in future ricerche su questo problema. 

Con la ovvia notazione: (ik|@|hl), indicheremo l’elemento di matrice di 
un generico operatore @ tra funzioni di tipo (3); analogo significato avranno 
i simboli (ik|@|hl),, @]@]|r)., (@]6|R),, in riferimento rispettivamente alle 
funzioni (4), (7), (8). Nel seguito affiggeremo di un indice s quantita riferen- 


(1°) Il calcolo perturbativo dianzi accennato degli Autori (1) (2) (3) & stato fatto 
in coordinate individuali sottraendo a conti fatti l’energia cinetica del baricentro; seb- 
bene questo procedimento sia abbastanza ragionevole non & tuttavia chiaro che genere 
di errori si compiano in esso. 
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tisi a calcoli eseguiti con la funzione S (tanto per il Li® quanto per ’He’), di 
un indice p quantita calcolate con la funzione P. L’operatore G puo essere in 
particolare l’unità, nel qual caso interverranno simboli come: (ik|hl)., (ik|hl),, 
(i|r)., (i{h),; 0 pud essere l’energia potenziale J(r,,) tra due nucleoni n, m; 
nel qual caso interverranno simboli come (ik|I(rmn)|hl), che abbrevieremo in 
(ik|mn|hl),, ece. L’energia coulombiana non sara nel seguito considerata. 


2:1. — Caleoli relativi al Li‘. 


Il valor medio dell’energia potenziale, V,, è una combinazione lineare degli 
elementi di matrice che compaiono nella prima colonna della seguente tab. I, 
il coefficiente che moltiplica ciascun elemento di matrice essendo la opportuna 
combinazione di w, m, b, h, determinata dalle altre colonne della tabella. 
(Esempio: (16|12|26), va moltiplicato per: %w + 8m + 4h + 4b). Una 
volta costruita la ora detta combinazione lineare, V, si ottiene moltiplicandola 
per — A e dividendola per il fattore di normalizzazione della funzione (3) Ns; 
questo & pari a: N, = (56/56), — 2(16|56), + (13]56).- 

La tabella IT contiene le espressioni esplicite degli elementi di matrice che 
compaiono nella prima colonna della tabella Me 

Il valor medio della energia cinetica T, & dato da: 


ide: 
uU) T= =i D'T,P, dr = 
L 


6 
= = S 46617156), —-20617:156), + (13]7:156) 
N ; 


Ss 
e, poiche à: 


He ED 11 ee 
; AT 


1 


A A) A 
O = 0 270 
—, + 160° — +80? || 
ae Co 00 


=: 2 2641284720 +38 = 36); , 


il calcolo di T, & ricondotto alla valutazione dei nove elementi di matrice che 
si trovano nella tabella III. 

Nel caso della funzione P il fattore di normalizzazione &: N, = (56|56), — 
— 2(16|56), + (13156), ed il valor medio dell’energia potenziale si ottiene in 
modo completamente analogo al caso della funzione S sostituendo N, con N, 
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e la tabella II con la tabella IV. L’energia cinetica media T, & data da: 


= 6 
(13) T, — =; | ¥ ,(56| 7,156), — 2(16| 7,156), + (13 
p 1 


T,|56), 


e la sua valutazione & facilitata dall’osservazione che é: 


A a 


+80" 5] 150), 


res fis? ks) ; > ) ER 

(14) > 2 :156 p—= ams B + 120 +- 126 + SP? ae 4- 160° + 

ed è percid ricondotta agli elementi di matrice che si trovano nella tabella V. 

Dalle formule precedentemente date si pud verificare, ad esempio, che al 

tendere di B a zero tanto E,= T, + V, quanto E, — T, + V, tendono a 
EH, + HE, con: 


S N? . 4 o se 
(15) FE, == M O0 6: (u => m) I wil 5 

> h? Ô We 
16 i) = — 4 Mate i == La —— _— R 
(16) Ey ai 30— A(w + m +b +h) F a = 


Queste espressioni coincidono con quelle date da altri Autori per l’energia 
di legame della particella x e del deutone. 


22. — Caleoli relativi all’He°. 


Il valor medio dell’energia potenziale U, è una combinazione lineare degli 
elementi di matrice che si trovano nella prima colonna della tabella VI, il 
coefficiente di ciascun elemento essendo la opportuna combinazione di w, m, b, h 
determinata dalle altre colonne della tabella. (Esempio: il coefficiente di 
(111414). è: w + 4m + 2h— 2b). Una volta costruita la ora detta combi- 
nazione lineare, U, si ottiene moltiplicandola per — A e dividendola per il 


fattore di normalizzazione n,; questo & pari a: (5|5),— (2]5).- 
La tabella VII contiene le espressioni degli elementi di matrice in questione. 
Il valor medio della energia cinetica: 


are 11 
C=) D] MER .Ae = ED, (6|P13).— @I » 
s 1 


RE 


si ottiene facilmente tenendo presente che à: 


he 


15 
k = ou 


ae re 
a 2001100 PS er: 15): ; 


4 


(17) Sy, Ti|5 
vi 


onde è ricondotto agli elementi di matrice che compaiono nella tabella VIII. 
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Nel caso della funzione P (8) tutto procede in modo analogo; il fattore 
di normalizzazione è: n, = (5|5),— (2]5),; la tabella VII va sostituita con 
la tabella IX; e la valutazione del valor medio della energia cinetica: 


Ty = I 57:15), — (217:15),} 
1 


Vp 
è facilitata dall’osservazione che é: 


2 


‚2 25 
T5), = 720 + 160? — + L Hp? — | 15}, 


S Na: 
(18) > le 


SH o 


a 


ed & percid ricondotto agli elementi di matrice che si trovano nella tabella X. 
Si pud controllare che tanto €,= U, + ©, quanto €, = U, + ©, tendono, 
al tendere di 6 a zero ad E, (15) come deve essere. 


3. — Valutazioni numeriche e discussione. 


Le valutazioni numeriche sono state fatte nel caso dell’He® e del Li® per 
un valore di o tale da rendere massima separatamente la energia di legame 
della particella «. In unità in cui il range delle forze & assunto eguale ad 1 
esso &: o = 0,33; ed il corrispondente valore calcolato dell’energia della parti- 
cella « è 25,21 MeV. Per 6 si è scelto pure il valore che massimizza l’energia 
di legame (16) del deutone; 6 risulta (casualmente) (1°) pure eguale a 0,33 e 
la corrispondente energia di legame del deutone è: 0,9 MeV. Questo basso 
valore dell’energia di legame del deutone dipende dal notorio (1%) fatto che 
una gaussiana è una cattiva funzione d’onda variazionale per il deutone. Per f 
si sono considerati, tanto nel caso dell’He® quanto in quello del Li® i valori: 
0, 1/24, 1/16, 1/8, 1/5. 

Per i parametri di saturazione si sono scelti i valori: 


(1) Sebbene tale stato P sia in contraddizione con il momento magnetico (vedi 
nota (?)) abbiamo svolto il caleolo anche per tale stato in modo da poter esaminare 
in dettaglio il meccanismo che da luogo a maggior legame nello stato S che nello 
stato P, al contrario di quel che avviene nell’ He’. 

(2) In effetti i valori di o e 6 che veramente massimizzano le rispettive energie 
di legame differiscono di qualche percento; per semplificare i calcoli numerici abbiamo 
attribuito a o e 6 lo stesso valore. 

(38) H. MARGENAU and K. CARROL: Phys. Rev., 54, 705 (1938). 
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in accordo con: 


(20) b+h=0,25 


(21) 4w + 2b—h—m —= 0. 


Questa condizione è la più favorevole, cioè quella che da più energia di 
legame, compatibile con le condizioni di saturazione dei nuclei pesanti. B facile 
constatare, con l’aiuto delle tabelle I e VI, ed avendo alla mano i valori nume- 
rici degli integrali, che una volta imposte le condizioni (2), (20), (21), le varia- 

zioni dei parametri 


S i compatibili con esse 
© Li z ; ‘ 
À 61) ENERGIA CINETICA € POTENZIALE non influiscono prati- 
© 

Selen funzione p camente sulla energia 
. 607 —— funzione s : 5 

x di legame del Li® e 
iat dell’ He’, 

c 

4 I valori dei para- 


metri (19) coincidono 
con quelli degli altri 
Autori. 

I risultati sono ri- 
portati nelle fige. 1, 
2, 3 seguenti, tanto per 


la funzione S, quanto 
per la funzione P. La 


fig. 1 da, al variare 
di Bi valori dell’ener- 


Fig. 1. 


gia cinetica e dell’ener- 
gia potenziale per il Lif, la fig. 2 gli stessi valori per l'He. E la fig. 3 da, 
tanto per il Li® quanto per l’He’ la variazione con ß dell’energia di legame. 

Si vede dalla fig. 3 anzitutto che la massima energia di legame si ottiene per 
ogni 5 con la funzione P nel caso dell’He® e con la funzione S nel caso del Lif. 
Cid è in accordo con il modello a shell. 

Ma cio che più si nota & che Venergia di legame del Lis decresce monoto- 
nicamente al decrescere della distanza. La massima energia di legame si ottiene 
quando il deutone è a distanza infinita dalla particella «. 

Che nell’He® la massima energia di legame si ottenga quando il neutrone 
è a distanza infinita dalla particella x & ragionevole in quanto l’Heÿ é instabile. 
Ma € spiacevole che nel Li® capiti la stessa cosa. La ragione sta, come si vede, 
nel fatto che la energia potenziale cresce (in valore assoluto) molto meno della 
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energia cinetica. Cid & dovuto essenzialmente alle condizioni di saturazione 
le quali danno luogo a troppo forti interazioni repulsive. 
Ci si pud ‘chiedere come si modificherebbero le curve delle figg. 1, 2, 3, nel 


S 

x He® 
S534 ENERGIA CINETICA E POTENZIALE 
© 

‚Oo . 

a ere funzione p 
= A — funzione s 
2 

i=.) 

L 

= 

Ww 


757.08 


(naW) 9/DIZU9}00 DI6J8U7 


caso che modificassimo i parametri w, m, b, h, rinunciando alla condizione 
di saturazione (21). E possibile, rinunciando a tale condizione, in particolare 


assumendo egual dose di 
forze di scambio e non di 


A 


Energie dilegame dell’He5e Li®con 


2 le funzionisep (in MeV) ie 
scambio (tipo Serber) spie- x a Calais 
gare simultaneamente la 2” - — Er — = : 
stabilita del Lis e l’insta- Sa] (He = ZHEDERER 
bilità dell’ Hes? La rispo- 5 | He® [252 |24.2| 238] 227 Me 
sta sembra possa essere 2 AS 
affermativa in virtù del a = 
diverso momento angolare 
orbitale nei due casi, che Li? 
fa si che l’energia poten- funzione p 
funzione s 

ziale cresca nell’He? assai 
piü lentamente che nel M 
Lie Lae a oe es Wap 

Speriamo comunque 
di poter presentare in fu- Fig. 3. 
turo più estesi calcoli nu- 
merici a partire dalle formule qui stabilite, si da poter rispondere alla 


domanda sopra posta. 
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Per il momento la conclusione, valida sotto Vipotesi naturalmente che le 
funzioni d’onda variazionali qui usate siano sufficientemente aderenti alla 
realta, &, che le condizioni di saturazione quali si scrivono di solito, sono 
incompatibili con Venergia di legame del Li‘; questo fatto, unito alla loro 
incompatibilita con le esperienze di scattering n-p a 90 MeV e un altro argo- 
mento a favore della idea che la saturazione dei nuclei pesanti debba essere. 
spiegata con qualche altro meccanismo. 


Desidero ringraziare il prof. B. FERRETTI per varie utili discussioni avute. 
Inoltre sono grato all’Istituto Nazionale per le Applicazioni del Calcolo, nelle 
persone del Direttore, prof. M. PICONE, e del dott. W. Gross, per l’esecuzione 
dei calcoli numerici a partire dalle espressioni algebriche contenute nelle tabelle. 


TABELLE RELATIVE AL Li® 


TABELLA I. 


Coefficienti di w, m, b, h, che moltiplicano à vari integrali tanto nel caso della funzione P' 
quanto in quello della funzione S. 


mn m | b | h | 
(12]12|12) +1 +1 dei | de ÿl | 
(2 'cpope .| a ag ee 
Lune] ee eal | © Peaster 
x (6112115) SF ; — 6 ETES ty ic — 2 ia 1.6 
Tran. | an eR ce 
(151226) OS AR alae el 
2112156) ee ene ee 
er ee ee ae u | ae hen 
 as1256) ES ee | 
(3441256) | ae re. | re | 
APR >, fu OUI ERP a ÉTAIT TSS N RONA 
(36] 12]56) 2% 26 | 0 | 0 | 


(56112156) + 6 + 6 0 | 0 
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TABELLA II. 


Integrali che compaiono in V, e, per «= 0 in N,. Tutte le espressioni qui date sono da 
moltiplicarsi per wn inessenziale fattore numerico. Il simbolo El (— 3/2)[] indica che 


LAS r Apres enero! parentesi va cry alee alla aes, — Su 


(12]12]56), = 


(34112156), 


(12/12/12), 


(56112156), 


(15112126), 


(16/12/15), 


(12112126), 


(6112116), 


(13112156), 


(36112156), 


(16112126), 


(16112156), 


I 
= 


as Be, (40 + 20)804 (x -| do + 20)(40 +78 + 2 aaa! 
2 4 er APN J 


3 : 4 3 A . 
— }|(se ! B - 3 ) sr| 24028 + 96020 + 18080 + 24oxù + 


Il 
= 


er] 
DETTE 
RE 

= u (perte à 29 

| 

| 


3 


Il 


pi (— ler le; poet =) 
3 3 9 2 i re 
(ses Fe); B+ 30 + se] | + (sr abe 


3 
= El (- 3) | 640? {64080 + Soßx + 320040 + 6uBd} | 


| 


— El (- a je + 26) 1 [320? + 6o« + 1600 4- vol ae + 24) — pl + 


: DEN 
EI - à [see + 2a) (ao +1 en : | (are + 1208 + = ») 

3 E SAME 3 Ge 
= pa (— 3) He (sr +5 = B+38+ (a6 + sa) (40450) 75 \ 


= El (- sl 2 “8 + : also = = 2]|(se + i; 6+ ») = = ae 
9 
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TABELLA III. 


Integrali che compaiono in T,; stessa avvertenza che per la tabella IT. 


| 


(56102(/00)|56), = — 128-72 7120408 


(56 | 8?(8/2B)|56), = — 128-24 r"/??00® 


HI 


(56 | 62(0/26) | 56), = — 128-24 r/28620° 


9 
(13 | o2(@/@o) | 56), = — 241-5/202(20 + 0) ce HAE a [60 + 6] + 20(20 + ö)- 


(160 + sp | 


(131B2(0/28)156), = — 481 12ß2o(20 + 6)? (50 a i 


9 
(13 | 6°(2/86)|56), = —- 481%0ö%(20 + 6) bt +1086 Te ps 
(16102(2/20)156), = — 23120? = + 320° | 9600 + 160° + 2208 + 24Bo + 7 By 
2 o 
3 | u. 9 
SEI an H = BO + 36° + 120° + 5 Bo 


3 27 
(16]62(2/86)|56), = — : ator aot | 2200 iat ais Te UO). pr} 


r = (64)?096B t= 80(40 + 26)? (+ + = ie pl 


2 = 64024 de 45 < ot oe : z 
2 = 040 ve ph +59) (Fa + 29)(40 + Lo) rl 
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TABELLA IV. 


Integrali che compaiono in V, e, per « = 0, in N,. Stessa avvertenza che per la tabella IT. 


(56112156), 


(1211212), 


(16/12/16), 


(34112156), 


(12/12/56), 


(13112156), 


(36112156), 


(15112126), 


(12|12]26), 


(16112156), 


(6112115), = 


(16112126), 


bol vw 


-32-256-020(80 + 2«)(56]12|56)>/2 
5 04: 256-0%(40 + 2x) (12] 12] 12)8/2 


-16:6402(6400 + 8ox + 6x6) (16|12|16)5/8 


38 4 
— 5:5: 1440(do + 20)%(80 + 22)(34]12]56)27 
3 8 1 
— 55'144 Bo%(4a 4 20 + 2a)(4o-+ 20) (12]12]56)2° 
3 ae ota a 32 
— 144: (40 4+ 26) 5 © {6x0 L 320? + da + 1600} + ao? + 1e 208 : 
(13 ] 12 |56)8/8 
ö 4 
— © .144-80(80 + 26) Ê (40 + 36) + 2 (36]12|56)5/3 
: tet 4 oO 4 8 : a 
©.144:80(40 + 26) |2(40 + 20) kr oF )- te (15| 12] 26)° 
 64-14408|(40 + 36) E a, 0, #00 (2112126358 
= À AES Te ee, = ae I 5 
MA 9 A VIRE : 
32 1 


-(16112156)5/8 


k 8 
-144 80 - (320? + 2406 + 4a0 + 20x) (5 o + 


8 
) + 


à 


+ 3(2a + sx (se ae | (16 | 12] 15)5/2 


; iy 4 N 
144 : 6402 É ar à (do + 36) + 5 a (16]12]26)3/ 


32 - Il Nuovo Cimento, 


486 


G. MORPURGO 


TABELLA V. 


Integrali che oempaiono in T,. Stessa avvertenza che per la tabella II. 


(56 | 0?(2/00)|56), — 


(56 | ö2(8/26) 156), 


(5618°(0/08)156), 


(13]02(0/80)]|56), — 


(13] 6°(0/20)| 56), 


(13]6°(0/0B) [56), 


(16]0?(0/dc)|56), = 


(16] 6°(0/06) | 56), 


3 
(16 62(2/0B) |56), = > 


Le espressioni di 


| 


3 004! 


wl © 


(56 | 0?(8/00) |56), 771: 00 + 3r-5/2g45 


ee 


| © 


(56192(2/00)156), 71003 + 8527-512 


w 


u 56] Bx(0/2B)156),r-1 | 
zy a al | 6°(6/0B)156).r i 


no] 


- 1024 | = 5 (13 | o2(0/d0) | 56),t-1(20 + 6)20? — o%t-5/2- 


-[20?(20 + 6) + o(20 + om | 


3 5 
Bo — 5 (13] 9%(2/00) 156) (20 + 0)2o2 — 04-5262 + 6) | 


I or 


; (13182(0/28)156) 420 + 8)°0° | 


3 5 
5250 2 (16|04(2/20) 56),2-107(400 + 362 + 1262) + 


+ 02275/2(6002 + 3620 + 24%) | 


K 


5 5 
64-16-f (16) 94(0/20)| 56) too =| 302 == 120?) =e 


+ 026%27312(36:+ 20) i 


64-16 | 3 (16/B2(2/08)156),2-102(400 + 302 + 1208) | 


» r, t sono date nella tabella III. 
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TABELLE RELATIVE ALL’ He’ 


TABELLA VI. 


Coeficienti di w, m, b, h, che moltiplicano i vari integrali tanto nel caso della funzione 8, 


quanto in quello della funzione P. 


| 
a m b | h 
Le Sei lees Beets 
(511415) 6 6 0 0 
(4}14]4) 4 | | 2, = 
| (211415) | —3 —3 | 0 0 
(411415) | — 6 = 0 | 0 | 
(1|14]4) =, va = ro | 


TABELLA VII. 


Integrali che compaiono in V, e (per x = 0) in N, . Stessa avvertenza che per la tabella I. 


(511415), = El — 5) [ows = B(80 + 2 


els, = mi 3) [sore + 20 far + 138) — 54 | 

3 & Te os Il 
anal). = #1 (— 5 { (0 5) ee 3°) 5 Bt) + Bond + à ol) 
(111414), = El (-3) to { (40 + a + = A) (= F set 


N.B: Non abbiamo riportato l’espressione di (4]14|4), perehè in base alla tabella VI 
e alla condizione di saturazione (21) va moltiplicata per zero. 
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TABELLA VIII. 


Integrali che compaiono in ©,; stessa avvertenza che per la tabella II. 


(5]02(9/°0)|5), — — 480?BT5/2(80)2 

(51B2(0/0B)15), = — 22T —5/2(80) 

(4|0?(2/20)|5), = 2.9408. 9-52 (40 + = ß (se + 2 Ip 
é 2 12 12 9 


(41 82(2/08)|5), — — > 28-3280): = seep 
owe 2 3 144 


B(80)* S = 640? (sr Je uv 8) — id pl 
12 9 


TABELLA IX. 


Integrali che compaiono in V,; stessa avvertenza che in tabella IT. 


(511415), = 64(80)?(80 + 2x) (5[1415)5/8 


(211415), — — 16(80°)(80 + 2a) (2]14]5)3/3 
(411415), — — 128-802(80 + 2«) (4]14]5)5/3 
(111414), = — 16-6402(80 — 3x) (1114/4)5/ 


Vedi N.B. alla tabella VII. 


TABELLA X. 


Integrali che compaiono in ©,; stessa avvertenza che per la tabella Il. 


(5]0°(8/20)|5), — a HOLSODDTD EEE + 12878 | 


(51PKEJOMI3), — 64(80)8 > T1 (5]842/28)19), 


(4102(2/00)|5), = — 16-6402 | : (4102(2/80)15),808-1 + 12 | 
5 
(41 °(8/°B)]5), == — 18507 S~? (4]8°(2/0)|5). 


Te § sono dati in fondo alla tabella VIII. 
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SUMMARY 


By making use of wave functions describing a deuteron and an « particle ((3) 
and (4)) we have calculated the dependence of the binding energy of Li® on the distance 
between the center of mass of the « particle and the center of mass of the deuteron. 
The largest binding energy is obtained (if the nuclear parameters are subjected to the 
saturation conditions) when the deuteron and the « particle are an infinite distance 
apart (see fig. 1 and 3); this is due to the fact that the absolute value of the potential 
energy increases much slower than the kinetic energy when the distance decreases; 
probably the saturation conditions are responsible for this slow increase. For purpose 
of comparison analogous calculations have been done for He? (fig. 2 and 3). In agree- 
ment with the magnetic moment, in Li® the S wave function (3) gives rise, for each 
distance, to an higher binding energy than the P state (4); the opposite is true for He’. 
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Electrons Associated with Negative : Mesons (*). 


NG, INA 


Department of Physies, University of Wisconsin, Madison, Wis. 


(ricevuto il 4 Febbraio 1953) 


Summary. — A study has been made in photographie emulsions of low 
energy electrons associated with stopped negative u-mesons from the Uni- 
versity of Chicago cyclotron. A total, of 1000 meson endings has been 
studied. Of these mesons, 358 decayed into a high energy electron and 
32 caused a star. The remaining 610 mesons stopped in the emulsion 
with no associated tracks other than tracks of low energy electrons 
(10 < H,< 100 keV). A low energy electron track is observed from 17 
of the ~y-e decays (5 percent). The energy distribution of the low energy 
electrons from the p-e decays is consistent with the assumption that the 
low energy electrons are monoenergetic. The average energy of these 
electrons is about 17 keV. The low energy electrons are thought to be 
due to « Auger» transitions of the negative u-meson from the atomic 
capture in the light elements of the emulsion. One or more low energy 
electron tracks are observed from 255 stopped u-mesons which did not 
decay cr cause a star. These electrons are thought to be due to mesonic 
«Auger» transitions from the heavy elements in the emulsion. The 
large difference in the probabilities of electron ejection from the mesonic 
capture in the light and heavy elements (5 percent as compared to 44 per- 
cent) is confirmed by a study of low energy electrons associated with 
negative p-meson stars. 


1. — Introduction. 


It is generally accepted, that a negative 1 meson which stops in a solid is 
captured in various mesonie « Bohr » orbits of an atom and goes by various 
transitions to the. lowest possible mesonic level. The meson may then be 


(*) Supported in part by the Research Committee of the graduate School from 
funds supplied by the Wisconsin Alumni Research Foundation. 
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either captured by the nucleus or it may decay into an electron. The nega- 
tive meson many undergo « Auger » transitions as well as radiative transitions 
in the process of going to the lowest state. A number of investigators have 
studied, in photographic emulsions (1?) and in cloud chambers (°°), the low 
energy electrons which accompany stopped mesons. These low energy electrons 
are interpreted as due to non-radiative « Auger » transitions of the meson. 

In the previous studies the mesons were obtained from cosmic rays. Both 
positive and negative x mesons as well as u mesons are found in the cosmic 
rays. In the case of the emulsions studies, the number of negative u. mesons 
was estimated from the stopping power of the various constituents of the 
emulsions for very low energy mesons, the probability of nuclear capture and 
the number of u-e decays, if an assumption is made regarding the ratio of 
negative to positive y mesons. From the number of u. meson endings which 
have only an associated low energy electron track and the estimated number 
of negative y mesons, the probability of low energy electron ejection from the 
capture in the heavy elements in the emulsion could be estimated. 

From previous studies (*7) it has been found possible to obtain, from the 
University of Chicago cyclotron, a relatively high flux of negative u. mesons 
with a low negative x meson contamination. By determining the percentage 
of the negative y mesons which decay in the emulsion (°), the percentage of 
the. mesons which stop in the light and heavy elements can be computed. 
Because of the possible uncertainty in the number of negative {1 mesons in 
the cosmic ray studies, it seemed worthwhile to redetermine the percentage 
of stopped negative u mesons which give rise to one or more low energy electrons. 
Also it seemed of interest to determine if the atomic capture of a ~y. meson 
in the light elements in the emulsion leads to the ejection of electrons in the 
energy range which can be detected by photographic emulsions. 

Aside from the study of the fundamental processes that are invelved in 
the electron ejection, the low energy electrons which occasionally accompany 
meson capture have been used as a criterion for distinguishing negative meson 
capture in the light elements from the capture in the heavy elements (8). The 
validity of this criterion can be readily verified by a study of the association 
of low energy electrons with other phenomena caused by stopped u. mesons. 


(4) M. G. E. Cosyns, C. C. DiLworTH, G. OCCHIALINI and M. ScHONBERG: Proc. 
wns Soc., A 62, 801 (1949). 


@) W. F. Fry: Phys. Rev., 79, 893 (1950); 83, 594 (1951). 

(3) A. Bonerrr and G. Tomasını: Nuovo Cimento, 8, 693 (1951). 

(4) W. Y. CHane: Rev. Mod. Phys., 21, 166 (1949). 

(5) G. ee L. LEDER and F. Rise: Phys. Rev., 81, 626 (1951). 
(6) W. F. Fry: Phys. Rev. 85, 676 (1952). 
- (7) H, MoriNAGA and W. F. Fry: Nuovo Cimento, 10, 324 (1953). ; 
. (8) M. Menon, H. Muirnpap and O. ROCHAT: Phil. Mag.. 41. 583 (1950). 
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2. — Procedure: 


Ilford G-5 emulsions 400 microns in thickness were exposed behind ab- 
sorbers to the 122 MeV negative 7 meson beam of the University of Chicago 
cyclotron. The 7 meson beam contains also a small percentage of u meson 
from the decay of x mesons while still in the vicinity of the target. The 
u mesons were separated from the z mesons of the same momentum by utilizing 
the greater range of the u mesons in a Cu absorber (%). About 10* meson 
endings per em? of emulsion were obtained in the plates with an exposure 
of two hours. The background of tracks in the plates is quite low. It is 
estimated that more than 50 percent of the.tracks in the plates are due to 
negative mesons. 

The emulsion was searched with an over all magnification of 250x. Only 
those meson tracks which stopped well inside the emulsion and were longer 
than 200 microns were accepted for this study. Bach meson ending was 
studied carefully with a magnification of 980 X. 


3. — Results. 


A total of 1007 meson endings has been studied. In 7 cases it was not 
possible to determine whether the u meson caused the emission of one or 
more low energy nuclear particles or low energy electrons. The results of the 
study of the remaining 1000 meson endings are summarized in Table I. 


TABLE I. — Phenomena Associated with Stopped "u Mesons. | 
arr VUS | 7 
Phenomena Number of Events | 
| | 
| : _| 
| —p-e decay . eae een do 341 , atomic capture | 
| ~u-e decay with an additional one or 358 _ in the light ele- | 
| more low energy electrons - 7 ments | 
u ending with no associated tracks . 355 | 
u ending with one low energy electron 
| BLACK. SS ee ed Deter 180 | atomic capture 
~u ending with two low energy electron 610 | mainly in the | 
| VLAGRB ERREUR PCR 57 2 255 heavy elements | 
| ~u ending with three low energy electron | 
| TAG 8 RE rer ace ee 18 | 
=RIMEEON. SUAS i.e EE RE ees 32 | : 


A minimum ionizing track associated with the “u meson ending has been 
considered to be due to a decay electron, if the energy of the particle is greater 


acne 


Fig. 1. — A photograph of a negative u-e decay with a low energy electron is shown 
in the above figure. The low energy electron track is indicated by the arrow. The 
energy of this electron is about 15 keV. 
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than 5 MeV, as determined from the multiple scattering of the track. The 
expression «low energy electron » refers to an electron in the energy range 
from 10 keV to 100 keV. Only in 5 cases, an electron from a u meson was 
found to have an energy in the interval from 100 keV to 5 MeV. These 
5 cases have been included among the u meson endings with low energy electrons. 

The percentage of negative u. mesons which decay in the emulsion is found 
to be (341+17) (100)/(1000) = 36 +2 
percent which is in good agreement with 
the previous results (5) (39 + 3 percent). 
The emulsion consists principally of two 
groups of elements; namely a heavy group, 
Ag and Br, and a light group O, N, © 
and H. Mesons which stop in the emuls- 
ion, stop either in silver bromide crystals 
or in the gelatine so there is essentially 
no competition in the atomic capture of 


number/5 keV 


oOo 


the mesons by the light and heavy groups 
of elements. The probability of nuclear 
capture has been calculated to be propor- 
tional (?) to Z* for the light elements, with 
equal probability of capture and decay 
in elements (1°) with Z= 11, therefore 
essentially all of the negative p me- 


0 10 20 30 
electron energy in keV 


Fig. 2. — The energy distribution of 

low energy electrons, which accom- 

pany negative u-e decays, is shown 
in the above histogram, 


sons which stop in the silver bromide crystals are expected to be captured 
while 92 percent of the negative u mesons that stop in the gelatine are ex- 
pected to decay, if it is assumed that the atomic «catching » probability is 
proportional (1) to Z for the light elements of the gelatine. Thus it was 
found that 39 percent of the mesons which stop in a photographic emulsion, 
stop in the gelatine and 61 percent stop in the silver bromide crystals. 

In 16 cases the -u-e decay is accompanied by one low energy electron and 
in one case the -u-e decay is accompanied by two low energy electrons. An 
example of a -u-e decay accompanied by a low energy electron is shown in 
Fig. 1. The energy of the electrons was estimated from the measurement of 
the range and using the range-energy relationship (122). For electrons of 
energy between 10 and 20 keV, the estimation of the energy from the range 
of the track, is necessarily quite poor but it is thought to be superior to the 
estimation of the energy from the number of grains. The energy distribution 


(9) J. W. WHEELER: Phys. Rev., 71, 462A (1947). 
(0) H. K. Tıcno: Phys. Rev., 74, 1337 (1948). 

. (4) E. Fermi and E. TELLER: Phys. Rev., 72, 399 (1947). 
(2) B. Zasac and M. Ross: Nature, 164, 311 (1949). 
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of the low energy electrons associated with -u-e decays is shown in Fig. 2. 
The distribution is consistent with the assumption that the electrons are mono- 
energetic, however the results do not prove by any means that this assumption 
is correct. In each of the 17 cases the energy of the high energy electron 
was estimated from the multiple scattering of the track and was found to be 
greater than 10 MeV. Therefore it is extremely improbable that the high 
energy electrons were ejected during the atomic capture of the -u meson in 
the heavy elements. The high energy electrons in these 17 cases are un- 
doubtedly decay electrons from the negative u mesons. The probability that 
the low energy electrons associated with the ~u-e decays are due to the acci- 
dental association of 6 rays produced by the decay electron at the u-e junctions 
has been experimentally determined and found to be negligeably small. The 
low energy electrons are most probably due to non-radiative meson « Auger » 
transitions from the atomic capture of the negative u mesons in the light 
elements of the gelatine. The probability that the capture of a negative 
u meson in the light elements will lead to the ejection of one or more low 
energy electrons of energy greater than about 10 keV is then found to be 
(17/358) (100) = 4.7 percent. 

In 610 cases the meson stopped in the emulsion without associated tracks 
other than low energy electron tracks. One or more low energy electron tracks 
are associated with 255 of these meson endings. The energy distribution of 
these electrons was not obtained since the energy distribution has been pre- 
viously studied by several investigators (12>). These electrons are considered 
to be due to non-radiative transitions from various mesonie «Bohr » orbits 
in the heavy elements of the emulsion. The probability that the capture in 
the heavy elements will lead to the ejection of an electron of energy greater 
than about 10 keV is then (255/[610 — (0.08)(358)])(100) = 44 percent. This 
percentage may be low because it is not always possible to distinguish a very 
low energy electron track from the end of the u meson track. Nevertheless 
the probability of electron ejection is much higher from the capture in the 
heavy elements than from the capture in the light elements. 

The large difference in the probabilities of electron ejection from the capture 
in light and heavy elements is confirmed by a study of negative u meson 
stars. It has been found (?) that the bulk of the one prong y meson stars are 
due to the nuclear capture in Ag and Br while a portion of the multiply 
pronged u meson stars can be attributed to the capture in C, N, and O. A 
total of 75 one prong stars has been studied of which 39 (52 percent) have 
one or more low energy electrons associated with the star. Likewise 30 nega- 
tive u meson stars have been studied where the characteristics of the stars 
establish with reasonable assurance that the mesons were captured by the light 
elements. Only in one of these 30 stars is a low energy electron track observed 
to originate from the star. Although the atomic capture of a negative meson 
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in the light elements can lead to the ejection of a low energy electron, as evi- 
denced by the electron from -u-e decays, the probability of electron ejection 
is much smaller from the light elements than from the heavy elements; there- 
fore the absence of electrons from stars caused by stopped negative mesons 
can be used to distinguish capture in the light elements from capture in the 
heavy elements. 


The author is greatly indebted to Professor H. L. ANDERSON for his per- 
mission to use the facilities of the University of Chicago cyclotron. The author 
is also indebted to Mr. ARTHUR ROSENFELD for many interesting discussions 
and for his help in making the exposures. 


RAPATS SIUENSRON(CE) 


Si sono studiati in emulsioni fotografiche gli elettroni di bassa energia assoclati 
con mesoni pz negativi in fine corsa prodotti dal ciclotrone dell’Universitä di Chicago. 
È stato studiato un totale di 1000 mesoni in fine corsa. Di questi mesoni, 358 di- 
sintegrarono in un elettrone d’alta erergia e 32 produssero una stella. I restanti 
610 mesoni si arrestarono nell’emulsione mostrando solo tracce associate di elettroni 
di bassa energia (10 < E,< 100 keV). In 17 dei decadimenti ~p-e (5%) é visibile la 
traceia di un elettrone di bassa energia. La distribuzione energetica degli elettroni 
di bassa energia originati nelle disintegrazioni u-e è compatibile con l’ipotesi che gli 
elettroni di bassa energia siano monoenergetici. L’energia media di questi elettroni 
à circa 17 keV. Si suppone che gli elettroni di bassa energia derivino da transizioni 
« Auger» del mesone y negativo originate dalla eattura atomica negli elementi leggeri 
dell’emulsione. 255 dei mesoni y arrestati che non disintegrarono o non produssero stelle 
mostrano una o piü tracee dovute a elettroni di bassa energia. Si suppone che questi 
elettroni siano dovuti a transizioni « Auger» derivanti dalle catture atomiche negli ele- 
menti pesanti dell’emulsione. L’elevata differenza fra le probabilità di emissione di elet- 
troni derivanti dalla cattura di mesoni negli elementi leggeri e pesanti (5% in confronto 
del 44%) à confermata dallo studio degli elettroni di bassa energia associati con le stel- 
le dovute ai mesoni u negativi. 


(*) Traduzione a cura della Redazione. 
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(ricevuto il 21 Febbraio 1953) 


In a recent paper (1), CORTINT et al. 
have analyzed data from nuclear explos- 
ions observed in emulsions e posed at 
high altitudes. One of their conclusions, 
arrived at after a comparison cf their 
results with theories of meson product- 
ion, was that «the plural theory is not 
sufficient to explain the experimental 
data». This conelusion would be of 
great importance in the present contro- 
versy concerning the actual mode of 
meson production, but we shall show 
that it does not follow from their expe- 
rimental results. 

ÜORTINI et al. have based their argu- 
ments on rather a rough comparison 
(the best available at the time) of their 
results with the plural theory as formul- 
ated by HEITLER and JANossy. Sub- 
sequently, detailed numerical results de- 
rived from various models cf plural 
production have been presented in a 
form which makes possible a better 
comparison of theory with experi- 
ment (2). 


() G. CORTINI, A. MANFREDINI, and G. SE- 
GRE: Nuovo Cimento, 9, 659 (1952). 

(2) H. Messer, R. B. Ports, and C. B. A. 
MCCUSKER: Phil. Mag., 48, 889 (1952). 

(2) H. MESSEL, R. B. Ports: Phys. Rev., 88, 
488 (1952). 


It will be sufficient for present pur- 
poses to consider the integral multipli- 
city spectra of shower particles as con- 
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Fig. 1. - Integral multiplicity spectra of shower 
particles. s=number of shower particles, P(s) = 
= number of stars having > s shower particles. 
© and straight line are taken from CCRTINI et 
al.; + values predicted by plural theory, 


sidered by CoRrTINI et al. Their expe- 
rimental results are reproduced in Fig. 1, 
the straight line being the one drawn by 
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them as a best fit to the experimental 
points. Op the same figure are plotted 
the results predicted by us assuming 
plural theory on the HJ(1) model (in 
which it is assumed that the primary 
proton initiates a cascade within the 
nucleus and in each individual nucleon- 
nucleon collision one meson is produced 
which escapes from the nucleus without 
further interactions). The agreement of 
the two sets of results is quite good. 

We have obtained our results as 
follows. For energies H, greater than 
the primary cut-off energy B, (~ 3 GeV 
for latitude 47° N), the average number 
of protons with energ >E, is 


(1) N(B,) = kB;>, 


where k is a constant and y is the pri- 
mary power law exponent (we shall 
take y = 1.1, the value used by CORTINT 
et al.). From Table IV, reference 3, we 
can determine from the first and eighth 
columns the average number of shower 
particles s (or n, in our notation) pro- 
duced in a collision of a proton of energy 
E, with a silver or bromine nucleus in 
the emulsion. Hence we can plot the 
absolute number P(s) of stars having 
>s shower particles against the number 
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of shower particles s, normalising our 
results to agree with those of CORTINI 
et al. at s = 4.9. 

This gives 


(2) P(s) & 75(15/B,)"-! 


and the points plotted in Fig. 1 cor- 
respond to 
E, = 20, 30, 50, 


100, 200. 


In interpreting this comparison of our 
results with the experimental ones it 
has to be borne in mind that (a) the 
experimental results do not definitely 
distinguish collisions with light nuclei 
from those with heavy nuclei, whereas 
the theoretical results are for heavy 
nuclei (b) the statistics of the experi- 
mental data for large s are very poor. 
Because of these and other uncertainties 
we do not elaim that the agreement of 
the results is proof that the plural theory 
is correet; what we do claim is that the 
plural theory is sufficient to explain 
the experimental data, contrary to the 
conclusion of Corrinı et al. It should 
be noted that other models of plural 
production (25) could explain the results 
equally well. 
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Un esempio elettrodinamico di teorema di ortogonalitä alla Van Hove. 


A. LOINGER 


Istituto di Fisica dell’ Universita - Pavia 


(ricevuto il 6 Marzo 1953) 


E stato dimostrato da VAN Hove (1) che per un campo scalare neutro in inte- 
razione Scalare con un certo numero di sorgenti ferme ed infinitamente pesanti lo 
spazio hilbertiano degli stati stazionari del campo interagente con le sorgenti & orto- 
gonale allo spazio hilbertiano degli stati stazionari del campo libero. Questo teorema 
di ortogonalita (a proposito del quale v. anche (?)) rende evidentemente priva di 
senso, in generale, l’applicazione della teoria delle perturbazioni al problema con- 
siderato. 

In questa lettera ci proponiamo di far vedere come un analogo teorema di orto- 
gonalita valga anche nel problema di un elettrone in interazione col campo elettro- 
magnetico, ove ci si metta nelle approssimazioni non relativistica e di dipolo ed 
inoltre si trascuri il termine dell’hamiltoniano contenente il quadrato del poten- 
ziale vettore. Sono queste, com’é noto (?), delle condizioni che permettono, grazie 
alla trasformazione canonica di Bloch e Nordsieck, di risolvere esattamente il pro- 
blema detto. 

Rinormalizzata la massa ed eseguita la trasformazione di Bloch e Nordsieck, 
l’equazione di Scliôdinger del nostro problema si scrive, nella rappresentazione 
in cui X e le Q sono diagonali (4): 


h> il 02 ow 
(1) - EN >(- a G2) nee WX. 9,0 = in 
; = 2 
Le autofunzioni della (1) sono: 


(2) Pip pore Mpegs p(X, Q) = exp F P-xIT] NP) ’ 
k 


(4) L. Van Hove: Physica, 18, 145 (1952). 

(?) A. LOINGER: Nuovo Cimento, 9, 1080 (1952); 10, 356 (1953). 

(?) Vedi, ad esempio, G. MORPURGO: Suppl. al Nuovo Cimento. 8, 109 (1951). Per tutto quel che 
riguarda il metodo di Blech e Nordsieck rinviamo a questa rassegna e alla bibliografia in essa eitata 
(e, in particolare, alla memoria di PAULI e FiERZ in Nuovo Cimento, 15, 167 (1938)). 

(4) Le notazioni da noi usate sono quelle della rassegna (°). 
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ove le In (Qr) sono le consuete funzioni (normalizzate) associate ai polinomi di 
Hermite. 

L’equazione di Schrödinger del problema elettrone + campo e.m. disaccop- 
piati é: 


he 1 a | Op 
3 lanes A, + = el cr x0,1) =ih—, 
(3) | 2m, 14) >A 0Q? : a) | TA Q,0 ot 

2e 1 7 h 1/2 
(raminentiamo che X = x +h D: ORNE ul 2 er; ove e, à il versore 
k m,he|k| \e|k| 
di polarizzazione). 
Le sue autofunzioni sono: 
i 
(4) Psp ses P (x, Q) = exp h Ps II log, (Qx) x 
k 


Vogliamo ora mostrare che lo spazio dei vettori di stato (2) & ortogenale allo 
spazio dei vettori di stato (4). 

Poichè il prodotto scalare di uno qualsiasi dei vettori di stato (2) per uno qual- 
siasi dei vettori di stato (4) differisce solo in un numero finito di fattori dal prodotto 
scalare dei corrispondenti stati di vuoto (5), & sufficiente dimostrare l’ortogonalità 
di questi ultimi. Avremo pertanto da calcolare Vintegrale 


| exp i P(x +») am) IT ho(Qx) exp E = px 
i % 


il quale si puö trascrivere 


TI #10, dx TT 49» 
ke k 


| exp : (P— pi dx [] | exp [iP-Q:Yr145(Qx) Ad ; 
ma è: 


a 
II | exp LiP-Qyle (Qu) AQx = TT exp E ar aa = 
k 


E = 
eP2 [d|k| 
= exp - | =0; ev.d. 


3amehe? 5 |k| 
0 


Ovviamente, l’interesse del presente problema (come anche dei problemi studiati 
in (»2)) & essenzialmente di carattere matematico, se non altro perchè la validità 
dell’hamiltoniano della (1) & verosimilmente limitata al primo ordine in e? (vedi (?), 
p- 137). Malgrado cid, il risultato ottenuto ci sembra egualmente non privo di un 
certo interesse, in quanto esso rafforza il sospetto (insorto in seguito alle conclu- 
sioni di (»2)) che le difficoltà cui dä luogo l’applicazione dei metodi perturbativi ai 
problemi di elettrodinamica quantistica (e, più generalmente, ai problemi della 
teoria quantistica dei campi) siano dovute, almeno in parte, all’esistenza di teoremi 
di ortogonalitä analoghi a quelli validi pei casi semplici studiati. 


(8) Escludiamo gli stati stazionari a energia infinita. 


IL NUOVO CIMENTO VOTERSEN 1° Aprile 1953 


On the Non Local Interaction and on the Statistical Interpretation 
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(ricevuto il 18 Marzo 1953) 


Recently several valuable contributions to the theory of non local fields were 
made (1). The purpose of the present remarks is to show the possibility of an 
approximate description of non local interactions, starting from a local hamiltonian 
formalism and using the covariant interaction representation. The equation of 
interaction representation can be transformed in a way as to obtain the non adia- 
batic formalism proposed by Tamm (?). The stationary Schrödinger equation, used 
successfully by Tamm, Dancorr and M. Livy, is suitable to analyse the statistical 
significance of the cut-off operators obtained covariantly from a set invariant form- 
factors. 

Let us start from the well known equation of the interaction representation 
(TOMONAGA): 


(1) ii — = HY, 


where H indicates a local, invariant, hamiltonian density. H depends on free field 
operators. In the approximation here considered, we shall assume for these free 
fields the validity of the usual local equations of motion and of the usual quantum- 
commutation relations. In the case of the charge-symmetrical pseudoscalar inter- 
action between meson and nucleon fields: 


H= ig > py*t'yf’ . 
Generally, H can be an invariant density: 


H= Ÿ Ta.a, Dale) --- Pa (x), 


(!) MARKOV, YuKAWA, MCMANUS, PEIERLS, RAYSKI, ©. BLOCH. KRISTENSEN, MÖLLER, PAULI, 
(2) I. Tamm: Journ. of Phys. of Sov. Un., 9, 449 (1945); Dancorr: Phys. Rev., 78, 382 (1950). 
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where &,, are the components of different interacting fields, 73, are operators de- 
pending on spin an charge coordinates A,, and x is the four-vector (otre) | Li, tor 
example, I’ depends also or, some of the differences p’ — p" of the momenta, H re- 
presents a non local interaction. 

Introducing an invariant form factor F(wz, ....x„) depending on the differences 
between n +1 variables: 2, —x = m; .… %n—# = Mm one can define (*) a new 
hamiltonian density (*) : 


(2) H’ (x) = >» ie ee | Dy lary) + Da (Ly) (ea, «- Gaia dr — 
= > a, | Dy (x + ma) + de + NE ni Nay dre anne 


Instead of chosing one universal F, we shall assume that Here 2 Ay An) 
depends also on other parameters. We shall need to decompose covariantly every 
field operator in positive and negative frequency parts, 6.9. 


1 j 
Of (x + m) = one dép, ap) (py + m2) OF (py) exp [tprl™ + mn) | 
7 


Let, @(prärex) be the Fourier transform of Frères), which we introduce now 
in the integrands of (2). Then one has: 


(3) H'(x) = >) Reis | dp, 1 dp, exp (ip: u. Pn)® |G (Pe Aas en) OF (p1) D (pu): 


# . . 
where ae hy are new constant operators depending on the signs: &, = + of px, and 


on the annihilation or creation operators aj and: 
(4) GPs Ans Ex) = | exp [i(pyny + - + Prtin) Fre» Ars Er) Ami ++ Ain + + 


Instead of defining F, we can chose @. And the operators @ can depend expli- 
citly not only on the invariants formed with p,, but also on the parameters as 
spins, isotopic spins, and signs of the frequency of the interacting fields. If G is 
Lorentz invariant so is also F. Let us consider a term of (3), and, in chosing @. 
assume the prescription suggested previously (*): one introduces the total mo- 


mentum P, = > Pa” of the ingoing particles and the four-veloeity: u, = PVP, pe, 
oO 


then one can put: 


n 
(5) G | H G Pru”, Prußr" x 
Kal 


() G. WaracHIn: Nuovo Cimento, 8, 592 (1951); 9, 208 (1952). 
(4) See also R. Ascott: Nuovo Cimento (to be publihed), 


- [I Nuovo Cimento. 
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In the center‘ of mass system of ingoing partieles the components u, are 
(1,0, 0, 0) and @ can be chosen in away that for |p.|> ©. it tends to 0 as 
|px|~* and for |p,|=0, @=1. This last condition is necessary in order to obtain 
in the limit of slowly varying fields the usual local theory. Usually we shall have 
pan“ = m?; in some cases p® can denote a difference between two momenta. 

In some cases the vector u, in @ is the same for all elements of H’ and 
F(n;A,g,) vanishes rapidly outside the domain D;: 


nu)? — Neunk SP (er): 


The chief difficulty of the proposed method consists in the condition (Dirac, 
Biocn, Tomonaca) of commutability of H'(x) and H'(z’) in a pair of points x and a’ 
satisfying: |x —x'|>|z,—.25|. If the absolute value of the interval between these 
two events is greater than the linear dimension / of the interaction domain D,, 
this condition is fulfilled. 

We think that, although for a problem of initial values given in a pair of points 
sufficiently near to each other (e.g. given in the same domain D,), the interaction 
representation is incompatible with the introduction of formfactors, an approximate 
description with form factors of the type specified above can be valid for certain 
classes of problems, e.g. for stationary collision problems threated with the non 
adiabatic method. 

As shown recently by Cryi (5) one can obtain covariantly from the equation (1) 
the Schrodinger stationary equation of the Tamm-Dancoff method. This method has 
a great advantage of permitting to solve the statistical equilibrium problem inde- 
pendently from the question of propagation and quantization of free fields. The 
conservation laws do hold exactly in the same way as in local theories. 

From (3) one can see that every matrix-element of the hamiltonian operator 
in the Schrödinger representation is now modified by the introduction of @ factors (5). 
The introduction of such invariant factors @ for each ingoing and outgoing particle 
in every matrix element of interaction matrix H’, does not alter the hermitian 
character of H'. It is also compatible with the usual quantization of free fields 
as adopted in the T.D. formalism. Let us consider the collision of two particles 
in the center of mass system, which is supposed at rest in the stationary state con- 
sidered in the Tamm-Dancoff method. Then the G factors appear in matrices corre- 
sponding to transitions between states for whict the space component of the total 
momentum of ingoing (or outgoing) particles is = 0. They become in this case 
particularly simple and, in the case of high evergy collisions, give rise to a spherical 
simmetrical distribution ie = G(Pos m? | of outgoing particles. 

The following statistical interpretation of the «cut-off» factors G, suggested 
in earlier papers (%) seems plausible: the G factors define the statistical weight of 
the states in the C-system of ingoing particles considered in the definition of G. 
This statistical weight of a non degenerated state is — 1, in the limit of low energies 


and is — G(Po, m?) in the case: [pl AT where 7 is a universal constant. 


(5) M. Cini: A covariant formulation of the non adiabatic method. in Nuovo Cimento (to be 
published), 


(%) Nature, 142, 393 (1938); ©. R.. 207, 358. 421 (1938). 
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If in a collision there is a competition between several possible results, the 
states with lower momenta and higher multiplicity or higher masses of created 
particles acquire greater probability. 

The G factors must satisfy the condition: if p = 0, 


(6) Gm, m?) — 1. 


This condition can be considered as an equation which determines the possible 
masses of the elementary particles. The statistical factors G can depend on spin 
and charge variables. Therefore one can expect to have different mass equations 
of the type (6) for different spins and charges. Examples of such mass equations 
will be discussed in a forthcoming paper. 


LIBRI RICEVUTI E RECENSIONI 


M. ViLLA — Repertorio di Matemati- 
che, Cedam, Padova, 1951. 


E ovvio che un buon insegnante non 
deve conoscere solo le nozioni che espone 
ai suoi allievi, ma questi devono inqua- 
drarsi in una visione ampia e profonda 
di tutta la materia del suo insegnamento. 
NEE TEE 
geometria euclidea, specie nei suoi fon- 
damenti, non pud riuscire veramente effi- 
cace se l’insegnante non conosce, almeno 
nelle linee generali, le geometrie nou 
euclidee. 

Pereiö molto utile ai professori di 
matematica delle Scuole Medie (e in par- 
ticolare ai giovani laureati che intendono 
prepararsi ai concorsi) il presente vo- 
lume. che, sorto per iniziativa del pro- 
fessor VILLA, raccoglie articoli di nume- 
rosi matematici italiani opportunamente 
coordinati dallo stesso prof. VILLA. 

Alcuni articoli del libro in esame si 
riferiscono ad argomenti che nelle loro 
parti elementari rientrano nel programma 
delle nostre scuole medie. Cosi G. Rıccı 
e P. BuzANO espongono in modo com- 
pleto la teoria dei numeri; F. CONFORTO 
mette in chiara evidenza i rapporti fra 
la geometria euclidea e le geometrie non 
euclidee; G. CAMPEDELLI ed R. CALAPSO 
trattano rispettivanıente dei metodi sin- 
tetici per risolvere problemi di geometria 
e problemi risolubili con la riga e il com- 
passo; U. CASSINA espone criticamente 
la teoria dell’equivalenza e le nozioni di 
lunghezza, di area e di volumi; S. CnE- 
RUBINO considera l’applicazione dell’al- 
gebra alla geometria; A. AGOSTINI Ti- 
corda gli elementi fondamentali della 
trigonometria piana. 

Altri articoli richiamano quelle no- 
zioni di matematica svolte d’ordinario 
nei corsi nniversitari (alcune anzi rien- 


trano nel programma delle ultime classi 
del liceo scientifico) indispensabili alla 
eultura di un insegnante di matematica 
delle scuole medie. Cosi A. Amato, G. 
CIMMINO € (. RIMINI espongono rispet- 
tivamente l’analisi algebrica dei nostri 
bienni, la teoria dei limiti e degli algo- 
ritmi finiti, il calcolo differenziale e inte- 
grale. M. VıLta richiama invece la geo- 
metria analitica e la teoria generale delle 
trasformazioni geometriche; L. CAMPE- 
DELLI la teoria delle curve € superfici; 
P. Buzano e C. Rımını la teoria dei vet- 
tori e dei calcoli approssimati; G. Cım- 
MINO gli elementi della teoria delle fun- 
zioni analitiche; G. ZAPPA i gruppi e le 
loro applicazioni alle equazioni algebriche. 
Importante l’articolo di F. SIBIRANI sulla 
probabilita e sulla matematica finan- 
ziaria, argomenti che rientrano nei pro- 
grammi di matematica degli Istituti tec- 
nici ma che non sono svolti in nessun 
corso fondamentale per la laurea in 
Matematica o in Matematica e Fisica. 
Chiude il libro un articolo di A. AGOSTINI 
sulla storia delle matematiche. Manca la 
trattazione delle equazioni differenziali 
perché non figurano nei programmi d 
concorso a cattedre di matematica nelle 
scuole medie. Perd gli estensori di quei 
programmi hanno dimenticato che il pro- 
fessore di matematica nelle scuole medie 
superiori insegna anche la fisica e che 
per una conoscenza che non sia del tutto 
superticiale di questa disciplina, è ne- 
cessaria la conoscenza delle equazioni 
differenziali. Ma non & qui il luogo di 
insistere su questa ed altre assurdita dei 
programmi di concorso (1). 


(*) Mi limito ad osservare che sarebbe ora 
di rifare il programma di fisica per i concorsi 
nelle scuole medie, programma veramente anti- 
diluviano, 
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Comunque & doveroso dar lode al 
VILLA ed ai suoi collaboratori che hanno 
arriechito la matematica italiana di un 
nuovo importante trattato. Naturalmente 
& spontaneo domandarsi perché non si fa 
altrettanto per la fisica, cioè perché non 
si prepara un libro che permetta di ap- 
profondire la fisica in relazione al suo 
insegnamento elementare (?). F ben vero 
che la redazione di questo libro offre 
notevoli difficoltà sia per la maggior com- 
plessitä della fisica, sia per il suo rapido 
sviluppo, sia per la mancanza di un la- 
voro critico analogo a quello che 1 mate- 
matici hanno fatto per le parti elemen- 
tari della loro scienza: sono certo pero 
chei fisici italiani, volendo, possono supe- 
rare questi ostacoli. Ritengo che ne 
avrebbe vantaggio non soltanto linse- 
enamento della fisica nelle scuole medie. 


DARIO GRAFFI 


(4) Aleuni anni fa, per iniziativa del com- 
pianto prof. GIOVANNI GENTILE jr. si cominciö 
a redigere un trattato di fisica del tipo accen- 
nato nel testo. Di questo trattato & uscito perd 
solo il 1° volume e non si ha notizia se verrà 0 
no completato. 


A. DAUVILLIER — La physique cos- 
mique. Bibliothéque de Philoso- 
sophie scientifique, pag. 846. Flam- 
marion, Parigi, 1951. 


La fisica cosmica è una scienza nuova 
di coordinamento e di sintesi che costi- 
tuisce, oggi, l’appassionato oggetto di 
studio di molti scienziati che sentono il 
problema di una sistemazione logica e 
coerente di tutte le acquisizioni teoriche 
e sperimentali in ogni campo dello scibile. 

La fisica cosmica noné geofisica né 
astrofisica, è invece una scienza generale 
di universale portata che unisce, inquadra 
e comprende non solo quelle discipline, 
ma anche la fisica propriamente detta, 
la chimica cosmica e la stessa biologia. 

In questo volume, della collana di 
«Filosofia Scientifica», Edizione Flam- 
marion, l’autore che ha buona esperienza 


di lavoro nel campo della fisica, della geo- 
fisica ed astrofsica, presenta in forma di- 
scorsiva ilrisultato di suoi personali pensa- 
menti e valutazioni e si sforza di dare una 
sistemazione di portata cosmica a leggi, 
fatti e fenomeni propri della fisica, della 
cosmologia, della geologia e della biologia. 

Cosi egli affronta il problema del- 
Vorigine del sistema solare e in parti- 
eclare del nostro pianeta, e propone una 
nuova ipotesi cosmogonica basata sulla 
generazione di coppie di pianeti gemelli, 
che fornirebbe la spiegazione di molte 
particolarità del sistema solare. 

Altri importanti argomenti, come l’ori- 
gine del rilievo terrestre, l’orogenesi, la 
glaciazione quaternaria ed il meccanismo 
delle glaciazioni, vengono dall’autore di- 
scussi e trattati con personali deduzioni. 

Vengono poi i tre capitoli concernenti 
il magnetismo dei corpi celesti, il geo- 
magnetismo e quello su i raggi cosmici; 
in essi l’autore espone più © meno suc- 
cintamente i fatti cercando di coordinarli 
e tentandone una interpretazione, per 
quanto espressamente ed onestamente 
dichiari che la frammentarieta e l’incom- 
pletezza delle nostre conoscenze non puod 
in questo campo certamente condurre a 
soluzioni definitive. 

In due successivi capitoli sulla « Ge- 
nesi della vita» e sul « Problema della 
evoluzione » l’autore espone alcune sue 
considerazioni sulle cause geochimiche 
e cosmiche che fornirebbero una spiega- 
zione della formazione dei continenti e 
dell’apparizione della vita sul pianeta. 
In particolare l’apparizione della vita 
sulla Terra sarebbe il risultato di una 
concomitanza straordinariamente fortuita 
di condizioni fisico-chimiche favorevoli, 
e il fattore determinante dell’evoluzione 
sarebbe poi lo psichismo. 

Il libro finisce con un capitolo su 
«Il valore sociale della Scienza » che pra- 
ticamente si ridurrebbe a fornire un 
passatempo alle persone più intelligenti 
ed a risolvere il problema dell’alimenta- 
zione della schiatta umana ! 


M. GIORGI 
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H.GARTMANN - Raumfahrt-Forschung, 
R. Oldenbourg Ver ag, München, 
1952, in-8°, pp. 200, fig. 57. 


L’ing. GARTMANN ha raccolto in que- 
sto volume una serie di saggi sul pro- 
blema del volo interplanetario aventi per 
autori alcuni dei pionieri tedeschi di 
questo campo di ricerche: dal prof. HER- 
MANN OBERTH, che si occupa questa 
volta delle « stazioni spaziali », a WILLY 
Ley, del quale sono noti anche in Italia 
gli seritti di carattere divulgativo sui 
razzi, e che qui racconta brevemente la 
storia delle fantasie, dei progetti e degli 
esperimenti relativi ale possibilita di eva- 
sione dalla Terra. 

Dopo che lo stesso CLARKE, organiz- 
zatore delle attivita della British Inter- 
planetary Society, si & dedicato alla 
compilazione di romanzi del genere 
«science fiction » ora di moda, i dubbi 
sulla serieta di alcune tra queste specu- 
lazioni sono cresciuti notevolmente. Bi- 
sogna riconoscere, perd, che i membri 
degli analoghi gruppi tedeschi, affiliati 
di quella Federazione Astronautica In- 
ternazionale cui & dedicato ultimo capi- 
tolo del presente volume, manifestano 
un grado di preparazione tecnica tale 
da tar apparire bene sviluppata, almeno, 
Vimpostazione teorica di questi problemi 
dal punto di vista della fisica. 

Particolarmente interessanti, infatti, 
risultano le applicazioni analitiche della 
meccanica celeste, al problema della na- 
vigazione interplanetaria, presentate dal 
prof. WERNER SCHAUB, e lo studio della 
fase preliminare costituita da una « son- 
da» percorrente un’orbita intorno alla 
Terra, sviluppato dagli ingegneri ENGEL 
e HanıscH tenendo conto delle possibi- 
lita attuali di un razzo a sei stadi da 
220 tonnellate. 

Meno sicure, invece, le considerazioni 
del dott. DIRINGSHOFEN sul comporta- 
mento del corpo umano in assenza di 
gravità, questione su cui solo ora si è 
cominciato a raccogliere dati sperimen- 
tali mediante il lancio di scimmie in 


razzi per alta quota, 
V. SOMENZI 
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K. E. ZIMEN - Angewandte Radio- 
aktivität, vim-+124 seiten mit 45 
Textabbildungen und einer Ta- 
fel, Springer, Berlin 1952. 


Questo volumetto di radioattivita 
«applicata » dovuto alla penna di uno 
specialista e presentato da O. Hany, rac- 
coglie in poco pit che un centinaio di 
pagine cid che & strettamente necessario 
conoscere per poter lavorare con isotopi 
radioattivi. 

La prima parte, di 56 pagine, con- 
tiene in un agile e scarno riassunto i fon- 
damenti scientifici dei metodi, descritti 
poi per mezzo di esempi nella seconda 
parte (28 pagine). Le ultime quaranta 
pagine sono dedicate a tabelle di dati 
d’uso corrente e ad una aggiornata e 
accurata bibliografia. 

Un singolare pregio di questa operetta 
consiste nel non perdere mai di vista lo 
scopo, che & quello di costituire una utile 
fonte di informazioni per chi si serve 
della radioattivita come strumento di la- 
voro senza preoccuparsi degli aspetti e 
delle questioni fondamentali ch’essa pre- 
senta. Di conseguenza, i resultati dell’in- 
dagine scientifica sono presentati in modo 
rapido e dogmatico, sempre dal punto di 
vista sperimentale, facendo posto alla 
formula empirica che risponde alle esi- 
genze pratiche assai più che alla descri- 
zione dell’esperienza concettuale che met- 
te in luce una difficoltà di principio o 
all’impostazione teorica del calcolo esatto 
di una determinata grandezza. 

Naturalniente, nessuno potra preten- 
dere di trovare qui tutti i dati che gli 
servono per la progettazione e l’esecu- 
zione delle sue esperienze; ma il chimico, 
il biologo, il medico, il tecnico e sotto 
certi aspetti anche il fisico, vi troveranno 
una guida chiara, e generalmente det 
tutto attendibile, per risolvere dubbi, 
chiarire possibilita e limiti dei vari me- 
todi, progettare esperienze e misure in 
modo sensato. 

Non manca tuttavia, qua e là, qualche 
frase che lascia perplessi, come alla 
pagina 16 dove si afferma, che la causa 
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(Ursache) di una disintegrazione radio- 
attiva è costituita dal fatto che l’energia 
totale della particella emessa e del nucleo 
residuo à minore dell’energia del nucleo 
di partenza. E neppure nıanca qualche 
errore, come alla pag. 28, ove dalla 
fig. 14-3 si ricava che la velocità di una 
particella « nell’aria resta costante fino 
a un paio di centimetri dalla fine del 
percorso; o alla pag. 53, ove a proposito 
della reazione C(x, p)N®, endotermica 
per 5 MeV, si afferma che per ottenerla 
sono necessarie particelle « di almeno 
5 MeV, senza tener conto che la soglia è 
pit elevata, dovendosi nel processo con- 
servare la quantità di moto totale. 

A parte questo, il libro indubhia- 
mente à assai utile e consigliabile a chiun- 
que si appresti a lavorare con radio- 
isotopi. 


M. AGENO 


R. ROTHE — Höhere Mathematik 
Teubner, Leipzig 1952; 4 voll. per 
cirea page. 1000 complessive; 
in-160. 


La nuova edizione, fatta a cura di 
W. ScHMEIDLER, di quest’opera gia ben 
nota, che trae origine dai corsi propedeu- 
tiei tennti dal Prof. Rotue alla Technische 
Hochschule di Berlino, differisce dalla 
precedente del 1949 per aleuni migliora- 
menti e poche aggiunte, restando inalte- 
rata nella sostanza. 

Sull’Opera stessa — senz’altro consi- 
gliabile come utile e pratico testo di con- 
sultazione per ingegneri e per chiunque 
applichi le scienze matematiche — ci limi- 
tiamo percid a riportare, oltre all’indica- 
zione degli argomenti trattati e ad una 
breve osservazione critica, le seguenti pa- 
role della prefazione che chiaramente illu- 
strano pregi e limiti del lavoro: , 

«... si è data molta importanza alle 
relazioni tra la matematica pura e le sue 
applicazioni e ai vari metodi della mate- 
matica applicata, senza perd trascurare le 
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cose più importanti della matematica 
pura stessa ». 

«Cresciuta dagli appunti per studenti 
del Prof. RoTHE, l’Opera costituisce una 
guida molto breve, con disposizione libera 
della materia e con molti esempi, appli- 
cazioni, esercizi. La brevita dello svolgi- 


_ mento non facilita lo studio. anzi richiede 


una collaborazione viva e coraggiosa che 
certamente non pud essere effettuata da 
tutti senza sospiri e sudore. Cid yale 
soprattutto per gli esercizi, per i quali 
sono dati solo una guida molto breve e 1 
risultati ». 

Nonostante questo indirizzo essenzial- 
mente pratico dell’Opera -- o forse prc- 
prio per questo — sarebbe talora deside- 
rabile maggior rizore critico. 

Per esempio, proprio in un’aggiunta 
di questa edizione (vol. I, pag. 88), al 
n. 11 del § 13, intitolato « Integrations- 
grenzen bei neuen Veränderlichen», si dice: 

«Sei z=qg(t) und umgekehrt t=y(), 
so wird die untere Grenze... t, = y(z,) = 
= y(a), entsprechend die obere 1, = 
= ylx,) = y(b); also, ganz allgemein 


b 
F= | f(x) dr = 


4, ty = po) 


= | jw) dx || flpb|p'(t) dt.» 


ay t= pla) 


Ora, se si bada alla definizione di fun- 
zione inversa data a pag. 15 (+), ed al 
fatto che per la sua esistenza non si pon- 
gono condizioni, si ottengono sofismi di 
tipo ormai elassico, come il seguente: 


(1) «Wenn in y = f(x), æ mit y vertauscht 
wird, æ = f(y), so bildet die Gesamtheit der 
Werte, die dadurch jedem Werte von x eines 
Bereiches zugeordnet werden können, die 
Umkehrungsfunktion y = g(x)”. 
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ponendo: # hit: 


E bensi vero che prima, a pag, 45, 
è dato il problema: 

«Ist y = f(z) für a <x <b eindeutig, 
stetig und beständig zunehmend (oder 
beständig abnehmend), ohne strecken- 
weise konstant zu sein, so hat die aus 
x — f(y) sich ergebende Umkehrungs- 
funktion y = g(x) in dem entsprechenden 
Bereiche dies2lben Eigenschaften ». 

Perö tale teorema non serve se non 
si avverte che in seguito si considere- 
ranno funzioni inverse solo in tali ipotesi. 
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La materia contenuta nei primi tre 
volumi € — grosso modo, e astrazion fatta 
dal molto diverso rapporto tra teoria e 
pratica — quella che generalmente viene 
svolta nei due corsi d’Analisi Matematica 
delle nostre Facolta di Scienze; ove perd 
si aggiungano numerose applicazioni alla 
geometria, alla meccanica e in generale 
alla tisica. 

Il quarto volume, costituito da tre 
quaderni staccati, & una varia e ricca 
raccolta d’esercizi. 

Non ci risulta che sia stata ristampata 
la Formelnsammlung zur höheren Mathe- 
matik dello stesso Autore, che costituisce 
un utile complemento dell’Opera. 


D. Fürst 
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